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Preface 


This  thesis  is  the  result  of  a  suggestion  by  Dr. 

Corrado  R.  Poli,  of  the  Department  of  Aeronautical  Engineer¬ 
ing  at  AFIT,  for  a  study  on  gravity-assisted  solar  probe 
trajectories.  The  need  for  this  investigation  was  estab¬ 
lished  at  the  NASA-Stanford  University  Summer  Training  Pro¬ 
gram  in  Systems  Engineering  which  he  attended  in  the  summer 
of  1966.  Results  of  this  program,  a  systems  analysis  c£  a 
solar  probe  mission,  are  reported  in  the  ICARVS  design  study 
(Ref  12). 

The  thesis  is  written  on  a  level  which  should  provide 
no  difficulty  to  a  reader  who  has  completed  a  first  course 
in  astrodynamics.  The  method  developed  for  gravity-assisted 
trajectories  is  useful  in  feasibility  studies,  and  can  be 
applied  to  any  ballistic  flyby  mission.  Gravity  assist  has 
been  found  to  be  useful  (Ref  17)  in  many  types  of  missions 
such  as:  deep-space  probes,  trips  to  the  outer  planets, 
grand  tours  or  trajectories  which  pass  several  planets  on  a 
single'  flight,  and  out-of-ecliptic  trajectories.  If  the 
reader  finds  an  interest  in  this  area,  I  can  suggest  two 
thesis  topics  which  would  supplement  this  work*  1)  applica¬ 
tion  of  the  gravity-assist  equations  to  continuous-thrust 
missions  and  2)  a  linear  error  analysis  of  gravity-assisted 
trajectories. 

A  discussion  on  computer  programming  is  included  with 
each  mission  study  in  the  thesis,  and  the  programs  with 
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sample  output  are  given  in  Appendix  E.  These  sections  are 
provided  for  the  reader  who  may  wish  to  continue  the  study, 
and  if  desired,  may  be  skipped  with  no  discontinuity  in  the 
subject  material, 

I  wish  to  express  my  appreciation  to  Dr,  Poli  for  his 
suggestions  and  continued  interest  in  my  work,  and  to  the 
computational  staff  at  Building  57,  Wright-Patterson  Air 
Force  Base,  for  their  assistance  in  programming  the  equations. 
Finally,  I  thank  my  wife  for  the  intangible  contributions  she 
ha3  made  to  this  study,  and  for  the  typing  of  the  manuscript. 

Kenneth  A.  Myers 
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Abstract 

An  investigation  of  direct-transfer  trajectories  for 
solar  probe  missions  indicates  that  expensive  launch  vehi¬ 
cles  such  as  the  $125-million  Saturn  V  are  required  to 
achieve  perihelia  less  than  0.25  AU  and  inclination  angles 
above  20°.  Methods  are  developed  for  two-  and  three-dimen¬ 
sional  gravity-assisted  trajectories  (trajectories  which 
pass  through  one  or  more  planetary  gravitational  fields) , 
and  are  applied  to  solar  probe  missions  in  an  attempt  to  re¬ 
duce  launch  vehicle  costs.  The  analysis  is  based  on  the 
pieced-conic  approximation  and  the  assumption  of  circular 
coplanar  planetary  orbits.  It  is  found  that  a  perihelion-  of 
0.16  AU  can  be  obtained  with  the  $16-million  Atlas/Centaur/ 
TE-364-3  by  using  a  Venus  assist.  Even  greater  reductions 
in  perihelia  are  attained  with  multiple  passes  at  Venus;  in 
addition,  these  missions  allow  the  exploration  of  several 
regions  near  the  Sun  with  a  single  launch.  Venus-Mercury 
combination-assisted  trajectories  are  of  little  value  in 
solar  probe  missions.  A  Venus  assist  is  useful  for  10°-out- 
of-ecliptic  trajectories  when  the  pre-assist  orbit  lies  in 
the  ecliptic.  For  a  pre-assist  orbit  out  of  the  ecliptic, 
an  inclination  angle  of  25°  and  a  perihelion  of  0.53  AU  can 
be  achieved  with  the  Atlas/Centaur/TE-364-3  launch  vehicle. 
Solar  impact  and  90°-out-of-ecliptic  trajectories  can  be 
attained  with  a  Jupiter  assist;  however,  traversal  of  the 
asteroid  belt  and  a  three-year  mission  time  decreases  space¬ 
craft  reliability. 


xv 


I 


GA/AE/67-4 


GRAVITY-ASSISTED  TRAJECTORIES  FOR 
SOLAR  PROBE  MISSIONS 

I,  Introduction 

aw  — »»—»—— mmm 

Background  Information 

The  current  NASA  Pioneer  program  (Ref  22:69)  is  designed 
to  explore  solar  effects  on  the  space  environment  between  0.8 
and  1,2  AU  from  the  Sun.  Plans  for  advanced  missions  in  this 
program  include  the  exploration  of  regions  as  close  as  0.5  AU 
from  the  Sun,  both  in  and  out  of  the  ecliptic  plane.  The  tra¬ 
jectories  for  these  missions  can  be  achieved  with  the  Thrust- 
Augmented  Delta  and  Atlas  launch  vehicles,  but  larger  and 
more  expansive  booster  systems  are  required  for  the  explora¬ 
tion  of  regions  closer  to  the  Sun. 

The  perihelion  which  can  be  attained  for  a  given  burnout 
velocity  is  shown  on  the  graph  in  Fig.  1  for  a  hyperbolic 
ascent  departure  and  Hohmann  transfer  orbit.  Burnout  veloc¬ 
ities  in  excess  of  41,000  ft/sec  are  necessary  for  solar 
probe  missions  to  regions  less  than  0.5  AU  from  the  Sun.  The 
reason  for  this  high  velocity  requirement  is  that  the  space¬ 
craft  velocity  must  not  only  exceed  the  velocity  for  Earth 
escape  (36,700  ft/sec),  but  also  counteract  the  high  orbital 
velocity  of  the  Earth  about  the  Sun  (97,700  ft/sec). 

For  missions  to  the  0,2  to  0.1  AU  region,  the  range  of 
required  burnout  velocities  is  from  55,000  to  67,000  ft/sec, 
respectively.  This  requirement  can  be  met  with  the  Saturn  V 
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launch  vehicle,  but  the  cost  of  this  system  ($125  million) 
and  the  small  payload  for  a  solar  probe  (200  lb)  warrant  an 
investigation  for  less  expensive  techniques. 

Purpose  of  the  Thesis 

The  purpose  of  this  thesis  is  to  investigate  the  use  of 
gravity  assist  as  a  means  of  reducing  the  cost  (see  Table  I) 
and/or  burnout  velocity  requirements  for  booster  vehicles  in 
solar  probe  missions.  Gravity  assist  is  defined  as  a  signif¬ 
icant  trajectory  perturbation  between  launch  and  target 
positions  caused  by  a  close  approach  to  an  intermediate 
planet.  This  concept  is  also  referred  to  as  a  planet  fly-by 
or  swing-by.  The  basic  purpose  of  gravity  assist  is  to  change 
the  spacecraft  velocity  vector  with  respect  to  the  Sun;  the 
effectiveness  of  the  maneuver  depends  primarily  on  the  mass 
and  velocity  of  the  assist  planet  and  on  the  DOCA  (distance 
of  closest  approach)  of  the  spacecraft  to  the  surface  of  the 
planet. 

The  primary  goal  is  to  achieve  the  smallest  possible 
perihelion  with  the  least  expensive  booster  vehicles;  for 
out-of-ecliptic  missions  it  is  also  desired  to  achieve  the 
largest  possible  inclination  angle. 

History  of  Fly-By  Missions 

The  first  planetary  fly-by  was  achieved  in  December,  1962, 
when  the  Mariner  2  spacecraft  (Ref  4:695-699)  passed  within 
21,600  mi  (18,800  NM)  from  the  surface  of  Venus.  The  purpose 
of  this  flight  was.  to  collect  data  on  the  Venusian  atmosphere, 
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so  the  mission  was  completed  ^fier  the  encounter,  and  the 
spacecraft  was  left  to  orbit  the  Sun.  Mariner  4  (Ref  24* 

Fall  1965)  completed  a  similar  mission  when  it  passed  within 
5400  mi  (4690  MM)  from  the  surface  of  Mars  in  July,  1965. 

Survey  of  the  Literature 

Minovitch  (Ref  10,11)  has  conducted  a  two-part  investi¬ 
gation  on  gravity-assisted  trajectories.  The  first  part  of 
his  investigation  is  devoted  to  the  development  of  theory 
and  the  calculation  of  gravity-assisted  trajectories  to  the 
inner  planets.  The  second  part  is  devoted  to  deep-space  and 
solar  probe  trajectories  in  which  Jupiter  is  used  as  an  assist 
planet  for  missions  in  and  out  of  the  ecliptic  plane. 

Niehoff  (Ref  13)  has  continued  the  work  done  by  Minovitch  for 
gravity-assist  missions  in  the  ecliptic  plane  by  using  Jupiter 
as  either  a  target  or  an  assist  planet.  The  results  of  these 
studies  show  that  90°-out-of-ecliptic  trajectories  with  per¬ 
ihelia  less  than  0.1  AU  can  be  attained  with  a  Jupiter  assist 
and  an  Atlas/Centaur/Kick  launch  vehicle  ($16  million). 

Thus,  Jupiter-assisted  solar  probe  trajectories  have 
been  thoroughly  investigated.  The  excellent  performance 
attained  in  these  miss  .s  is  due  to  the  large  mass  of  Jupiter 
(318  times  the  mass  of  Earth) ,  but  these  missions  have  two 
disadvantages:  first,  the  flight  time  is  about  three  years, 
and  second,  the  vehicle  must  pass  through  the  asteroid  belt 
twice  on  its  path  to  the  Sun. 

Ross  (Ref  15)  has  presented  a  detailed  study  on  round 


4 


ESTIMATED  COST  OF  LAUNCH  VEHICLES  (1967) 


GA/AE/67-4 


00 

m 

in 

00 

•• 

•• 

H 

n 

CM 

CM 

cm 

CM 

H 

IM 

4-i 

»W 

<1) 

<1} 

0) 

to 

e 

& 

e 

0 

0 

0 

u 

u 

u 

to 

to 

to 

' 

* 

■K 

« 


5 


GA/AE/67-4 


trips  to  the  inner  planets  and  an  introduction  to  Venus- 
assisted  trajectories  for  solar  probes.  The  Jupiter-assist 
disadvantages  are  avoided  by  using  the  inner  planets  for 
gravity  assist  in  solar  probe  missions. 

Plan  of  Development 

The  over-all  plan  of  the  thesis  is  to  investigate  inner- 
planet  gravity-assist  missions, and  to  compare  them  with 
direct-transfer  and  Jupiter-assist  missions.  The  body  of 
the  report  contains  a  discussion  of  the  theory  and  results 
of  the  mission  studies.  Detailed  derivations  of  the  equa¬ 
tions  are  given  in  the  appendices. 

Chapter  II  is  an  investigation  of  direct-transfer  solar 
probe  trajectories.  This  chapter  serves  to  illustrate  the 
departure  orbit  technique  used  in  the  report,  and  it  provides 
a  standard  with  which  gravity-assisted  solar  probe  trajec¬ 
tories  may  be  compared. 

The  method  of  analysis  for  two-dimensional  gravity- 
assisted  trajectories  is  discussed  in  Chapter  III.  This 
method  is  applied  to  the  following  solar  probe  missions: 

1,  Venus  Assist  in  the  Ecliptic  (Chapter  IV) 

2,  Venua-Mercury  Combination  Assist  (Chapter  V) 

3,  Multiple  Venus  Assist  (Chapter  VI) 

a.  Double  Pass  at  Venus 

b.  Triple  Pass  at  Venus 

The  method  in  Chapter  III  is  extended  to  the  case  of 
three-dimensional  gravity  assist  in  Chapter  VII,  and  is 
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applied  to  the  following  out-of-ecliptic  solar  probe 
missions: 

1,  Out-of-ecliptic  Venus  Assist  (Chapter  VIII) 

a.  Pre-Assist  Orbit  in  the  Ecliptic  Plane 

b.  Pre-Assist  Orbit  out  of  the  Ecliptic  Plane 

2.  Jupiter  Assist  (Chapter  IX) 

It  has  been  shown  (Ref  11:39)  that  Mars  offers  little 
or  no  benefit  as  an  assist  planet  for  .olar  probe  trajec¬ 
tories,  so  these  missions  have  not  been  investigated  in 
this  study. 

Conclusions  are  made  in  Chapter  X,  and  favorable 
missions  are  recoir  ended  for  further  consideration. 

The  derivation  of  the  equations  for  gravity-assisted 
trajectories  is  included  in  Appendix  A,  B,  and  C,  and  a 
method  for  the  generation  of  launch  opportunities  is 
described  in  Appendix  D.  Computer  programs  used  in  the 
numerical  calculations  are  given  in  Appendix  E. 

Mathematical  'todel 

The  assumptions  and  limitations  used  in  this  study 
are  as  follows: 

1,  Planetary  orbits  are  assumed  to  be  circular  and 
coplanar,  except  for  Mercury  which  is  assumed  to  have  a 
coplanar  orbit  with  an  eccentricity  of  0.2. 

2.  The  pieced-conic  approximation  is  used,  so  that 
the  calculation  of  each  gravity-assisted  trajectory  is  a 
sequence  of  two-body  problems. 
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3.  1_  '  assumed  that  a  planet  and  its  SOI  (sphere 

of  influence)  may  be  approximated  as  a  point  mass  on  the 
planetary  orbit  when  viewed  on  a  heliocentric  scale. 

Thus,  a  given  pre-assist  trajectory  (the  trajectory  from 
Earth  departure  to  the  SOI  of  the  assist  planet)  may  be 
assumed  to  intersect  the  SOI  at  any  desired  point.  The 
actual  point  of  intersection  may  be  controlled  by  mid¬ 
course  guidance  maneuvers  or  by  small  variations  in  the 
launch  time.  The  validity  of  this  assumption  is  analyzed 
in  Chapter  III. 

4.  The  length  of  time  which  the  spacecraft  spends  in 
the  SOI  of  the  assist  planet  is  assumed  to  be  negligible 
compared  to  the  length  of  the  total  mission  time.  Thus, 
the  spacecraft  position  vector  with  respect  to  the  Sun  can 
be  approximated  by  the  planetary  position  vector  at  the 
time  of  closest  approach  to  the  assist  planet. 

5.  Ballistic  trajectories  with  a  single  impulse  at 
Earth  launch  are  employed  in  all  missions.  The  gravity- 
assist  techniques  which  are  developed  can  also  be  applied 
to  continuous-thrust  missions;  the  idea  of  applying  an 
impulse  during  planetary  encounter  has  been  investigated 
by  Ross  (Ref  15jl49), 

6.  The  departure  orbit  is  a  direct-ascent  hyperbola 
from  Earth;  thus,  the  HEV  (hyperbolic  excess  velocity)  for 
the  departure  orbit  is  a  function  only  of  a  specified 
burnout  velocity.  The  determination  of  azimuth,  flight 
path  angle,  right  ascension,  and  declination  at  burnout 
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is  left  for  more  detailed  investigations  (Refs  5,  9). 

7.  The  trajectories  are  designed  so  that  perihelion 
lies  inside  the  Earth’s  orbit  about  the  Sun.  That  is,  the 
heliocentric  flight  path  angle  (angle  between  spacecraft 
velocity  vector  and  inner  radial  line  from  Sun  to  Earth)  at 
the  time  of  Earth  launch  is  always  equal  to  or  less  than  90°, 
Modification  of  the  equations  for  trajectories  outside  the 
Earth's  orbit  is  illustrated  in  the  Jupiter-assist  missions 
of  Chapter  IX. 

8,  Only  Type  I  trajectories  (heliocentric  transfer 
angle  from  Earth  to  the  assist  planet  612  ^  180°)  are 
investigated  in  this  report.  Additional  circuits  about 
the  Sun  before  interception  of  the  assist  planet  3erve 
only  to  increase  the  time  and  reliability  requirements  of 
the  spacecraft. 

General  Method 

In  general  the  thesis  is  a  collection  of  preliminary 
mission  studies.  Trajectories  are  calculated  for  each 
mission  by  selecting  four  arbitrary  parameters:  a  burnout 
velocity  V^,  departure  flight  path  angle  <j>0,  inclination 
angle  i,  and  distance  of  closest  approach  d. 

For  convenience,  the  equations  derived  in  the  appen¬ 
dices  are  programmed  on  the  IBM  7094  digital  computer,  and 
trajectories  are  generated  for  the  desired  range  of  V^,  <|>0, 
i,  and  d.  Each  chapter  includes  a  discussion  of  the  program 
for  the  particular  mission  under  consideration,  and  results 
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are  presented  in  graphical  form. 

A  list  of  solar  system  data  used  for  the  computations 
is  given  in  Table  II.  Gravitational  parameters  and  the 
astronomical  unit  are  taken  from  Reference  8,  Chapter  III, 
pp.  2  and  3;  equatorial  radii  are  from  Reference  16,  p,  371; 
all  other  data  are  from  Reference  2,  p.  378,  The  perihelion 
and  aphelion  for  Mercury,  0.30750  and  0.46670  AU,  respec¬ 
tively,  are  taken  from  Reference  23,  p.  14. 
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II.  Direct-Transfer  Trajectories  for  Solar  Probe  Missions 
Purpose  and  Scope 

The  purpose  of  this  chapter  is  threefold:  first,  to 
investigate  the  velocity  requirements  for  solar  probe 
missions;  second,  to  introduce  the  equations  which  are  used 
for  Earth  departure  orbits;  third,  to  develop  a  set  of 
reference  trajectories  which  may  be  compared  with  gravity- 
assisted  trajectories  in  later  chapters. 

The  trajectories  for  direct-transfer  missions  are  bro¬ 
ken  into  two  phases: 

1,  Geocentric  escape  orbit  -  direct-ascent  hyperbolic 
orbit  from  the  surface  of  the  Earth  to  the  boundary  of  its 
SOI. 

2.  Heliocentric  transfer  orbit  -  an  elliptical  orbit 
which  originates  at  a  point  on  the  Earth’s  orbit  about  the 
Sun. 

Geocentric  Escape  Orbit 

The  initial  phase  of  a  solar  probe  mission  involves 
escape  from  the  Earth's  gravitational  field.  A  direct- 
ascent  hyperbola  is  employed  in  all  the  trajectories;  it  is 
assumed  here  that  parking  orbits,  used  primarily  for  storage 
dumps  and  pre-departure  checkout,  are  seldom  used  in  solar 
probe  missions. 

Velocity  on  a  Direct-Ascent  Hyperbola.  From  two-body 
theory,  it  is  known  that  parabolic  velocity  is  the  limiting 
velocity  for  escape  from  the  Earth's  gravitational  field. 
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However,  hyperbolic  velocity  is  usually  employed  in  order 
to  avoid  the  long  time  interval  required  for  escape  on  a 
parabolic  orbit.  Also,  the  hyperbolic  excess  velocity 
(HEV)  for  a  parabolic  orbit  is  zero,  so  it  is  necessary 
that  the  escape  orbit  be  hyperbolic  if  it  is  desired  to 
achieve  a  perihelion  less  than  1.0  AU. 

The  velocity  of  a  spacecraft  on  a  hyperbolic  orbit 
is  given  by  (Ref  16:96) 


V*  = 


(2-1) 


where  r  is  distance  between  the  spacecraft  position  and 
the  center  of  mass  of  the  planet  P1#  When  the  spacecraft 
approaches  the  SOI  of  Pj,  ,  r  +  »,  and  the  velocity  becomes 


(2-2) 


which  is  the  HEV.  Thus,  for  any  point  on  the  hyperbolic 
orbit, 


V2 


(2-3) 


Selection  of  V^.  The  vehicle  is  launched  in  an 

easterly  direction,  as  shown  in  Fig.  2,  to  take  advantage 
of  the  Earth's  rotational  velocity.  Burnout  usually  occurs 
from  50  to  100  mi  above  the  surface,  but  the  radius  of  the 
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Earth  (Bp1~  3960  mi)  serves  as  a  good  approximation  for  the 
actual  burnout  radius.  Thus, 


H  =  V-Pl+  (2-4) 

RPl 

The  use  of  Eq  (2-4)  avoids  the  sometimes  conflicting 
definitions  of  ideal,  mission,  and  total  velocities.  Vb  is 
the  actual  velocity  of  the  vehicle  with  respect  to  the  Earth 
when  burnout  occurs.  The  velocity  required  to  accomplish 
the  mission  is  greater  than  Vjj  because  of  gravity  and  drag 
losses,  engine  inefficiencies,  and  the  Earth's  rotational 
velocity  (Ref  19*136-138.) . 

The  values  of  Vb  used  in  Eq  (2-4)  may  be  compared 
directly  with  data  for  existing  launch  vehicles  in  Pig.  3. 
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PAYLOAD -POUNDS 


Fig.  3 

Launch  Vehicle  Performance 

(From  Ref  25:Chap.  II,  p.  46) 

In  turn,  the  cost  of  these  vehicles  may  be  found  in  Table  I. 
An  expanded  view  of  the  region  of  interest  (solar  probe  pay- 
loads  less  than  500  lb)  is  given  on  the  chart  in  Fig.  4. 

Selection  of  fixes  the  value  of  the  HEV,  Voop^,  (see 
Fig.  5)  which  is  then  used  in  calculating  the  heliocentric 
transfer  orbit. 

Heliocentric  Transfer  Orbit 

After  the  escape  phase  of  the  mission  is  completed,  the 
vehicle  enters  a  heliocentric  elliptic  orbit.  The  actual 
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departure  trajectory,  shown  as  a  heavy  line  in  Fig.  6,  is 


a  combination  of  the  motion  of  the  vehicle  on  its  escape 
hyperbola  and  the  motion  of  the  Earth  about  the  Sun.  The 
thin  lines  represent  the  hyperbolic  escape  orbit  with 
respect  to  the  Earth  for  a  sequence  of  observations  after 
launch. 

Selection  of  <J>0,  i.  As  stated  in  the  introduction, 

the  trajectories  investigated  in  this  report  depend  upon 
the  choice  of  a  burnout  velocity  Vb,  departure  flight  path 
angle  <j>  0 ,  and  an  inclination  angle  i .  The  parameter 
specifies  a  departure  orbit  which,  in  turn,  fixes  the 
magnitude  of  The  selection  of  <f>0  and  i  defines  the 

geometry  of  a  heliocentric  orbit  which  is  assumed  to  orig~ 
inate  from  a  point  on  the  Earth's  orbit  about  the  Sun. 

Calculation  of  70.  The  vector  geometry,  shown  in  Fig. 
7 ,  is  used  to  calculate  W0 ,  the  velocity  of  the  spacecraft 

with  respect  to  the  Sun  after  Earth  escape.  The  XjYjZj 
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Fig.  8 

Heliocentric  Transfer  Orbit 
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coordinate  system  is  located  with  its  origin  at  the  center 
of  mass  of  P1#  The  Xj  axis  lies  along  7^  (the  velocity  of 
Pj  about  the  Sun) ,  the  Yj  axis  points  toward  the  Sun,  and 
the  Zj  axis  completes  the  right-hand  system,.  Prom  the  figure, 


VQ  sin  *0  cos  i  ex  +  VQ  cos  *0  ey  + 


V0  sin  *0  sin  i  ez^ 


(2-5) 


where 


7P  —  Vp  ey 
M  M  Ai 


(2-6) 


Combination  of  Eqs  (2-5)  and  (2-6)  ,  and  use  of  the  dot 
product  to  find  an  expression  for  the  square  of  the  magnitude 
of  Vo,Pi ,  gives 

2  2 

V“Pi  c(Vo  sin  *0  cos  1  -  VpJ  +  (Vo  cos  ^o)  + 

(v0  sin  4>0  sin  i)  2  (2-7) 

The  value  of  VQ  must  satisfy  Eq  (2-7)  which  reduces  to 
the  following  quadratic  equation: 

Vj  -  2VPj  sin  $0  cos  i  V0  +  V*  -  vSp  «  0  (2-8) 


The  solutions  to  this  equation  are  given  by 
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V0  =  Vp  sin  cos  i  + 


sin24>0  cos2i  + 


(2-9) 


The  plus  sign  in  Eq  (2-9)  applies  to  orbits  which  lie 
outside  the  Earth's  orbit  about  the  Sun  (V  is  perihelion 
velocity  for  <J>0  «  90°)  j  the  negative  sign,  for  orbits 
which  lie  inside  the  Earth's  orbit  (VQ  is  aphelion  velocity 
for  $0  =  90°)  ,  is  used  for  most  solar  probe  trajectories. 

If  the  selected  values  of  <j>0  and  i  make  V0  a  complex 
number,  then  Vb  (or  V^p^)  is  not  large  enough  to  attain  a 
heliocentric  orbit,  and  the  vehicle  is  unable  to  escape  the 
Earth's  gravitational  field. 

Orbit  Characteristics .  Now  that  the  V0,  <J>0,  and  rPj 
(the  distance  of  the  Earth  from  the  Sun)  are  known,  the 
geometrical  characteristics  of  the  heliocentric  transfer 
orbit  (see  Pig.  8)  may  be  calculated. 

The  true  anomaly  at  launch  is  given  by  (Ref  21:62) 


tan  f0 


(rPj  vo  /Vo)  cos  *0  sin  *0 

1  “  (rPj  v5  /v0)  sin2$0 


(2-10) 


where  f0  is  measured  from  perihelion,  and  is  a  negative 
angle  for  elliptic  orbits  insida  Earth's  circular  orbit, 
that  is,  -180°  $  f0  $  0°, 

The  semi-major  axis  of  the  heliocentric  orbit  (Ref 
21:73)  is  given  by 
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a 


l 


2  “  rp.  V| 


(2-11) 


and  the  eccentricity  is 

el  =  ( Vo  /uo  ~  1]  *  sin2$o  *  ««2*0  (2-12) 


The  transfer  time  on  an  elliptic  orbit  is  found  by 
integrating  the  moment  of  momentum  equation  for  the  orbit¬ 
ing  body  (Ref  21j73).  The  result  of  this  integration  is 


/ - 

(  i -  \ 

/si 

2  arc  tan 

'  /l  -  e  tan  f  | 

J  M 

1  V  1  +  e  7  J 

1 

_ 

\v  / 

e  / 1  -  e2  sin  f 
1  +  e  cos  f 


(2-13) 


which  gives  the  time  measured  from  periapsis  to  the  true 
anomaly  f.  Inspection  of  this  equation  reveals  that  neg¬ 
ative  true  anomalies  (-180°  $  f  <  0°)  produce  negative 
time  values,  so  a  negative  sign  must  be  placed  in  front 
of  Eq  (2-13)  when  the  spacecraft  is  traveling  from  aphelion 
to  perihelion. 

For  convenience,  Eq  (2-13)  is  written  in  functional 
form,  that  is, 


t  »  t (a,  y,  e,  f) 


(2-14) 
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Thus,  for  direct-transfer  trajectories  in  this  chapter 


-  t(a1#  u0,  elt  f0) 


(2-15) 


which  gives  the  time  from  launch  to  perihelion. 

The  aphelion  of  the  heliocentric  transfer  orbit  is 
given  by 

rA  “  ai  d  +  ei)  (2-16) 

and  finally,  for  the  perihelion, 

qm  =  (1  -  ej)  (2-17) 


Computer  Programming 

As  shown  on  the  flow  chart  in  Fig.  9,  Program  1  (see 
Appendix  E)  consists  of  a  parametric  variation  of  Vfa,  <j>0, 
and  i,  respectively.  Characteristics  of  the  orbits  result¬ 
ing  from  the  choice  of  these  parameters  are  printed  in  rows 
of  data  which  vary  with  i  for  each  selected  value  of  and 
$0.  A  line  of  asterisks  across  the  page  of  output  is  used 
to  separate  the  results  for  each  choice  of  V^,  as  shown  in 
the  sample  printout. 

Control  provisions  are  made  in  the  program  for  cases 
in  which  the  values  of  V^,  $0,  or  i  are  such  that  escape 
from  Pj  is  not  possible.  The  increments  for  the  independent 
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variables  may  be  selected  as  desired,  and  in  this  program 
are 

Ai  =  5° 

A$  =  10° 
ro 

AVb  =  1000  ft/sec 
where  the  initial  values  are 

i  =  0° 

♦o  »  9°° 

Vb  =  36,000  ft/sec 

Special  precautions  must  be  taken  to  ^ind  the  true 
anomaly  f0  from  Eq  (2-10)  on  the  computer,  since  tan  f0  is 
a  double-valued  function.  The  arc  tangent  function  sub¬ 
routine  in  the  IBM  7094  system  is  designed  so  that  the 
output  angle  for  a  given  argument  is  less  than  +  ir/2  and 
greater  them  -  u/2.  Thus,  by  considering  the  double¬ 
valuedness  of  the  tangent  of  any  angle,  it  is  reasoned 
that 

f0  =  arc  tan  (tan  f0) ,  if  tan  f0  <  0 

(2-18) 

f0  =  arc  tan  (tan  f0)  -  it,  if  tan  f0  >  0 


An  arithmetic  statement  function  is  used  to  find  the 
time  on  an  elliptical  orbit.  The  function  is  defined  like 
Eqs  (2-13)  and  (2-14) ,  above,  and  it  appears  at  the  beginning 
of  each  program  in  the  thesis. 

In  general,  the  variable  names  in  the  programs  and  the 
output  correspond  to  the  notation  used  in  the  report.  The 
following  list  relates  the  notation  on  the  computer  output 
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with  symbols  defined  in  the  text: 

VB  *  Vb  XI  ^  i  El  ~  ex 

VlNFPl  *  V„p  FO  *  fQ  TlMOP  ^  Top 

PH  ^  4>0  A1  ^  3j  RP  *  qm 

/ 

Output  quantities  in  all  programs  are  given  in  the 
same  units  used  in  the  text. 

Conclusions 

The  graph  in  Fig.  10,  plotted  from  computer  output, 
illustrates  the  perihelion  which  can  be  achieved  with  a 
given  for  several  values  of  <{>0.  Fig.  11  is  plotted 
with  the  same  data,  except  that  mission  time  is  plotted 
as  a  function  of  $0  for  various  burnout  velocities. 
Examination  of  these  two  graphs  (both  for  transfers  in 
the  ecliptic)  indicates  some  interesting  properties  of 
direct-transfer  orbits. 

The  primary  reason  for  using  non-tangential  (<)>0  <  90°) 
transfers  is  to  decrease  the  mission  time.  Fig.  11  illus¬ 
trates  that  mission  time  is  strongly  influenced  by  small 
changes  in  <j>0  for  low  values  of  V^.  For  instance,  at  Vj-,  = 
37,000  ft/sec,  a  change  in  4> 0  from  90°  to  88°  decreases 
the  mission  time  from  140  to  120  days.  However,  at  = 
50,000  ft/sec,  this  same  change  in  $0  produces  a  mission 
time  decrease  of  only  0.2  days. 

As  shown  in  Fig.  10,  perihelion  is  significantly 
increased  when  <j>0  departs  from  90°,  For  example,  at  «= 
41,000  ft/sec  a  decrease  in  4> 0  from  90°  to  80°  reflects 
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an  increase  in  perihelion  from  0.50  AU  to  0.66  AU.  It 
appears  that  the  advantage  of  shorter  mission  time  for 
non-tangential  transfers  is  outweighed  by  the  disadvantage 
of  increased  perihelion  distance. 

Fig,  12  gives  the  Vb  requirements  for  direct-transfer 
out-of-ecliptic  missions  with  =  90°.  It  illustrates 
that  velocity  requirements  are  heavily  imposed  for  inclina¬ 
tion  angles  above  10°. 

From  Fig.  10,  tangential  transfers  to  0.5  AU  in  the 
ecliptic  require  burnout  velocities  less  than  41,000  ft/sec 
(Atlas/Centaur).  However,  more  sophisticated  boosters  must 
be  used  for  missions  closer  to  the  Sun?  that  is,  Saturn  V 
is  required  for  a  direct  mission  to  0.1  AU. 
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III,  Two-Dimensional  Gravity-Assisted  Trajectories 
Purpose  and  Scope 

The  burnout  velocity  and  cost  requirements  for  direct- 
transfer  solar  probe  trajectories  were  studied  in  the 
previous  chapter.  It  is  hoped  that  these  requirements  may 
be  relaxed  by  using  the  gravitational  attraction  of  one 
or  more  planets  to  reduce  the  spacecraft  perihelion. 

The  purpose  of  this  chapter  is  to  outline  the  method 
which  is  used  for  gravity-assisted  solar  probe  trajectories 
in  the  ecliptic  and  to  clarify  the  assumptions  made  in  the 
analysis.  A  detailed  derivation  of  the  equations  is  given 
in  Appendix  A. 

Two-Dimensional  Gravity  Assist 

It  is  assumed  that  a  gravity-assisted  trajectory  can 
be  separated  into  a  series  of  conic  orbits.  The  accuracy 
of  this  assumption  was  investigated  by  Sturms  and  Cutting 
(Ref  18)  who  computed  a  precise  trajectory  and  midcourse 
guidance  requirements  for  an  unmanned  probe  to  Mercury 
with  a  gravity  assist  at  Venus.  In  their  report  it  is 
shown  that  conic  trajectories  are  an  excellent  approxima¬ 
tion  to  the  actual  trajectories,  and  percentage  differences 
for  the  orbital  elements  are  within  2%. 

The  approach  taken  here  is  to  divide  the  entire  trajec¬ 
tory  into  four  conic  orbits : 

1.  Departure  Orbit  -  planetoeentric  hyperbola  at 
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departure  planet  P l . 

2.  Pre-Assist  Orbit  -  heliocentric  ellipse  from  Pj 
to  the  assist  planet  P2 . 

3.  Assist  Orbit  -  planetocentric  hyperbola  at  P2, 

4.  Post-Assist  Orbit  -  heliocentric  ellipse  after 
the  spacecraft  passes  P2 , 

Departure  Orbit  and  Pre-Assist  Orbit.  The  procedure 
which  was  used  to  determine  a  geocentric  escape  orbit  and 
a  heliocentric  transfer  orbit  in  Chapter  II  is  used  again  in 
this  chapter;  however,  the  heliocentric  transfer  orbit, 
defined  by  the  selection  of  and  $0,  is  now  called  a  pre¬ 
assist  orbit.  Since  this  analysis  is  two-dimensional,  the 
inclination  i  is  zero. 

The  point  of  intersection  A  of  t'he  pre-assist  orbit 
and  the  orbit  of  P2 ,  shown  in  Fig,  13,  is  assumed  to  be 
the  point  at  which  the  pre-assist  orbit  intersects  the 
SOI  of  P2 .  This  assumption  (no,  5  in  the  Introduction)  is 
more  meaningful  when  the  SOI  radius  of  P2  is  compared  with 
the  distance  of  P2  from  the  Sun.  For  example,  Venus  has  a 
SOI  radius  of  383,000  mi,  and  a  mean  orbital  radius  of 
67,100,000  mi.  If  the  SOI  radius  is  assumed  to  be  the 
chord  of  a  circular  sector,  it  subtends  a  heliocentric 
angle  of  only  0.33°.  A  similar  calculation  for  Jupiter 
reveals  that  its  SOI  radius  subtends  an  angle  of  3.5°, 

The  geometry  of  the  pre-assist  orbit  is  used  to 
calculate  true  anomaly  fj,  velocity  Vj ,  and  flight  path 
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angle  <f>1  at  point  A.  These  quantities  are  taken  as  pre¬ 
assist  orbit  conditions  at  the  time  the  spacecraft  enters 
the  SOI  of  P2. 

Assist  Orbit,  The  velocity  of  the  spacecraft  with 
respect  to  P2  as  it  enters  the  SOI  is 

v.4  =  v  -  Vp  (3-1) 

This  velocity  (at  an  infinite  distance  from  P2)  is  the 
inbound  HEV  for  a  hyperbolic  orbit  at  P2  called  the  assist 
orbit. 

From  considerations  on  hyperbolic  orbits  (Ref  16:96) 
it  is  known  that  the  outbound  HEV  vector  is  equal  to 
in  magnitude,  but  differs  by  an  angle  of  2v  in  direction 
(see  Fig.  14) .  This  angle  is  determined  by  choosing  a 
value  for  d,  the  DOCA;  after  this  choice  has  been  made, 
the  geometry  of  the  assist  orbit  is  completely  specified, 

Post-Assist  Orbit.  Point  A  in  Fig,  13  is  assumed  to 
be  fixed  in  heliocentric  space  during  the  time  that  the 
spacecraft  is  within  the  SOI  of  P2  (assumption  no,  4  of  the 
Introduction) .  It  has  been  shown  (Ref  10)  that  a  spacecraft 
spends  approximately  one  day  on  an  assist  orbit  in  the  SOI 
of  Venus.  The  above  assumption  is  reasonable  when  it  is 
considered  that  the  mean  daily  motion  of  Venus  is  1,602°. 

Thus,  neglecting  small  changes  in  the  position  and 
velocity  of  P2  about  the  Sun,  the  post-assist  orbit  can  be 
assumed  to  originate  at  point  A.  Also,  by  using  the  same 
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arguments  which  ware  applied  to  the  pre-assist  orbit,  point 
A  is  assumed  to  be  the  point  at  which  the  post-assist  orbit 
intersects  the  SOI  of  P2 . 

In  accordance  with  the  above  assumptions,  the  following 
expression  is  written  for  the  velocity  of  the  spacecraft 
with  respect  to  the  Sun  as  it  exits  the  SOI  of  P2  t 

V3  =  Vp  +  V.  (3-2) 

This  expression  and  the  geometry  of  the  assist  orbit  are 
used  to  find  the  true  anomaly  f3,  velocity  V3,  and  flight 
path  angle  (ij  on  the  post-assist  orbit  at  P2 .  Finally, 
these  conditions  are  used  to  determine  the  geometry  and 
perihelion  of  the  post-assist  orbit. 

Total  mission  time  T0p  is  defined  as  the  elapsed  time 
from  launch  at  Pj  to  perihelion  of  the  post-assist  orbit. 

The  time  spent  in  the  SOI  of  P1  and  P2  is  neglected,  since 
it  is  much  less  than  T0p.  Thus,  use  of  the  functional 
notation  in  Eqs  (2-13)  and  (2-14)  gives 

Tqp  13  T12  -  t  (a3 ,  e3 ,  f3)  (3-3) 

where  T12  is  the  time  from  Pj  to  P2 : 

T 1 2  "  t(al'  vO'  eir  fl>  “ 

b(a1 ,  Mq,  ej  ,  f0)  (3-4) 


Corollary 

As  explained  in  the  previous  section,  the  pre-assist 
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orbit  (determined  by  the  choice  of  and  <j>0)  defines  the 
geometry  of  the  velocity  vector  diagram  in  Fig.  15(a)  at 
the  time  of  entry  to  the  SOI  of  P2.  The  actual  point  of 
entry  is  not  specified  until  d  is  chosen.  Thus,  an  implicit 
result  in  the  above  analysis  is  that  a  given  pre-assist 
orbit  may  be  assumed  to  enter  the  SOI  at  any  point.  This 
corollary  .is  the  basis  for  the  third  assumption  in  the 
Introduction,  and  will  be  of  more  importance  in  the  three- 
dimensional  discussion  of  Chapter  VII. 

In  precision  trajectory  calculations  the  point  of 
entry  and  DOCA  are  determined  by  a  pre-assist  orbit.  During 
an  actual  flight  the  point  of  entry  must  be  controlled  by 
midcourse  guidance  and  small  deviations  in  the  launch  time. 

Geometry  of  the  Gravity-Assist  Maneuver 

Locus  of  Outbound  HEV  Vectors.  Since  iv»o|  =  | VM^ j 
a  circular  locus  of  vectors  may  be  drawn  for  a  given 
pre-assist  orbit.  As  shown  in  Fig.  15(b),  the  center 
of  this  circle  lies  at  the  tip  of  the  vector  ,  The 
selection  of  d  automatically  fixes  the  angle  v,  and  the 
resulting  V«0  vector  is  drawn  in  Fig.  15(c).  Finally,  in 
Fig.  15(d),  is  drawn  from  the  base  of  to  the  tip 
of  VMq. 

It  should  be  noted  that  physical  limitations  are 

imposed  on  the  angle  v ,  so  that  only  an  arc  of  the  circle 

in  Fig.  15(b)  represents  possible  V-  vectors.  The  maximum 

o 

turn  angle  2vmax  for  a  given  pre-assist  orbit  occurs  when 
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Fig.  15 

Gravity-Assist  Geometry 
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d  =  0  (grazing  pass  at  ?2)  .  This  can  be  seen  analytically 
by  inspection  of  Eqs  (A-18)  and  (A-19). 

Effect  of  Gravity  Assist.  The  effect  of  gravity  assist 
is  to  rotate  7ooi  through  an  angle  of  2v  in  the  plane  of 
motion.  The  direction  of  rotation  can  increase  or  decrease 
the  spacecraft  velocity  with  respect  to  the  Sun,  and  is 
controlled  by  the  side  of  the  planet  which  is  passed  by  the 
spacecraft  (dark-side  or  light-side  passage) . 

Some  interesting  observations  may  be  made  in  the 
velocity  vector  diagram  of  Fig.  15(d).  First  of  all,  the 
minimum  | V 3 |  for  a  given  pre-assist  orbit  is  attained  when 
V»0  is  directed  opposite  to  Vp^.  Similarly,  the  maximum 
1 V 3  i  is  attained  when  7«0  is  in  the  same  direction  as  \Tp  . 
Also,  the  maximum  heliocentric  velocity  change  |V3  -  ^  j 
is  achieved  when  v  =  vmaX  (or  d  =  0) . 

Another  observation  is  that  if  >  t^p2 I /  and  if 

v  (or  d)  can  be  chosen  within  its  range  of  physically 
possible  values  (d  >  0)  such  that  Vco0  lies  opposite  Vp  , 
then  V3  also  lies  opposite  to  Vp^ ,  and  a  retrograde  orbit 
is  obtained.  This  is  found  to  be  the  case  for  some  Jupiter- 
assisted  trajectories  (Ref  11)  . 

Visualization.  The  gravity-assist  maneuver  may  be 
visualized  by  separating  the  velocity  vector  diagram  of 
Fig.  15(d)  into  two  reference  frames:  one  with  respect  to 
the  Sun,  and  the  other  with  respect  to  P2.  This  procedure 
is  best  illustrated  with  a  hypothetical  example. 
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Suppose  that  a  gravity-assisted  trajectory  has  been 
calculated  and  that  the  above  procedures  have  been  followed 
to  obtain  the  velocity  vector  diagram  in  Fig.  16(a).  Since 
all  vectors  in  this  figure  are  drawn  relative  to  heliocen¬ 
tric  space,  they  may  be  referenced  to  the  P2  -  Sun  line  and 
translated  into  two  separate  diagrams  as  shown  in  Figs.  16(b) 
and  (c) .  All  angles  in  Fig.  16(a)  have  been  preserved  in  the 
translation,  and  now  the  entire  maneuver  may  be  easily  visu¬ 
alized  by  viewing  it  from  two  different  frames  of  reference. 

Comparison  with  Other  flethods 

The  method  used  for  gravity-assisted  trajectories  in 
this  study  is  similar  to  that  used  by  Minovitch  and  Niehoff 
(Refs  10,  11,  13).  Their  investigations,  however,  are  pri¬ 
marily  concerned  with  interception  of  a  target  planet  after 
a  close  approach  to  one  or  more  intermediate  planets.  The 
target  planet  rendezvous  imposes  a  constraint  which  need 
not  be  considered  in  solar  probe  missions. 

Minovitch  has  conducted  a  very  thorough  study  or  gravity- 
assisted  trajectories,  and  his  approach  includes  the  inclina¬ 
tions  and  eccentricities  of  planetary  orbits.  However,  the 
analysis  may  yield  physically  unrealizable  trajectories  with 
negative  DOCA  values. 

Niehoff  has  derived  expressions  for  the  maximum  velocity 
and  energy  change  at  the  assist  planet,  but  satisfaction  of 
these  expressions  constrains  the  geometry  of  the  assist 
orbit  which  also  fixes  a  DOCA  for  a  given  pre-assist  orbit. 
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Thus,  trajectories  which  have  maximum  velocity  or  energy 
change  are  not  necessarily  those  trajectories  which  have 
minimum  perihelia. 

Ross  (Ref  15)  has  devised  a  systematic  approach  to  the 
study  of  non-stop  interplanetary  round-trip  trajectories. 
His  analysis  includes  the  inclinations  ar.d  eccentricities 
of  planetary  orbits,  and  results  are  presented  in  design 
manuals  (Planetary  Flight  Handbook,  Vol.  Ill,  Parts  1,  2, 
and  3,  NASA  SP-35)  for  manned  missions  to  Venus  and  Mars. 

He  has  conducted  a  short  investigation  on  Venus-assisted 
solar  probe  missions  for  the  1965-66  time  period. 

The  method  developed  in  this  thesis  lends  itself 
nicely  to  feasibility  studies,  and  is  designed  so  that 
gravity-assisted  trajectories  may  be  calculated  as  a 
function  of  specified  burnout  conditions  at  Earth  and 
DOCA  at  the  assist  planet.  The  equations  are  derived  in 
a  general  form,  and  may  be  applied  not  only  to  solar  probe 
missions,  but  also  to  any  of  the  target-planet  or  round- 
trip  missions  which  were  studied  by  Minovitch,  Niehoff, 
and  Ross. 
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IV.  Venus  Assist  in  the  Ecliptic 

A  method  of  analysis  for  gravity-assisted  trajectories 
was  described  in  the  preceeding  chapter.  It  is  the  purpose 
of  this  chapter  to  apply  that  method  to  the  specific  problem 
of  a  Venus  assist  for  solar  probe  missions  in  the  ecliptic. 
This  problem  has  been  investigated  in  other  reports  (Ref3 
3,  11) ,  but  is  included  here  to  serve  as  a  foundation  for 
the  development  of  other  Venus-assist  missions  which  are 
studied  in  later  chapters. 

Also,  this  chapter  is  designed  to  outline  the  feasibil¬ 
ity  cf  gravity-ass > st  missions  and  some  of  the  additional 
considerations  which  must  be  made  in  the  selection  of  a 
trajectory  for  more  precise  design  study. 

Additional  Considerations  for  Gravity-Assist  Missions 

Guidance  Requirements.  The  most  significant  difference 
between  direct-transfer  and  gravity-assist  missions,  as  far 
as  hardware  is  concerned,  is  guidance  requirements,  A  grav¬ 
ity-assist  mission  would  have  no  hope  of  success  without 
some  form  of  midcourse  maneuvering  capability. 

This  problem  is  not  as  formidable  as  it  appears,  for  it 
hast  been  shown  that  guidance  requirements  for  the  gravity- 
assist  maneuver  are  modest,  and  are  within  the  capabilities 
of  present-day  technology,,  Sturms  and  Cutting  (Ref  18)  have 
shown  that  180  lb  of  additional  midcourse  guidance  equipment 
and  fuel  are  required  to  accomplish  a  Venus  assist  for  a 
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1120-lb  Mercury-probe  payload,,  A  design  study  (Ref  3)  on 
Venus-assist  missions  for  solar  probes  indicates  that  an 
80-lb  guidance  system  is  sufficient  for  a  460-lb  payload. 

From  this  discussion  it  is  assumed  that  the  payload 
for  Venus-assisted  solar  probe  missions  in  this  thesis 
V7eighs  approximately  250  to  300  lb  (200  lb  for  the  probe 
and  50  to  100  lb  for  additional  guidance  equipment) •  The 
probe  weight  is  based  on  findings  in  the  ICARVS  design 
study  (Ref  12*103),  where  a  150-  to  200-lb  spacecraft  is 
designed  for  a  direct-transfer  mission  to  0.1  AU. 

DOCA.  It  is  found  ir  this  and  all  of  the  following 
chapters  on  gravity-assist  missions  that  the  post-assist 
orbit  perihelion  decreases  with  a  decrease  in  DOCA  at  the 
assist  planet.  As  mentioned  in  the  previous  chapter,  a 
grazing  pass  represents  the  physical  lower  limit  for  DOCA, 
but  a  planetary  atmosphere  may  represent  an  additional 
limitation. 

Little  is  known  about  the  Venusian  atmosphere,  but  it 
has  been  estimated  (Ref  4*680)  that  the  atmospheric  pressure 
at  an  altitude  of  120  km  (65  NM)  above  the  surface  of  Venus 
is  about  2  x  10" 6  atm.  For  Earth  this  pressure  is  found  at 
an  altitude  of  about  107  km  (58  NM)  from  the  1959  ARDC  Model 
Atmosphere  (Ref  25*Chap.  II,  p,  3),  The  Earth’s  atmosphere 
effectively  ends  at  400,000  ft  or  65.8  NM  (reentry  altitude 
in  Ref  23*104).  From  this  data  it  should  be  safe  to  assume 
for  Venus  that  atmospheric  effects  are  insignificant  for  a 
DOCA  above  100  NM. 
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Communications ,  It  is  shown  in  the  two  design  studies 
mentioned  above  (Ref  3,  12)  that  a  solar  probe  must  have 
data  handling  and  storage  capability  for  the  transmission 
of  information  at  a  suitable  time  after  perihelion.  Thus, 
the  actual  mission  time  may  be  longer  than  that  defined  in 
this  report  (time  from  launch  to  perihelion) . 

The  communications  necessary  to  effect  a  midcourse 
maneuver  for  a  planetary  flyby  have  been  demonstrated  in 
the  Mariner  missions.  For  many  of  the  gravity-assisted 
trajectories  in  this  study  the  spacecraft-Sun  elongation 
angle  from  Earth  is  large  enough  to  allow  communication 
with  the  spacecraft  at  the  time  of  planetary  encounter. 

Computer  Programming 

As  shorn  on  the  flow  chart  in  Fig.  17  the  program 
for  gravity  assist  in  the  ecliptic  (Program  2  in  Appendix 
E)  consists  of  a  parametric  variation  in  Vb,  $0,  and  d, 
respectively.  The  first  portion  of  the  program  is  designed 
to  calculate  and  store  the  departure  orbit  and  pre-assist 
orbit  conditions  at  launch  for  various  values  of  Vb  and  <j>0. 

The  value  of  d  is  parametrically  varied  for  each  value 
of  Vb  and  $0,  and  a  line  of  printed  output  (sample  follows 
Program  2)  consists  of  trajectory  characteristics  for  a 
given  value  of  d. 

As  in  Program  1,  control  provisions  are  made  for  cases 
in  which  escape  from  Pj  is  not  possible.  Also,  provision 
is  made  for  the  cases  in  which  the  pre-assist  orbit  does 
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not  intersect  the  orbit  of  P2.  That  is,  the  perihelion  of 
the  pre-assist  orbit  must  be  at  least  as  small  as  the  orbit¬ 
al  radius  of  P£  (see  Eq  (A-6)). 

The  comments  above  Eq  (2-18)  regarding  the  arc  tan  sub¬ 
routine  and  the  calculation  of  fQ  in  Program  1  are  also 
applicable  in  Program  2.  These  considerations  apply  again 
to  the  calculation  of  fj  and  f3,  the  true  anomalies  of  the 
pre-  and  post-assist  orbits  at  P2,  respectively,-  That  is, 
the  programmed  expressions  for  f x  and  f3  may  be  obtained  by 
replacing  fQ  in  Eq  (2-18)  by  fj.or  f3. 

The  calculation  of  the  angle  ^  in  Eq  (A  -22  )  requires 
a  procedure  similar  to  that  used  in  the  calculation  of  f0. 
The  arc  sin  subroutine  on  the  IBM  7094  is  designed  so  that 
the  output  angle  for  a  given  argument  is  less  than  +  w/2 
and  greater  than  -  n/2,  Prom  Pig.  64  and  Eq  (A- 22),  the 
following  condition  must  hold: 

(ij>  +  2v)  >  v/2 t  if  Vi  sin  >  V_ 

Pl  (4-1) 

(*  +  2v)  <  it/2,  if  Vj  sin  <  VPi 

and  thus,  for  the  program, 

$  *  ir  -  arc  sin  [sin  (\|)  +  2v))  -  2v,  if 
Vj  sin  >  VPj 

$  «  arc  sin  [sin  ($  +  2v)J  -  2v,  if 
V  j  sin  4  j  *  Vp 


(4-2) 
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Notation  on  the  printed  output  for  Program  2  is  the 
same  as  that  used  in  Program  1  of  Chapter  II.  Additional 
and  revised  notation  is  as  follows* 


VO  ^  V 

VINF 

*  V«:. ,  V„0 

NU  *  v 

A  *  a 

DOCA 

'll  d 

E2  %  e2 

ET  ^  n 

F3  % 

f3 

PSI  *  9 

TH12  *  812 

A3  * 

a3 

TIM12  *  T 

FI  *  f  1 

E3  ^ 

e3 

TOP  'v  Tqp 

VI  %  Vj 

V3  * 

V3 

QM  -v  qm 

Mission  Profile 

A  trajectory  from  the  computer  data  is  illustrated  in 
the  mission  profile  of  Pig.  18.  An  Atlas/Centaur/TE-364-3 
launch  vehicle  is  selected  for  the  mission;  from  Fig.  4  it 
is  found  that  the  burnout  velocity  of  this  vehicle  is 
53,000  ft/sec  for  a  250  lb  payload.  A  tangential  departure 
is  used  (<j>0  «  90°)  in  order  to  obtain  the  smallest  possible 
perihelion  with  the  given  burnout  velocity. 

The  launch  date  for  this  mission  (30  August  1970)  is 
selected  froswoemputer  printout  of  Program  8  in  Appendix  E. 
Heliocentric  longitudes  of  Earth  and  Venus  at  the  time  of 
launch  are  taken  from  Reference  1.  The  angles  in  Fig.  18 
have  been  drawn  carefully  with  a  protractor,  but  are  not 
indicated  because  they  serve  only  to  clutter  the  drawing. 

Conclusions 

From  the  success  of  the  Mariner  Program  and  the  results 
of  the  reports  discussed  above  (Refs  3,  18),  it  can  be 
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concluded  that  gravity-assist  missions  to  Venus  are  within 
the  capability  of  present-day  technology. 

The  example  mission  indicates  that  a  perihelion  of 
0.161  AU  can  be  achieved  with  the  Atlas/CentaurA1E-364-3 
launch  vehicle  ($16  million}  and  a  Venus  assist;  whereas 
a  Saturn  V  ($125  million)  or  a  Saturn  I/Centaur/Kick  (no 
cost  data,  but  Saturn  I/Centaur  alone  costs  $41  million) 
is  required  for  a  direct-transfer  to  this  perihelion. 

A  graph  of  perihelion  versus  burnout  velocity  for 
varying  DOCA  in  Fig.  19  indicates  that  the  performance 
(achieved  perihelion)  of  gravity  assist  at  Venus  changes 
only  slightly  for  DOCA  variations  below  1000  NM.  Thus, 
the  performance  results  for  a  DOCA  of  100  NM  are  fairly 
representative  of  the  results  for  any  DOCA  less  than  1000 
NM. 

A  plot  for  direct-transfer  missions  (dot-dashed  line) 
is  shown  in  the  figure.  A  comparison  of  the  curves  reveals 
that  a  DOCA  above  100,000  NM  offers  no  significant  advantage 
in  solar  probe  missions.  A  horizontal  measurement  on  the 
graph  for  any  given  perihelion  below  0,35  AU  indicates  that 
the  reduction  in  V^  is  approximately  6500' ft/sec,  that  is 

Vi.  -  Vb  a  6500  ft/sec. 

Direct  DOCA»0 

Fig.  19  also  illustrates  that  the  most  significant 
reductions  in  perihelion  are  achieved  with  the  lower 
values  of  Vfa  (less  than  55,000  ft/sec).  Thii  point  is 
further  illustrated  with  the  plot  of  perihelion  versus 

, 
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DOCA  in  Fig.  20,  The  horizontal  dashed  lines  in  this 
figure  indicate  perihelion  for  direct-transfer  missions t 
these  lines  are  asymptotes  for  the  solid  lines  which  are 
for  gravity-assist  missions. 

Mission  time  for  any  and  DOCA  can  be  found  in 
Fig.  21,  This  plot  indicates  that  for  a  given  the 
mission  time  is  relatively  invariant  with  DOCA.  A  cross¬ 
plot  cf  this  graph  is  found  in  Fig.  22  for  a  DOCA  of  100 
NM.  Also  shown  in  Fig.  22  is  mission  time  for  direct- 
transfer  orbits.  It  can  be  seen  that  the  mission  time 
requirements  for  a  Venus  assist  are,  for  all  practical 
purposes,  about  the  same  as  direct  transfer. 
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Perihelion  vs.  DOCA  for  Various  Burnout  Velocities 


10,000  20,000  30,000 

Fig.  21  DOCA  (NM) 

Mission  Time  vs.  DOCA  for  Various  Burnout  Velocities 
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Fig,  22  Vjj  (ft/sec  x  10 3 ) 

Mission  Time  for  Venus  Assist:  Compared  with  Direct  "Transfer 
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V .  Venus-Mercury  Combination  Assist 

The  theory  developed  in  Appendix  A  and  the  Venus-aesist 
mission  in  the  preceeding  chapter  is  easily  extended  to  the 
case  of  a  second  assist  at  Mercury.  The  post-assist  orbit 
for  the  Venus  assist  is  treated  as  a  pre-assist  orbit  for 
the  Mercury  assist,  and  the  procedures  for  a  single  assist 
are  repeated.  Thus,  each  Venus-assisted  trajectory  and  a 
chosen  DOCA  at  Mercury  determines  a  combination-assist 
trajectory. 

General  Method 

The  general  method  may  be  described  by  reference  to 
the  equations  derived  in  Appendix  A.  First  of  all,  a 
value  is  selected  for  Vj-,,  $0,  and  d  at  the  first  assist 
planet  P2.  This  determines  the  post-assist  orbit  charac¬ 
teristics  a3  and  e3;  if  the  post-assist  orbit  intersects 
the  orbit  of  a  second  assist  planet,  P3,  it  is  treated  as 
a  pre-assist  orbit  at  P3,  and 

r  •«  rp  (5-1) 

p3 


Eqs  (A-8)  through  (A-13)  are  used  to  find  the  orbit 
conditions  at  intercept  of  P3: 


vi 


(5-2) 
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I 

i 


Sin  $5 


~  «f  (1  -  e j) 

5?3  ^2a3  “  rP3* 


(5-3) 


tan  fs 


jrP3  vf  /»*0)  cos  $5  sin  $5 
1  "  (rP3  V1  /v0)  sin2*s 


(5-4) 


The  heliocentric  transfer  angle  and  time  from  P2  to  P3, 
respectively,  are  given  by 

023  a  fs  “  f3  (5-5) 


*23  *=  t  <a3#  U0»  e3#  U)  “ 

t  <a3,  e3,  f3)  (5-6) 


To  determine  the  assist  orbit  at  P3,  as  shorn  in  Fig, 
23,  use  is  made  of  Eqs  (A-14)  and  (A-15)  as  follows* 

V2.  ■  V?  +  V|  - 

»i5  5  p3 

(5-7) 

2V5  Vp  sin  $5 
3 

and  a4  a  ypg  /v£ig  (5-8) 

After  selection  of  d3,  the  -BOCA 
at  P3,  the  eccentricity  of  the 
second  orbit  is  calculated  as 
in  Eg  (A-18) t 
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e 


4  “  ~  (*P 


a4  i  s3 


Rr>  +  a,.  +  d, 


(5-9) 


and 


sin  v4  «  i. 
-  6 1 


(5-10) 


The  final  post-assist  orbit  conditions  at*P3  found 
with  Eqs  (A-22)  through  (A-26) : 


«  arc  sin  V s  ...  cos  4> s 

-  2v4 

(5-11) 

Lv«is  J 

v„  =  v„. 
o?  is 

(5-12) 

V£3  +  V*o7  ‘  SVP3  v-o7  008  *» 

(5-13) 

V* 

cos  $7  *=  — sin  4>4 

(5-14) 

tan  f7 


rP3  v$  /»*»)  cos  *7  sin  *7 
1  -  (rp  V2  /*J  sin2$7 


(5-15) 


Finally,  the  post-assist  orbit  characteristics  after 
the  second  assist  are  determined  from  Eqs  (A-27)  through 
(A-31) , 


a5 


(5  16) 


e|  «  (rp3  vf  /p0  -  l)  2  sin2$7  +  cos2+7 


(5-17) 
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T0P  ~  +  T*3  “ 

(5-18) 

ra5  -  a5  (1  +  e5> 

(5-19) 

and  for  the  final  perihelion. 

’  rp5  "  as  <1  -  *s> 

(5-20) 

Computer  Programming 

The  program  for  combination-assisted  trajectories 
(Program  3)  is  illustrated  with  the  flow  chart  in  Fig.  24. 
The  first  half  of  Program  3  is  identical  to  Program  2,  so 
that  single-assist  trajectories  are  generated  by  a  paramet¬ 
ric  variation  of  Vfaf  <pQ,  and  d.  Each  trajectory  calculated 
in  the  first  half  of  the  program  is  used  to  generate  a  set 
of  combination-assisted  trajectories  by  a  parametric  varia¬ 
tion  in  d3,  the  DOCA  a>.  P3. 

In  addition  to  the  control  provisions  made  in  Program 
2,  a  provision  is  made  for  the  case  in  which  the  post-assist 
orbit  from  P2  does  not  intersect  the  orbit  of  P3.  The  per¬ 
ihelion  of  the  post-assist  orbit  from  P2  must  be  at  least 
as  small  as  the  orbital  radius  of  P3. 

A  sample  of  the  computer  printout  is  given  in  Appendix 
E  after  Program  3.  The  first  list  of  results  for  each  Vjj 
and  is  a  set  of  Venus-assisted  trajectories  for  varying 
d.  Below  this  is  a  listing  of  Venus-Mercury-assisted  tra¬ 
jectories  for  varying  d3  for  each  value  of  d. 

The  notation  on  the  printed  output  is  the  same  as  that 
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used  in  Programs  1  and  2  with  the  following  additions  or 
changes : 


N  *  V 

T23  *  T23 

E5  *  e5 

T12  ^  T12 

QMM  *  q_ 
mm 

V7  *  V7 

TH23  ^  023 

D3  *  d3 

N4  *  v4 

F5  *  f5 

F7  -v  f? 

PSI4  * 

VINF5  -v  V.ig,  V»05 

A5  'x*  ag 

THOP  *  0 

i 

TOP  'v  T0p  (Venus-Mercury  Assist) 

TTOT  *  T0p  (Venus  Assist  with  no  Mercury  Assist) 


Eccentricity  of  Mercury’s  Orbit 

When  Program  3  is  applied  to  the  Venus-Mercury  combina¬ 
tion-assist  problem,  it  is  necessary  to  consider  the  eccen¬ 
tricity  of  Mercury's  orbit,  since  it  is  not  nearly  circular. 

In  order  to  obtain  a  range  of  possible  results,  two 
cases  are  studied:  Mercury  at  perihelion  r^  and  Mercury 


at  aphelion  r. 


Thus,  for  the  perihelion  case, 


rP  =  r« 

p3  (5-21) 

and  for  the  velocity  cf  Mercury  (Ref  8:Chap.-  II,  pp.  11,12, 
Eqs  (1-134)  and  (L-153) ) , 


iilsLfsSL 


r  +  r  ] 


For  the  aphelion  case, 


a 


(5-22) 


(5-23) 
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and  for  the  velocity, 


2yo  rPg 

ag  (rag  +  ^ 


(5-24) 


Conclusions 

The  three  missions  which  have  been  studied  thus  far 
(direct-transfer,  Venus-assist,  and  Venus-Mercury  assist) 
are  co?npared  in  Pig.  25  on  a  plot  of  perihelion  versus  burn¬ 
out  velocity.  Results  for  both  the  Venus-Msrcury-at-aphelior 
and  Venus-Mercury- at-perihelion  missions  are  shown.  Of  the 
two  cases,  smaller  perihelia  are  achieved  when  Mercury  is 
at  aphelion,  bttfc  the  difference  is  almost  insignificant a 

It  can  be  concluded  immediately  from  the  figure  that 
the  Venus-Mercury  combination-assist  missions  are  not 
worthwhile.  Not  only  is  the  reduction  in  perihelion  insig¬ 
nificant  when  compared  to  the  single  Venus-assist  missions, 
but  also  the  added  problems  of  guidance  at  a  second  assist 
planet  and  launch  opportunities  make  these  missions  more 
trouble  than  they  are  worth. 

Thus,  the  only  useful  result  of  this  chapter  is  that 
a  method  has  been  outlined  for  the  calculation  of  Venus- 
Mercury-probe  trajectories,  and  equations  have  been  derived 
for  duel  combination-assisted  trajectories.  However,  the 
analysis  has  contributed  little  for  solar  probe  missions. 


Comparison  of  Direct-Transfer,  Venus 
Mercury-Venus  Combination  Assist  M 
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VI.  Multiple  Venus  Assist 

The  theory  for  gravity-assisted  trajectories  in  Appendix 
A  was  extended  in  the  preceeding  chapter  for  the  case  of  a 
combination  of  gravity  assists  at  several  planets.  In  this 
chapter,  the  same  theory  is  extended  for  another  case  - 
multiple  gravity  assists  at  a  single  planet. 

The  problem  here  is  to  determine  a  DOCA  at  the  assist 
planet  (for  given  values  of  burnout  velocity  V^  and  departure 
flight  path  angle  <j>0)  which  will  cause  the  spacecraft  to 
intercept  P2  a  second  time  at  the  same  point  in  space.  As 
shown  below,  this  requires  a  post-assist  orbit  whose  period 
is  m/n  times  the  orbit  period  of  P2;  the  quantities  m  and  n. 
are  integers,  where  m  is  the  number  of  solar  circuits  of  P2 
and  n  is  the  number  of  solar  circuits  of  the  spacecraft  be¬ 
fore  the  second  intercept. 

This  procedure  may  be  repeated  for  additional  passes  at 
P2  until  the  spacecraft  aphelion  is  too  small  for  interception. 
Of  course,  the  DOCA  for  a  final  pass  if,  arbitrary,  just  as  it 
was  for  a  single  pass. 

The  theory  and  computer  programs  in  this  hapter  are 
developed  for  multiple  assist  at  any  planet;  but  application 
is  made  only  to  a  double  and  a  triple  assist  ah  Venus,  It 
has  already  been  shown  that  little  benefit  is  derived  with  a 
Mercury  or  Mars  assist;  a  multiple  assist  at  Jupiter  is  imprac¬ 
tical  because  of  excessive  flight  time  and  also  because  Sun 
impact  can  be  achieved  with  just  a  single  assist. 
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Double  Assist  at  Venus 

As  stated  above,  the  pre-assist  orbit  is  determined  by 
selecting  Vb  and  <j>0.  The  next  parameter  fcr  a  double  assist 
is  the  orbit  ratio  m/n.  If  the  spacecraft  and  P2  are  to 
intercept  a  second  time  at  the  same  point  in  space,  then 
necessarily, 


n  x 


sc 


=  ra  t, 


(6-1) 


where  tsc  is  the  first  post-assist  orbit  period  and  tp^  is 
the  orbit  period  of  P2. 

Thus,  for  the  orbit  ratio  selected,  the  required  post¬ 
assist  orbit  period  is 


T 


SC 


(6-2) 


where  xp^  is  known  from  planetary  data.  From  Kepler’s  Laws, 
this  also  specifies  the  semi-major  axis  of  the  post-assist 
orbit: 


a3  = 


(6-3) 


The  heliocentric  velocity  at  P2  on  the  post-assist 
orbit  is  found  with  the  vis-viva  integral: 


vi  -  \ 


(6-4) 
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As  shown  in  Fig.  26  there  are  two  possible  directions  for  ^3. 
The  ^3  case  (dotted  line)  is 
discarded,  because  it  serves  to 
increase  rather  than  decrease 
the  perihelion  of  the  post¬ 
assist  orbit.  If 

IV3 1  <  l*P2  "  *»ii  (6-5) 

then  the  trajectory  cannot  be 
realized  and  another  orbit 
ratio  must  be  selected. 

Conservation  of  energy  of 
the  flyby  orbit  gives 

V„0  =  V„.  (6-6) 

and  the  flight  path  angle  $  3 ,  is 
determined  from  the  Law  of  Cosines: 


sin  $  3 


2V,  VT 


v?  +  vi2 


(6-7) 


The  angle  i>2  in  the  figure  is  given  by 

V3 

sin  <|>2  =  _ _  cos  §  a 

and  from  the  Law  of  Sines,  the  turn  angle  is 


2v?  *»  arc  sin 


Vi 


COS 


-  h 


(6-8) 


(6-9) 
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Thus,  Eqs  (A-15)  and  (A-I9)  gi"e 

Up2 

&2  33  (6-10) 

V«, 

x 

e2  a  _JL— .  (6-11) 

sin  v2 

and  finally,  from  Eqs  (A-16)  and  (A-17) ,  the  required  DOCA 
at  P2  is 

d  «  a2  (=2  -  1)  -  Rp  (6-12) 

If  a  negative  value  is  obtained  for  d,  then,  again,  another 
value  of  m/n  must  be  chosen. 

The  remaining  characteristics  of  the  post-assist  orbit 
are  calculated  with  Eqs  (A-26)  through  (A-31) . 

If  only  minor  quidance  corrections  are  made,  the  first 
and  second  intercepts  may  be  assumed  to  occur  at  the  same 
point  in  space.  Thus,  the  post-assist  orbit  for  the  first 
pass  becomes  a  pre-assist  orbit  for  the  second  pass  with  the 
new  pre-assist  conditions  V3  and  $3  (see  Fig.  27).  The 
remainder  of  the  problem  is  now  identical  to  the  single 
assist  mission,  and  the  DOCA  for  the  second  pass,  d2 ,  may  be 
selected  arbitrarily.  The  post-assist  conditions  (V5,  $5) 
and  the  final  post-assist  orbit  characteristics  are  computed 
with  Eqs  (A-18)  through  (A-31) ,  The  final  perihelion  is 
given  by 
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^  =  a5  (1  -  e5) 

It  should  be  noted  here  that 
subscript  notation  in  this  chapter 
has  a  different  meaning  than  that 
used  in  Chapter  V. 

Computer  Program,  The 

computer  program  for  a  double 

assist  (Program  4) ,  illustrated 

in  the  flow  chart  of  Fig,  28, 

does  not  entirely  follow  the 

method  outlined  in  this  chapter. 

The  values  of  V,  and  $  are 
b  o 

parametrically  varied  as  before, 
but  d  is  found  by  iteration 
rather  than  direct  calculation 
for  given  values  of  m  and  n. 

A  value  is  selected  for  d,  and  the  first  assist  orbit 
is  computed?  next,  m  is  selected,  and  n  is  varied  until 
Eq  (6-1)  is  satisfied.  If  no  value  of  n  is  found  to  satisfy 
the  equation,  the  next  highest  integer  is  selected  for  m, 
and  lo  on  until  all  values  of  m  and  n  up  to  a  limiting 
integer  have  been  tried  for  the  value  of  d»  If  no  combina¬ 
tions  of  m  and  n  produce  the  desired  orbit,  then  d  is 
increased  by  an  increment,  and  the  procedure  is  repeated 
until  all  combinations  of  m  and  n  are  found  with  a  corre¬ 
sponding  value  for  d. 


(6-13) 
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This  procedure  demands  more  computer  time  than  would 
be  required  for  the  method  developed  in  the  chapter,  but  it 
also  avoids  a  major  rewrite  of  the  basic  program  which  is 
used  in  the  two  previous  chapters. 

In  this  program,  m  was  allowed  to  vary  from  1  to  5,  and 
n  from  1  to  9.  Obviously,  the  trajectories  with  low  values 
of  m  have  the  more  satisfactory  mission  times.  It  is  found 
in  the  computations  that  Venus-assisted  trajectories  are  not 
physically  realizable  for  m  or  n  equal  to  one. 

Notation  changes  or  additions  on  the  output  of  Program 
4  are  j 

D0CA1  -V  d  NSP  *  n  VINF3  -v  ,  V»n 

13  0  3 

TAUP2  -v.  Tp2  TAUP  ^  Tgc  T1M0P  *  Mission  Time 

MPl  *  m  D0CA2'v  d2  QM2  *  Final  Perihelion 

T  'v  Time  from  Launch  to  Perihelion  of  First  Pre- 
Assist  Orbit 

THTOT  -v  Heliocentric  Angle  from  Launch  to  Final 
Perihelion 

Mission  Profile.  The  mission  profile  in  Fig.  29  illus¬ 
trates  the  trajectory  for  a  burnout  velocity  of  42,000  ft/sec 
and  a  tangential  departure.  From  Fig.  4  this  trajectory  may 
be  achieved  by  an  Atlas/Centaur  launch  vehicle  with  a  300-lb 
payload.  This  payload  allows  for  additional  guidance  equip¬ 
ment  which  may  be  required  for  a  second  pass  (Ref  3) . 

A  significant  advantage  of  this  trajectory  is  that  the 
semi-major  axis  of  the  first  pre-assist  orbit  (0.728  AU)  is 
very  nearly  that  of  Venus’  orbit  (0.723  AU) .  The  value  of 
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Fig.  28 

Flow  Chart  for  Program  4 
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Pig.  29 

Mission  Profile  for  Double  Venus  Assist 
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V^  which  gives  exactly  aj  ^  0.723  AU  may  be  found  with  a  few 
calculations,  but  suppose  that  this  has  been  done.  Then  by 
Kepler's  Third  Law,  the  period  of  the  pre-assist  orbit 
matches  the  period  of  Venus.  Thus,  if  midcourse  maneuvers 
fail  to  establish  the  initial  intercept  at  Venus,  a  second 
attempt  may  be  made  one  Venus  year  later. 

Conclusions.  In  the  double-assist  mission  profile  it 
was  found  that  a  perihelion  of  0,201  AU  can  be  established 
with  the  Atl:vs/Centaur  ($11.9  million).  A  direct-transfer 
to  this  perihelion,  from  Figs.  4  and  10,  requires  a  much 
larger  booster  such  as  the  Saturn  I/Centaur  with  an  added 
kick  stage  or  a  Saturn  V. 

The  double-assist  missions  not  only  reduce  booster 
vehicle  costs,  but  also  provide  the  opportunity  to  explore 
two  regions  of  the  solar  system  with  a  single  launch.  For 
the  mission  profile  above,  data  can  be  obtained  from  the 
0.3  AU  region  after  the  first  pass  at  Venus  and  from  the 
0.2  AU  region  after  the  second  pass. 

Double-assist  missions  with  different  orbit  ratios  are 
compared  in  Fig.  30.  The  left-hand  end  of  each  curve  repre¬ 
sents  the  smallest  value  of  V^  with  which  the  desired  orbit 
ratio  can  be  achieved/  the  right-hand  end  represents  an 
upper  limit  on  d  (5000  NM)  in  the  computations.  If  d  is 
increased,  the  curves  may  be  extended  as  far  to  the  right 
as  desired. 

From  this  graph  the  most  promising  orbit  ratio  is  2/3, 
because  it  offers  short  mission  times,  and  gives  physically 


GA/AE/67-4 


realizable  trajectories  for  almost  the  entire  spectrum  of 
Vb  which  is  used  in  Venus-assist  missions.  The  3/5  ratio 
is  also  attractive  because  of  the  small  perihelia  which 
can  be  achieved;  however,  Vb  for  these  missions  must  be 
greater  than  55,000  ft/sec  which  automatically  demands  a 
more  expensive  booster.  The  5/8  missions  are  discounted 
for  their  unreasonable  time  requirements. 

One  additional  note  of  interest  on  this  graph  is  that 
perihelion  distance  decreases  (for  a  given  Vb)  with  the 
value  of  the  orbit  ratio. 

Results  of  the  2/3  missions  are  illustrated  on  the 
graphs  in  Figs.  31  through  34.  In  Fig.  31  double-assist 
missions  are  compared  with  direct-transfer  missions  for 
various  values  of  d2.  The  most  significant  reductions  in 
perihelia  are  achieved  for  the  lower  values  of  V^. 

The  value  of  d  required  to  establish  a  post-assist 


value  of  Vb.  From  this  graph,  it  can  be  seen  that  d 
increases  with  Vbj  therefore,  the  effect  of  the  first 
assist  diminishes  as  Vb  increases,  and  only  the  second 
assist  serves  to  decrease  the  perihelion.  This  trade-off 
property  is  also  evident  in  Fig.  33,  where  the  double¬ 
assist  missions  (solid  lines)  are  compared  with  the  single¬ 
assist  (dot-dashed  lines)  and  direct-transfer  missions 
(dashed  lines)  for  a  few  values  of  Vb#  From  this  graph, 
if  *  55,000  ft/sec,  then  a  single-assist  mission  is 
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more  beneficial  than  a  double-assist  mission  for  d2  > 

4000  NM. 

An  inherent  disadvantage  in  double-assist  missions 
is  the  time  requirements,  as  illustrated  in  Fig,  34  for 
m/n  =»  2/3  and  d2  *  100  NM,  Double-assist  missions  require 
approximately  550  days  or  1,15  years,  whereas  direct 
missions  require  less  than  160  days. 


O 
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Comparison  of  Venus  Double-Assist  Missions  with  Various 
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Required  dj  for  Double-Assist 
Missions  to  Venus 
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V,  (ft/sec  x  10 3 ) 
Fig.  34  b 

Time  Requirements  for  Double-Assist  Missions 
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Triple  Assist  at  Venus 

The  analysis  for  a  double  assist  extends  easily  to  the 
case  of  a  triple  assist.  The  method  used  to  determine  the 
first  DOCA  for  a  second  intercept  is  again  used  to  determine 
the  second  DOCA  for  a  third  intercept. 

For  a  triple  assist#  two  relations  similar  to  Eq  (6-1) 
or  (6-2)  must  be  satisfied 


T 


SC} 


(6-14) 


x 


SC' 


(6-15) 


where  m^  is  the  number  of  solar  circuits  of  P2 ,  and  n^  is 
the  number  of  solar  circuits  of  the  spacecraft  after  the 
ith  intercept  (i  =  1,2) ,  From  the  double-assist  procedures 
then,  it  follows  that 


(6-16) 


(6-17) 


for  the  semi-major  axes  of  the  first  and  second  post-assist 
orbits,  respectively.  The  heliocentric  velocity  at  P2  on 
the  first  and  second  post-assist  orbits,  respectively,  is 
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v! 


<6-18) 


V|  =■  Mo 


(6-19) 


The  HEV  is  the  same  for  each  pass  due  to  conservation  of 
energy.  Thus,  the  flight  path  angles  $3  and  $5  may  be 
calculated  with  relations  similar  to  Eq  (6-7) , 

The  hyperbolic  assist  orbit  at  each  pass  is  found  with 
Eqs  (6-8)  through  (6-11) ,  which  give 


e2  =»  — ; - 

sin  v. 


(6-20) 


sin  V4 


(6-21) 


Finally,  the  required  DOCA  for  the  first  two  passes,  respec* 
tively,  are  given  by 


dj  a  a2  (e2  -  1)  -  Rp  (6-22) 

r2 

d2  =  a4  (e4  -  1)  -  Rp  (6-23) 


Again,  if  a  negative  value  is  obtained  for  either  DOCA,  then 
another  combination  of  or  n^  must  be  selected. 

The  DOCA  for  the  final  pass,  d3,  may  be  chosen  as 
desired,  and  the  problem  again  takes  the  form  of  a  single 
assist.  Thus,  the  selection  of  d3  determines  the  third 
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post-assist  orbit  conditions  V?  and  <{ >7  shown  in  Fig.  35. 
Finally,  the  post-assist  orbit 


characteristics  are  found  with 

,  To 

Eqs  (A-26)  through  (A-31) ,  and 

W  ♦  5  J  Sun 

the  final  perihelion  is 

u 

qm  =  a7  (1  -  e7)  (6-24) 

3 

V 

\ 

V\  V 

Computer  Program.  Program  5 

7  3 

(for  a  triple  assist)  is  summa- 

\ 

rized  in  the  flow  chart  of  Fig. 

Fig.  35 

36.  The  program  is  identical  to 

Program  4,  except  that  a  section 

Vector  Geometry  for 

has  been  added  to  calculate  the 

Third  Assist 

second  DOCA  for  a  third  intercept. 

After  the  first  post-assist  is  computed  (dj,  mi,  ni 
determined) ,  a  value  is  selected  for  d2 ,  and  m2  and  n2  are 
varied  up  to  some  limiting  integer  (5  and  9,  respectively) 
until  Eq  (6-15)  is  satisfied.  If  no  orbit  is  found  for  all 
combinations  of  m2  and  n2 ,  then  d2  is  increased  by  an  incre¬ 
ment,  and  so  on  as  in  Program  4.  After  Eqs  (6-14)  and  (6-15) 
have  been  satisfied,  a  final  post-assist  orbit  is  computed 
for  various  values  of  d3. 

Notation  on  the  output  is  as  follows: 


D0CA2  'v-  d2 

MP2  ^ 

m2 

N6  'v*  v6 

TAUP  -v  tSCi 

NS2  ~ 

n2 

QM3  'v  qm3 

TAUPP  -v  rgC2 

D0CA3 

*  d3 

TIM02  Mission  Time 
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Fig.  36 

Flow  Chart  for  Program  5 


GA/AE/67-4 


THT02  ^  Heliocentric  Angle  from  Launch  to  Final 
Perihelion 

Mission  Profile.  A  triple-assist  mission  is  illustrated 
in  Fig.  37  for  a  tangential  departure  and  a  burnout  velocity 
of  53,000  ft/sec.  An  orbit  ratio  of  (2/3)  (5/9)  is  selected 
so  that  the  second  assist  may  be  completed  early  in  the 
flight;  thus,  a  partial  success  ip  assured  if  capability  for 
a  third  intercept  is  lost  during  the  mission. 

The  Atlas/Centaur/TE-364-3  is  used  for  this  mission  with 
a  payload  of  250  lb.  It  is  assumed  here  that  this  payload 
weight  can  be  realized  for  a  triple-assist  by  the  1972  time 
period.  A  heavier  package  would  require  either  a  larger 
booster  or  an  additional  kickstage  (Ref  6,  Ref  12:89-101). 

Conclusions.  In  the  triple-assist  mission  profile  it 
was  found  that  a  perihelion  of  0.096  AU  can  be  attained  with 
a  burnout  velocity  of  53,000  ft/sec.  A  direct  transfer  to 
this  perihelion  requires  a  burnout  velocity  of  68,000  ft/sec; 
thus,  triple  assist  has  reduced  the  booster  vehicle  require¬ 
ment  from  Saturn  V/Kick  Stage  to  Atlas/Centaur /TE-364-3.  This 
mission  also  allows  the  exploration  of  three  regions  (0.183  AU, 
0,134  AU,  0,096  AU)  near  the  Sun  with  a  single  launch. 

Various  orbit  ratios  are  compared  in  Fig.  38,  and  again 
the  left-hand  end  of  the  curves  represents  a  physical  lower 
limit  on  for  each  orbit  ratio,  and  the  right-hand  end  is 
limited  by  the  highest  values  of  d1  and  d2  in  the  computa¬ 
tions  (5000  NM  for  each) .  The  (2/3)  (5/9)  orbit  ratio  is 
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Launch  Date*  6  Apr  1972  (JD  2441412.8) 
Launch  Vehicle*  Atlaa/Centaur/TE-364-3 
Vjj  »  53,000  ft/sec,  $0  »  90° 

Payload  ®  250  lb 
T12  «  52.7 

Tqp  a  1659.2  days  a  4 .54  years 

Orbit  Ratio  „  1  £ 

3  9 


di  =  2315  NM 
d2  =  540  NM 
d3  •  100  NM 
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Mission  Profile  for  Triple  Venus  Assist 
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the  most  favorable  for  two  reasons:  mission  time  is 
comparitively  short,  and  trajectories  are  realizable  for  a 
wide  range  of  V^.  From  the  figure  there  are  three  orbit 
ratios  which  yield  smaller  perihelia,  but  these  cases  are 
discarded  because  of  long  mission  times  ( (5/8)  (5/9)  and 
(3/5)  (5/9))  or  high  burnout  velocity  requirements  ((4/7) 
(1/2)). 

Results  for  the  (2/3)  (5/9)  missions,  illustrated  in 
Fig,  39  through  42,  indicate  the  same  trends  which  were 
found  in  the  2/3  double-assist  missions. 

Triple-assist  and  direct-transfer  missions  are  compared 
in  Fig,  39;  it  is  found  that  the  perihelion  may  be  reduced 
by  as  much  as  0.28  AU  with  a  triple  assist. 

The  required  values  of  the  first  and  second  DOCA  are 
shown  in  Fig.  40,  The  value  of  dx  increases  with  Vb,  where¬ 
as  d2  decreases  with  Vb.  If  it  turns  out  that  the  multiple- 
assist  navigation  system  is  limited  to  one  particular  DOCA 
for  each  intercept,  then  a  burnout  velocity  of  48,500  ft/sec 
must  be  used  and  dj  =  d2  =  d3  -*  1000  NM.  This  trajectory  has 
a  perihelion  of  0,125  AU, 

The  trade-off  property  for  double-assist  missions  is 
also  found  in  the  triple-assist  missions.  This  is  illustrated 
in  Fig.  41,  where  the  triple-assist  missions  (solid  lines) 
for  any  value  of  d3  are  compared  with  double-assist  missions 
(dot-dashed  lines)  with  d2  =  100  NM.  At  a  burnout  velocity 
of  45,000  ft/sec,  the  perihelion  for  a  triple  assist  is 
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smaller  than  the  perihelion  for  a  double  assist  only  if 
d3  <  11,000  NM. 

The  time  requirements  for  triple-assist  missions,  shown 
in  Fig.  42,  indicate  that  these  trajectories  are  of  question¬ 
able  value.  Mission  time  is  about  1670  days  (4,6  years) 
which  demands  a  spacecraft  reliability  well  beyond  present- 
day  capability. 

It  is  probable  that  the  mission  time  for  a  multiple 
assist  can  be  shortened  by  effecting  one  or  more  of  the 
assists  at  the  point  where  the  post-assist  orbit  makes 
its  first  crossing  of  the  assist  planet  orbit  (rather 
than  the  point  of  the  second  crossing  which  was  used  for 
each  assist  in  the  above  investigation) .  This  possibility 
is  left  for  further  study;-  however,  it  is  doubtful  that 
the  perihelia  obtained  in  these  missions  will  be  smaller 
than  those  found  above. 
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VII.  Three-Dimensional  Gravity-Assisted  Trajectories 
Purpose  and  Scope 

The  method  outlined  in  Chapter  III  for  two-dimensional 
gravity-assisted  trajectories  was  applied  to  specific 
missions  in  Chapters  IV  through  VI.  The  purpose  of  this 
chapter  is  to  outline  an  extension  of  that  method  for  the 
case  of  three-dimensional  gravity-assisted  trajectories. 

This  development  will  then  be  used  in  the  out-of-ecliptic 
mission  studies  of  Chapters  VIII  and  IX. 

Coordinate  Systems 

The  coordinate  systems  (?*ee  Fig.  43)  which  are  used 
in  the  three-dimensional  analysis  are  as  follows: 

1.  The  XjYjZj  coordinate  system  is  a  planetocentric 
system  with  its  origin  Oj  located  at  the  center  of  mass  of 
the  departure  planet  P1#  The  Xj  axis  lies  along  Vp^  (the 
velocity  of  Pj  about  the  Sun),  the  Yj  axis  points  toward 
the  Sun,  and  the  Zj  axis  forms  a  right-hand  system. 

2.  XYZ  is  an  ecliptic  coordinate  system  with  origin  0 
at  the  Sun.  The  -Z  axis  is  directed  towards  the  north-pole 

of  the  ecliptic j  the  X-Y  plane  lies  in  the  ecliptic  such  that 
the  Y  axis  points  toward  Pj  at  launch,  and  the  X  axis  forms 
a  right-hand  system. 

3.  The  pre-assist  orbit  system,  xyz,  is  an  orbital 
reference  system  located  at  the  Sun,  This  system  defines 
the  heliocentric  longitude  of  the  spacecraft  between  the 
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time  of  launch  at  P1  and  intercept  at  P2.  The  z  axis  points 
to  the  north  pole  of  the  pre-assist  orbit  plane;  the  x-y 
plane  lies  in  the  pre-assist  orbit  plane  so  that  the  y 
axis  lies  along  the  Sun-spacecraft  line,  and  x  forms  a 
right-hand  system. 

4.  The  X'Y'Z'  coordinate  system  is  another  ecliptic 
system  with  its  origin  O'  at  the  Sun.  The  X'-Y*  plane  lies 
in  the  ecliptic  so  that  the  X*  axis  points  in  the  direction 
of  a  line  drawn  from  P2  to  the  Sun  at  intercept  and  the  Y* 
axis  lies  opposite  to  the  direction  of  Vp  .  The  Z’  axis  is 
perpendicular  to  the  ecliptic  and  forms  a  right-hand  system. 

5.  The  axes  of  the  X2Y2Z2  coordinate  system  are  parallel 
to  the  respective  axes  of  the  X'Y’Z '  system,  but  the  origin 

02  lies  near  P2  at  the  point  of  intersection  of  the  inbound 
and  outbound  asymptotes  of  the  assist  hyperbola.  The  system 
moves  with  P2  about  the  Sun  and  the  X2-Y2  plane  lies  parallel 
to  the  ecliptic  (not  necessarily  in  the  ecliptic) . 

6.  The  origin  of  the  x2y2z2 system  is  also  located  at 
02,  but  the  x2-y2  plane  lies  in  the  plane  of  the  assist  orbit. 
The  y2  axis  points  opposite  to  the  V*,^  vector,  the  x2  axis 
lies  in  the  direction  of  the  x2  component  of  V »o,  and  the  z2 
axis  forms  a  right-hand  system. 

A  General  Approach  to  the  Problem 

The  problem  of  three-dimensional  gravity  assist  is  best 
introduced  with  a  discussion  of  the  underlying  motivations 
for  the  technique  presented  in  the  next  section.  The  plan 
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here  is  to  construct  two  example  problems ,  and  to  describe 
a  method  of  attack  for  each. 

For  the  first  example  (see  Fig,  44) ,  consider  an  assist 
orbit  which  lies  in  a  plane  perpendicular  to  the  ecliptic, 
such  that  the  spacecraft  enters  the  SOI  above  the  ecliptic 
(with  Vcoi  in  a  plane  parallel  to  the  ecliptic)  ,  passes  P2 
along  a  north-south  meridian,  and  exits  the  SOI  below  the 
ecliptic.  Addition  of  to  the  planet's  velocity  Vp  , 
as  in  Eq  (3-1),  gives  \flt  the  required  velocity  of  a  pre¬ 
assist  orbit  at  P2. 

From  the  corollary  presented  in  Chapter  III  it  is 
argued  that  the  pre-assist  orbit  may  be  assumed  to  enter  the 
SOI  at  any  point.  Thus,  it  is  safe  in  this  example  to  choose 
a  pre-assist  orbit  which  lies  in  the  ecliptic.  The  assump¬ 
tion  of  circular  and  coplanar  planetary  orbits  contradicts 
the  possibility  of  this  scheme *  however,  minor  velocity 
corrections  (which  have  little  influence  on  the  burnout 
conditions  at  Pj)  serve  to  establish  the  slight  pre-assist 
orbit  inclination  necessary  to  achieve  the  desired  entry 
point. 

For  a  more  general  example,  consider  the  assist  orbit 
in  Fig.  45  which  has  a  given  DOCA  and  orientation  within 
the  SOI  of  P2.  The  plane  of  the  assist  orbit  may  be  related 
to  the  ecliptic  by  three  Euler  angle  rotations  ( o ,  6 ,  0  as 
described  in  Appendix  B) ,  from  the  x2y2z2  orbit  system  to 
the  X2y2Z2  planet  system,  where  the  origin  of  these  systems 
is  located  at  the  intersection  of  the  inbound  and  outbound 
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Inbound  Asymptote 
(Parallel  to  Ecliptic) 
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asymptotes  of  the  assist  hyperbola. 

The  assist  orbit  again  defines  and  addition  of 

this  vector  to  V.,  gives  V,  which  deterxaines  the  conditions 

r  2 

of  a  pre-assist  orbit  at  P2.  In  general,  this  orbit  is 
inclined  at  some  angle  i  which  can  be  determined  from  the 
components  of  V,  in  the  XYZ  ecliptic  system.  If  this  inclin¬ 
ation  is  small  (a  fraction  of  a  degree  as  shown  in  the  next 
section) ,  the  pre-assist  orbit  may  be  assumed  to  lie  in  the 
ecliptic  like  the  previous  example,  but  if  it  is  not  small 
then  either  a  large  midcourse  correction  (broken-plane 
transfer)  or  a  180°  transfer  (from  P2  to  P2)  is  necessary. 

In  a  similar  fashion,  the  given  assist  orbit  defines 
VmQ  which  is  added  to  V1^  for  the  post-assist  velocity  V3. 
Components  of  this  vector  give  the  post-assist  orbit  inclina¬ 
tion  i3. 

The  inclination  of  planetary  orbits  serves  only  to 
enhance  the  feasibility  of  the  above  examples  when  the  pre¬ 
assist  orbit  lies  nearly  in  the  ecliptic  plane.  A  launch 
date  may  be  selected  such  that  the  spacecraft  enters  the 
SOI  at  a  time  when  P2  is  sufficiently  below  the  ecliptic. 
However,  the  assumption  of  circular  coplanar  planetary 
orbits  is  retained,  because  it  is  desired  to  investigate  the 
effects  of  gravity  assist  alone. 

Three-Dimensional  Gravity  Assist  Method 

A  three-dimensional  gravity-assisted  trajectory  is 
broken  into  the  same  conic  orbits  which  were  defined  in 
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Chapter  III:  departure  orbit,  pre-assist  orbit,  assist 
orbit,  and  post-assist  orbit. 

Departure  Orbit.  The  departure  orbit  is  defined,  as 
before,  with  the  selection  of  a  burnout  velocity  V^.  The 
choice  of  a  departure  flight  path  angle  $0  and  an  inclina¬ 
tion  angle  i  for  the  pre-assist  orbit  determines  V  ,  the 
heliocentric  velocity  at  Pl4  These  quantities  are  illustrated 
in  the  system  in  Fig.  7. 

Pre-Assist  Orbit.  A  necessary  condition  for  the  pre¬ 
assist  orbit  is  that  it  must  intercept  the  SOI  of  P2.  This 
condition  is  satisfied  for  either  of  the  following  classifi¬ 
cations  : 

1.  Category  I  Trajectories  -  the  pre-assist  orbit  for 
these  trajectories  lies  in  the  ecliptic  or  slightly  out  of 
the  ecliptic,  as  explained  below. 

2.  Category  II  Trajectories  -  the  pre-assist  orbit  for 
these  trajectories  lies  out  of  the  ecliptic  and  the  heliocen¬ 
tric  transfer  angle  612  is  180° . 

For  Category  I  trajectories  the  specified  inclination 
angle  must  be  small  enough  to  allow  the  spacecraft  to  pass 
no  further  than  a  distance  p_  above  P2 ,  where  is  the 

SOI  radius.  As  shown  in  Chapter  III,  the  SOI  radius  for 
Venus  subtends  a  heliocentric  angle  of  about  0,33°;  this 
angle  is  also  the  maximum  allowable  inclination  for  a 
transfer  of  812  «  90°.  The  maximum  inclination  increases 
as  the  transfer  angle  increases  or  decreases  from  90°. 
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For  Category  II  trajectories,  when  812  =  180°,  the 
required  inclination  angle  must  be  computed  for  given 
values  of  and  <p0.  Details  of  this  computation  are  found 
in  Appendix  C, 

The  velocity  Vx  and  flight  path  angle  at  entry  to 
the  SOI  of  P2  (r  =  rp^)  are  computed  from  the  pre-assist 
orbit  characteristics.  In  order  to  find  the  compor''nts  of 
V-!  in  the  X2Y2Z2  system,  a  coordinate  transformation  is 
made  from  the  xyz  pre-assist  orbit  system  to  tbe  XYZ  ecliptic 
system,  and  from  the  XYZ  to  the  X’Y'Z1  ecliptic  system. 
Finally,  the  components  of  Vj  in  the  X'Y'Z'  system  are  trans¬ 
lated  to  the  X2y2Z£  system. 

Assist  Orbit.  The  inbound  HEV  at  P2  is  found  by  adding 
the  components  of  ^  and  in  the  X2Y2Z2  system  with  the 
usual  relation 

VMi  =  -  Vp  (3-1) 

The  orientation  of  with  respect  to  P2  is  given  by  the 
right  ascension  a  and  declination  6  shown  in  Fig.  46. 

As  mentioned  in  the  previous  section,  three  Euler 
angles  must  be  used  to  relate  the  position  of  the  *2y2z2 
system  to  the  X2Y2Z2  system.  Two  of  these  angles  are 
found  above  (a  and  6),  but  the  third  (8)  must  be  specified 
and  is  used  as  an  arbitrary  parameter  in  the  calculations. 
This  procedure  allows  the  computation  of  all  possible 
i si-Gt  orbits  for  a  given  pre-assiot  orbit,  and  subsequently 
the  r a*-ge  of  gravicy-aesisfc  performance  for  any 
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planet  P2 . 

The  selection  of  a  DOCA,  as  for  the  two-dimensional 
case,  determines  the  turn  angle  2v  which  gives  the 
direction  of  V®0  in  the  x2 y2z2  system. 

Post-Assist  Orbit.  The  components  of  V«0  are  trans¬ 
formed  from  the  x2y2z2  to  the  X2Y2Z2  system,  and  the  helio¬ 
centric  velocity  of  the  spacecraft  as  it  leaves  the  SOI  is 

7,  =  (3-2) 

3  2  0 

The  components  of  V3  in  the  X2Y2Z2  frame  (see  Fig.  47)  are 
used  to  determine  the  flight  path  angle  <{>3  and  the  inclina¬ 
tion  i3  of  the  post-assist  orbit  with  respect  to  the  ecliptic. 

Visualization  of  Three-Dimens ional  Gravity  Assist 

The  velocity  vector  diagrams  in  Fig.  16  for  two- 
dimensional  trajectories  extend  logically  to  three-dimensional 
trajectories;  the  circular  locus  of  outbound  HEV  vectors  here 
becomes  a  spherical  locus.  A  complete  summary  of  the  velocity 
vectors  at  P2  is  shown  in  Fig,  70  of  Appendix  B.  As  in  Fig. 

16  the  diagram  may  be  broken  into  two  frames  of  reference: 
one  with  respect  to  the  Sun  (Fig.  43)  and  one  with  respect 
to  P2  (Fig.  45). 
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VIII.  Out-of-Ecliptic  Assist  at  Venus 

In  this  chapter  the  three-dimensional  method  is  applied 
to  Venus-assisted  solar  probe  trajectories.  The  purpose  here 
is  to  investigate  the  effectiveness  of  a  Venus  assist  for 
both  a  reduction  in  perihelion  and  an  increase  in  the  inclina¬ 
tion  angle.  The  trajectories  are  divided  into  the  two  clas¬ 
sifications  defined  in  Chapter  VII :  Category  I  trajectories 
and  Category  II  trajectories. 

Category  I_  Trajectories 

As  mentioned  in  the  previous  chapter,  the  pre-assist 
orbit  for  Category  I  trajectories  may  lie  no  more  than  0.33° 
out  of  the  ecliptic  for  812  =  90°.  As  812  increases  or 
decreases  from  90°  (but  does  not  approach  180°)  the  maximum 
allowable  inclination  increases  by  less  than  about  0.2°. 

It  was  found  in  the  investigation  that  a  0.1°  change 
in  pre-assist  orbit  inclination  produces  less  than  a  0,1% 
change  in  the  characteristics  of  the  post-assist  orbit. 
Therefore,  only  the  results  for  zero  inclination  angle  are 
included,  since  they  are  sufficiently  representative  of  axl 
trajectories  in  this  category. 

Also,  the  results  for  non-tangential  departure  are  not 
reported,  since  time  requirements  are  of  less  concern  here 
than  burnout  velocity  requirements. 

Computer  Program.  The  program  for  Category  I  trajec¬ 
tories  (see  Fig.  48)  includes  a  parametric  variation  of  five 
quantities  in  the  order  listed:  burnout  velocity  V^, flight 
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Fig,  48  (Continued) 
Flow  Chart  for  Program  6 
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path  angle  4*0 ,  pre-assist  orbit  inclination  i,  assist 
orbit  orientation  angle  8,  and  DOCA  d. 

The  pre-assist  orbit  is  calculated  by  the  same  tech¬ 
niques  used  in  Program  1  for  out-of-ecliptic  direct-transfer 
trajectories.  If  rT  >  pp  {Eq  B-5)  at  the  assist  planet, 
then  the  pre-assist  inclination  is  too  large  for  an  inter¬ 
cept  and  a  new  value  of  <j>0  is  selected. 

The  transformation  matrix  C,  denoted  by  E  in  Program  6, 
is  used  to  find  the  components  of  in  the  X2Y2Z2  reference 
system.  The  right  ascension  a  and  declination  6  are  calcu¬ 
lated,  and  a  value  is  selected  for  the  angle  8.  These 
angles  determine  the  transformation  matrix  D,  denoted  by  F 
in  the  program,  which  is  used  with  a  selected  DOCA  to  find 
components  of  V«o  and  rm  in  the  X2Y2Z2  system.  Finally,  the 
components  of  V3  give  i3  and  the  post-assist  orbit  character¬ 
istics. 

Notation  on  the  printed  output  is  as  follows: 

VQX,VOY ,VOZ  *  Components  of  VQ  in  X1Y1Z1  system 

V1X,V1Y,VXZ  'v.  Components  of  Vj  in  X2Y2Z2  system 
R1XP , RlYP , R1ZP  *  Components  of  rJ  in  XYZ  system 
V3X,V3Y,V3Z  *  Components  of  V3  in  X2Y£Z2  system 
DX,DY,DZ  *  Components  of  rm  in  X2Y2Z2  system  in  units 

of  P„  radii 
2 

13  *  Inclination  of  Post-Assist  Orbit 

A2  *  'Jemi-Major  axis  of  Assist  Orbit  in  P2  Radii 

V1NFXY  *  V.a  ALP  *  6 

x2y2 

PSI  -V  a  BETA  *  8 
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Mission  Profile.  The  mission  profile  in  Fig.  49  is 
designed  for  an  Atlas/Agena  launch  vehicle  and  a  200 “lb  pay- 
load  with  a  burnout  velocity  of  40,000  ft/sec. 

It  is  assumed  that  the  guidance  and  navigation  require¬ 
ments  for  missions  out  of  the  ecliptic  are  the  same  as  those 
for  missions  in  the  ecliptic,  so  that  no  additional  payload 
is  required  for  a  three-dimensional  assist. 

The  assist  at  Venus  has  tilted  the  heliocentric  orbit 
of  the  spacecraft  to  an  inclination  of  10.6°,  and  has 
increased  the  perihelion  from  0.540  AU  for  the  pre-assist 
orbit  to  0.552  AU  for  the  post-assist  orbit. 

Conclusions.  The  graph  in  Fig.  12  indicates  that  a 
perihelion  of  0.552  AU  and  an  inclination  angle  of  10.6° 
for  a  direct-transfer  orbit  requires  a  burnout  velocity  of 
a  little  more  than  43,000  ft/s er.  An  Atlas/Agena  ($7.9 
million)  was  used  in  the  mission  profile  above,  but  an 
Atlas/Centaur  ($11.9  million)  with  a  small  kick  stage  is 
required  for  the  direct-transfer.  Thus,  a  few  million 
dollars  can  be  saved  by  using  a  Venus  assist  rather  than 
a  direct  transfer  for  this  mission. 

Fig,  50  is  a  plot  of  post-assist  orbit  inclination 
angle  i3  and  perihelion  qm  for  various  values  of  burnout 
velocity  Vfe  (solid  lines)  and  orientation  angle  S  (dashed 
lines) .  These  curves  are  plotted  directly  from  computer 
output  for  a  DOCA  of  100  NM.  The  dotted  line  on  the  right- 
hand  side  of  the  graph  is  a  locus  of  maximum  i3  for  any  Vj-,. 
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This  line  indicates  that  the  highest  inclination  angles  are 
achieved  with  either  low  or  high  burnout  velocities  (approx¬ 
imately  40,000  or  70,000  ft/sec  respectively).  It  can  be 
concluded  from  this  graph  that  a  Venus  assist  produces,  at 
best,  an  inclination  of  about  10°. 

Pig.  50  also  indicates  that  a  large  gain  in  i3  reflects 
a  substantial  increase  in  perihelion.  For  instance,  at 
Vb  =  40,000  ft/sec,  perihelion  increases  slightly  from  0.337 
AU  to  0.357  AU  for  a  change  in  i3  from  0°  to  5°,  respectively; 
however,  the  perihelion  increases  rapidly  to  0.550  AU  for 
i3  =  10.55c . 

A  plot  of  maximum  inclination  angle  versus  burnout 
velocity  is  given  in  Fig.  51  for  various  DOCA  values.  These 
•curves  are  obtained  from  a  plot  of  3  versus  i3  for  each 
value  of  Vb  and  d;  the  maximum  value  of  i3  is  then  transferred 
to  Fig.  51  and  a  locus  is  obtained  for  each  value  of  d.  The 
double  maximum  observed  in  Fig.  50  is  again  evident  in  Fig. 

51  for  each  value  of  d. 

The  perihelion  for  maximum-inclination  missions  and 
ecliptic  direct-transfer  missions  (dashed  line)  is  compared 
in  Fig.  52.  From  these  curves  it  is  found  that  perihtlion 
may  be  significantly  reduced  for  the  lower  range  of  V^.  For 
instance,  at  =  41,000  ft/sec  and  d  =  100  NM,  perihelion 
is  reduced  from  0.71  AU  for  direct  transfer  (10°  inclination) 
to  0.50  AU  (10.45°  inclination)  for  Venus  assist. 

Figs.  51  and  52  also  indicate  that  variations  in  L'OCA 
below  1000  NM  do  not  strongly  affect  the  post-assist  orbit. 
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A  change  in  d  from  0  to  10  0  NM  decreases ’.the  inclination  by 
about  2°  and  increases  the  perihelion  by  less  than  0.01  AU. 

A  comparison  of  Fig.  53  with  Fig.  21  indicates  that  the 
mission  time  for  Category  I  trajectories  is  about  the  same 
as  for  single  Venus-assist  missions  in  the  ecliptic. 

In  summary,  Venus-assisted  Category  I  trajectories  are 
useful  for  10°  out-of-ecliptic  solar  probe  missions  with 
perihelia  of  about  0.3  AU.  The  best  gains  for  these 
trajectories  (decrease  in  perihelion  and  increase  in  inclina¬ 
tion)  are  realized  for  the  lower  range  of  burnout  velocities 
around  40,000  ft/sec. 
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Fig.  52  (ft/sac  x  10 3 ) 

Perihelion  for  Maximum  i,  as  a  Function  of  Burnout  Velocity 


Maximum 
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Category  IX  trajectories 

The  analysis  for  Category  II  trajectories  is  identical 
to  that  used  for  Category  I  trajectories ,  except  that  the 
pre-assist  orbit  inclination  must  now  be  computed,  rather 
than  specified,  for  the  given  values  of  burnout  velocity  Vb 
and  ddparture  flight  path  angle  <j>0.  This  inclination  angle 
is  plotted  in  Fig.  54  as  a  function  of  for  a  tangential 
departure;  derivation  of  the  equation  for  this  curve  is 
given  in  Appendix  C. 

As  defined  in  the  previous  chapter,  the  heliocentric 
transfer  angle  e12  for  these  missions  is  180°.  For  a  given 
value  of  <p0,  the  pre-assist  orbit  inclination  increases 
with  Vfo;  therefore,  any  additional  burnout  velocity  above 
the  minimum  Vb  for  a  Venus  intercept  (Hohmann  transfer  which 
requires  37,700  ft/sec)  is  used  to  tilt  the  pre-assist  orbit 
out  of  the  ecliptic. 

Computer  Program.  The  program  for  Category  II  trajecto¬ 
ries  (Program  7  in  Appendix  E)  is  identical  to  Program  6, 
except  that  the  pre-assist  orbit  inclination  is  calculated 
for  each  value  of  and  A  flow  chart  for  this  program 
can  be  obtained  by  replacing  the  i  loop  in  Fig.  48  with  a 
box  for  the  calculation  of  i.  The  same  transformations  are 
made  and  the  same  notation  is  used  for  both  programs. 

Mission  Profile.  The  Atlas/Centaur/TE-364-3  launch 
vehicle  is  used  for  the  mission  profile  in  Fig.  55.  From 
Fig.  4  the  burnout  velocity  for  this  vehicle  for  a  250 -lb 
payload  is  53,000  ft/sec.  This  velocity  provides  a  pre- 
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assist  orbit  inclination  of  23.0°,  and  the  assist  at  Venus 
increases  this  inclination  to  24,9°. 

The  perihelion  of  the  pre-assist  orbit,  since  this  is 
a  cotangential  transfer,  is  necessarily  equal  to  the  mean 
orbital  radius  of  Venus  (0.723  AU) .  After  the  assist  this 
perihelion  is  reduced  to  0.529  AU. 

Conclusions.  In  the  mission  profile  above,  it  was 
found  that  a  perihelion  of  0.529  AU  and  an  inclination  of 
24.9°  can  be  achieved  with  an  Atlas/Centaur/TE-364-3  ($16 
million)  by  using  a  Venus  assist.  Fig.  12  indicates  that  a 
burnout  velocity  of  55,700  ft/sec  is  necessary  to  attain  the 
same  perihelion  and  inclination  with  a  direct  transfer? 
thus,  a  Saturn  I  ($29.8  million)  is  required  for  this  mission. 

Fig.  56,  similar  to  Fig.  50,  is  a  plot  of  perihelion 
versus  post-assist  orbit  inclination  for  various  values  of 
burnout  velocity  Vb  (solid  lines)  and  orientation  angle  B 
(dashed  lines) .  The  heavy  line  which  divides  the  right  and 
left  sides  of  the  plot  represents  a  V^  of  51,322  ft/sec. 

This  particular  burnout  velocity  and  its  associated  pre-assist 
inclination  combine  to  give 


That  is,  the  component  of  Vj  in  the  direction  of  motion  of 
V&nus  is  equal  to  the  velocity  of  Venus  about  the  Sun, 

Burnout  velocities  to  the  left  of  the  heavy  line  give 
Vx  >  vPf,  so  that  for  $0  ®  90°  the  right  ascension  o  of 
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the  inbound  HEV  vector  is  0°;  that  is,  the  spacecraft  enters 
the  SOI  behind  the  planet  and  exits  the  SOI  in  front  of  the 
planet.  Likewise,  the  right  half  of  the  plot  is  for  burnout 
velocities  which  give  <  Vp^ ,  and  o  =  180°  so  that  the 

spacecraft  enters  the  SOI  in  front  of  P2  and  exits  to  the 
rear. 

Thus,  the  angle  a  switches  from  0°  to  180°  along  the 
dividing  line.  This  transition  must  be  accompanied  by  a 
change  in  sign  of  the  orientation  angle  8,  so  that  aphelion, 
rather  than  perihelion,  is  achieved  at  P2. 

Fig.  56  indicates  that  the  minimum  perihelion  and  maxi¬ 
mum  inclination  angle  for  a  given  occurs  for  6  =  +90°  on 
the  left  side  of  the  plot  and  8  ®  -90°  on  the  right  side. 
Thus,  8  =  ±90°  represents,  the  most  desirable  Category  II 
trajectories,  and  further  results  of  these  missions  are 
presented  in  Figs.  57  and  58. 

Fig.  57  is  a  plot  of  post-assist  orbit  inclination  i3 
as  a  function  of  for  various  values  of  DOCA,  A  plot  of 
pre-assist  orbit  inclination  (dashed  line)  from  Fig.  54  is 
superimposed  on  the  graph  so  that  i3  may  be  compared  with 
inclination  before  the  assist.  It  is  found  that  gravity 
assist  does  not  significantly  increase  the  inclination  for 
Category  II  trajectories..  The  greatest  increase,  for  the 
higher  range  of  V is  only  about  3  to  4  degrees. 

The  perihelion  for  Category  II  trajectories  is  shown 
in  Fig.  58  as  a  function  of  Vb,  Significant  reductions  in 
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perihelia  (below  the  0,723  AU  perihelion  of  pre-assist 
orbits)  can  be  obtained  in  these  missions.  The  smallest 
perihelion  is  0,439  AU  at  a  burnout  velocity  of  41,500 
ft/sec. 

Mission  times,  shown  in  Fig,  59,  are  slightly  longer 
for  the  Category  II  trajectories  than  for  direct-transfer 
trajectories,  but  are  well  within  reason  for  solar  probe 
missions. 

In  conclusion,  Venus-assisted  Category  II  trajectories 
are  useful  for  20 °-  to  30®-out-of-ecliptic  missions  with 
perihelia  of  about  0.45  to  0.55  AU.  A  Venus  assist  offers 
little  improvement  in  pre-assist  inclination,  but  serves 
to  reduce  the  perihelion  by  as  much  as  0.28  AU, 
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IX,  Jupifcer-Assist  Missions 

■  hi  ■  i  ii  1 1 1  iiiiiiii  wmui  w—  m mm  mi  c— —— — m  ifunw 

As  stated  in  the  introduction,  Jupiter-assisted  trajec¬ 
tories  have  been  thoroughly  investigated  in  the  literature. 
However,  a  brief  study  of  these  missions  is  conducted  here 
so  that  the  results  may  be  directly  compared  with  other 
missions  in  the  thesis. 

Only  Category  I  trajectories  are  used  for  Jupiter 
missions  because  of  the  high  launch  energy  required  for 
an  out-of-ecliptic  pre-aasist  orbit. 

Computer  Programming 

The  program  for  Jupiter-assisted  missions  is  identical 
to  Program  6,  except  that  a  few  modifications  are  necessary 
to  account  for  the  fact  that  the  pre-assist  orbit  lies  out¬ 
side,  rather  than  inside,  the  Barth's  orbit.  The  first  of 
these  modifications  is  a  sign  change  in  Eq  (2-4)  for  the 
heliocentric  departure  velocity  V0. 

Eq  (2-18),  for  the  true  anomaly  at  Earth  launch,  must 
be  modified  to  read  as  follows: 

£0  «  arc  tan  (tan  f0),  if  tan  f0  >  0 

(10-1) 

f0  «  arc  tan  (tan  f0)  +  w,  if  tan  f0  <  0 

A  similar  change  is  required  for  ilt  the  true  anomaly  at  P2* 

Many  of  the  post-assist  orbits  for  Jupiter-assist  missions 
are  hyperbolic  (solar  system  escape  with  Jupiter  assist  is 
investigated  in  Refs  11  and  13) .  These  trajectories  are 
avoided  as  follows:  if  the  quantity  rp2  v|  /m0  is  less  than 
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2  (Ref  21; 72),  the  post-assist  orbit  is  elliptic,  and  the 
trajectory  results  are  printed;  otherwise,  the  traje  tory  is 
hyperbolic,  and  another  set  of  parameters  is  selected. 

Conclusions 

Pig.  60  is  a  plot  of  perihelion  versus  burnout  velocity 
for  a  Jupiter  assist  in  the  ecliptic  (i3  «  0  »  0°).  Solar 
impact  can  be  achieved  with  burnout  velocities  between  51,000 
and  53,000  ft/sec  for  a  DOCA  of  500,000  NM  (13  Jupiter  radii). 
Smaller  DOCA  values  can  be  used  to  provide  retrograde  post¬ 
assist  orbits  with  small  perihelia. 

Out-of-ecliptic  performance  of  Jupiter-assist  missions 
is  shown  in  Fig.  61,  where  perihelion  qm  and  post-assist 
inclination  i3  are  plotted  (with  DOCA  «  500,000  NM)  for 
various  values  of  orientation  angle  6  (dashed  lines)  and 
burnout  velocity  Vb  (solid  lines) »  From  this  graph  it  is 
found  that  90°-out-of-ecliptic  orbits  with  perihelia  leas 
than  0.1  AU  can  be  achieved  with  burnout  velocities  above 
52,000  ft/sec.  The  best  solar  probe  missions  are  attained 
when  3  is  about  10°. 

A  summary  of  results  for  90°-out-of-ecliptic  missions 
is  given  in.  Table  III,  These  missions  are  obtained  from 
a  plot  of  i3  versus  DOCA  for  each  burnout  velocity  at 

3  «  10°,  The  DCCA  for  i3  «  90°  on  this  plot  (second 
column  in  Table  III)  is  used  to  interpolate  the  data  on 
the  computer  printout  for  qm  and  Top. 

The  values  of  TOP  and  Vb  in  Table  III  are  plotted  on 
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the  graph  Pig.  62,  Also  included  on  the  graph  is  a 
plot  of  mission  time  for  trajectories  in  the  ecliptic  with 
a  DOCA  of  5;‘0,000  NM.  Little  difference  is  found  between 
the  time  requirements  for  Jupiter-asuist  missions  in  or 
out  of  the  ecliptic  plane. 

In  summary,  it  should  be  noted  that  Jupiter-assisted 
trajectories  can  be  achieved  with  present-day  launch  vehicles. 
An  Atlas/Centaur/TE-364-3  is  capable  of  providing  a  burnout 
velocity  of  52,000  ft/see  with  a  300-lb  payload.  This  launch 
velocity,  combined  with  a  Jupiter  assist,  can  produce  a 
perihelion  of  about  0.05  AU  and  an  inclination  of  90°. 

The  disadvantage  of  Jupiter  assist,  however,  is  the 
three-year  mission  time  involved?  also,  depending  on  orbital 
inclinations,  the  spacecraft  may  be  required  to  pass  through 
the  asteroid  belt  which  significantly  reduces  the  reliability 
of  these  missions. 
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Jupiter-Assist  Missions 
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Table  III 


90 9 -OUT-OP-ECLIPTIC  JUPITER-ASSISTED  TRAJECTORIES 


Vb  (ft/sec) 

DOCA  (NM) 

q  (AU) 

TQp  (years) 

51,000 

380,000 

0.059 

3.17 

52,000 

500,0.00 

0.048 

2.73 

53,000 

530,000 

0.056 

2.51 

54,000 

540,000 

0.044 

2.36 

3 
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X.  Conclusions  and  Recommendations 

Solar  Probe  Missions  in  thfe  Ecliptic 

Direct-transfer  trajectories,  characterized  by  high 
reliability  and  launch-date  insensitivity,  are  preferred 
for  solar  probe  missions  to  about  0.5  AU  in  the  ecliptic. 
However,  the  burnout  velocities  required  for  smaller 
perihelia  make  it  necessary  to  use  expensive  launch 
vehicles  such  as  th-  Saturn  I  or  Saturn  V. 

It  has  been  shown  in  the  thesis  that  a  gravity 
assist  at  Venus  can  be  used  to  significantly  reduce  the 
perihelion  of  solar  probe  trajectories.  A  few  extra 
pounds  of  planetary-approach  guidance  equipment  act  as  a 
substitute  for  large  and  expensive  booster  vehicles. 

The  multiple-assist  missions  at  Venus  offer  an  even 
greater  savings  in  launch  vehicle  costs.  The  Atlas/ 
Centaur/TE-364-3  ($16  million)  can  achieve  smaller  perihelia 
with  multiple  assist  than  a  Saturn  V  ($125  million)  with 
direct  transfer.  An  additional  advantage  of  these  missions, 
despite  their  long  time  requirements,  is  that  solar  data 
can  be  collected  over  several  regions  of  space  with  a 
single  launch. 

It  is  recommended  that  planetary-approach  guidance 
techniques  and  space-vehicle  reliability  be  improved  for 
an  attempted  triple-assist  mission.  The  probability  of 
accomplishing  the  first  passage  at  Venus  is  near  certainty 
with  present-day  technology.  If  the  second  or  third 
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passages  fail,  the  mission  can  be  classed  as  a  partial 
success,  for  data  will  have  been  collected  over  at  least 
one  solar  region. 

Multiple-assist  missions  are  the  most  reliable  of  any 
of  the  gravity-assist  missions  in  the  thesis.  Trajectories 
can  be  designed,  as  shown  in  the  double-assist  mission 
profile,  such  that  the  pre-assist  orbit  period  is  equal  to 
the  period  of  Venus'  orbit.  This  provides  the  opportunity 
for  further  attempts  to  pass  Venus,  if  one  should  fail. 

The  Mercury-Venus  combination-assisted  trajectories 
were  found  to  be  more  trouble  than  they  are  worth.  It  is 
also  concluded  from  that  investigation  that  Mercury  is 
of  little  value  as  an  assist  planet. 

Jupiter  offers  the  best  performance  for  gravity-assisted 
trajectories;  in  fact,  retrograde  orbits  can  be  obtained. 
However,  the  time  requirements,  communication  difficulties 
during  the  flyby,  and  penetration  of  the  asteroid  belt  raeke 
these  missions  prohibitive.  In  spite  of  these  drawbacks, 
a  Jupiter  assist  is  necessary  to  achieve  perihelia  less 
than  about  0.05  AU  with  current  launch  vehicles. 

Out-of-Ecliptic  Solar  Probe  Missions 

The  launch  energies  for  direct-transfer  trajectories 
more  than  20°  out  of  the  ecliptic  require  the  use  of  the 
Saturn  I  or  Saturn  V.  It  has  been  shown  that  Venus-assisted 
Category  1  and  II  trajectories  are  useful  in  reducing  this 
launch  requirement. 
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•  The  Category  I  trajectories  are  recommended  for  10 °- 
out-of-ecliptic  missions.  The  perihelia  for  these  trajec¬ 
tories  is  about  0.2  AU  less  than  that  obtained  with  direct- 
transfer. 

It  is  suggested  that  Category  II  trajectories  be  used 
l?r  20°-  to  30°-out-of-ecliptic  missions.  Most  of  the 
launch  energy  for  these  missions  is  devoted  to  the  required 
pre-assist  inclination ,  but  a  Venus  assist  serves  to  reduce 
the  perihelion  by  as  much  as  0.28  AU. 

A  Jupiter  assist  is  necessary  to  achieve  inclination 
angles  above  30°.  Direct-transfer  trajectories  with  these 
inclinations  require  launch  velocities  much  greater  than 
present-day  technology  can  provide. 

The  attractiveness  of  the  90®-out-of-ecliptic  trajec¬ 
tories  obtained  with  a  Jupiter  assist  cannot  be  overlooked. 
These  missions  are  recommended  for  future  consideration 
when  spacecraft  can  be  designed  to  withstand  a  three-year 
mission  time ,  and  when  planetary- approach  guidance  techniques 
have  become  more  perfected. 
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Appendix  A 

Two-Dimensional  Gravity- Assist  Equations 
Departure  Orbit 

The  burnout  velocity  on  a  direct-ascent  hyperbola  at 
the  departure  planet  Pj  is  given  by 

2Ppi 

V?  =  v£p  +  - -  (2-4) 

b  Pl  RP! 


Selection  of 
HEV  at  P  becomes 


l 


After  choosing  a  value  for  V^,  the 


/ 


2"M 


1/2 


(A-l) 


Pre-Assist  Orbit 

Characteristics  of  the  pre-assist  orbit  are  calculated 
with  the  same  equatiohs  which  were  used  for  the  heliocentric 
transfer  orbit  in  Chapter  II.  However,  the  initial  steps 
for  finding  70  are  more  direct,  since  the  orbit  is  now  two- 
dimensionax. 

Selection  of  6.  A  departure  flight  path  angle  $0  is 

MM  MM 

selected,  and  the  Law  of  Sines  is  used  to  determine  the  angle 
n1  shown  in  Fig.  63: 

Vpj  cos 

sin  n1  *»  . (A-2) 

V 

v«p 


\ 
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Calculation  of  VQ.  Alsc^  from  Fig.  63 


n  »  it  -  arc  sin  (sin  nl)  (A-3) 


and  a  =  n/2  -  n  +  $  (A-4) 

0 


Finally,  from  the  Law  of  Sines, 


V  a  v  „  SAfi-JL 
o  »p , 

l  cos  $ 


(A-5) 


Pre-Assist  Orbit  Characteristics . 
The  equations  for  pre-assist  orbit 
characteristics  derived  in  Chapter  II, 
are  listed  here  for  convenience: 


V2  /y  cos  $  sin  <p 


tan  f0  = 


'Px  'o'  r©/ - To  - To 

1  -  K,  Vi  /“o) 


(2-10) 


'Pi 


al  ~ 


2  "  rPi  Vo  /vo 


(2-11) 


el  **  (rp j  Vq  /»©  "  X)  2  sin2$0  +cos2^o  (2-12) 


A  necessary  constraint  on  the  pre-assist  orbit  is  that 
its  perihelion  must  be  at  least  as  small  as  the  assist  planet 
aphelion.  According  to  the  assumption  of  circular  and 
coplanar  planetary  orbits,  with  the  exception  cf  Mercury,  the 
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assist  planet  aphelion  is  taken  as  the  mean  distance  of  the 
planet  from  the  Sun.  This  constraint  may  be  expressed  as 

r?2  aj  (1  -  ej)  (A-6) 

where  rp^  is  the  mean  distance  of  the  assist  planet  from  the 
Sun  (aphelion  of  Mercury,  if  Mercury  is  used  as  the  assist 
plane -)  and  ax  (1  -  ex)  is  an  expression  for  the  perihelion 
of  the  pre-assist  orbit. 

Pre-Assist  Orbit  Conditions  at  the  Assist  Planet.  As 
explained  in  Chapter  III,  the  pre-assist  orbit  is  assumed  to 
end  at  the  point  where  it  first  intersects  the  orbit  of  the 
assist  planet.  Thus,  the  heliocentric  radius  of  the  pre-assist 
orbit  at  the  time  of  intercept  is  given  by 

r  =  rp2  (A-7) 

Eqs  (A-7)  and  (2-11)  are  used  to  compute  the  spacecraft  veloc¬ 
ity  with  respect  to  the  Sun  at  intercepts 


V? 


(A-8) 


from  the  vis-viva  integral.  Also,  for  the  flight  path  angle 
(Ref  16:83)  , 


(A-9) 
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The  true  anomaly,  using  Eqs  (A-8)  and  (A-9) ,  is  given  by 


(rP2  vi  /»»0)  cos  sin  ^ 
1  *  (rP2  v?  ^0)sinZ^i 


(A-10) 


where-,  f-j  =  arc  tan  (tan  fx)  (A-ll) 

The  heliocentric  transfer  angle  from  P,  at  launch  to  P 

l  2 

at  intercept  is 


«12  -  *i  “  *0  <A-12) 

and  the  transfer  time  is 

t  (a^,  pQ,  e^,  fj)  -  t  (a^,  ©j,  fQ)  (A-13) 


Assist  Orbit 

Fig.  64  which  was  discussed  in  Chapter  III  is  used  to 
write  a  vector  expression  for  the  inbound  HEV  of  the  space¬ 
craft  with  respect  to  the  assist  planet 


V  B  V  ••  V 
•i  1  p, 


(3-1) 


The  magnitude  of  this  velocity  is  found  by  using  the  Law  of 
Cosiness 


v°*i 


+  V*  -  2V,  V-  sin  $ 


*2 


P2 


(A-14) 


An  expression  for  the  HEV  on  a  hyperbolic  orbit  (Ref 
16s 96)  is  used  to  find  the  semi-major  axis  of  the  assist 
hyperbolas 
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Now  an  additional  parameter 
must  be  selected  in  order  to 
completely  determine  the  geometry 
of  the  assist  hyperbola  -  the 
DQCA  which  is  denoted  by  d  in 
the  mathematical  expressions. 

Selection  of  d.  The 
distance  between  the  vertex  and 
focus  of  the  assist  hyperbola  is 
given  by 

rm  °  a2  ^e2  "  D  (A-16) 


from  orbital  mechanics  (Ref  2:152).  The  eccentricity,  e2, 
is  undetermined  in  Eq  (A-16) ,  but  can  be  found  by  assigning 
a  value  to  d,  where 

r_  “  R_  +  d  (A-17) 

m  i?2 

Eqs  (A-16)  and  (A-17)  are  then  combined  as  follows: 


e2  “  ~  +  a2  +  d)  (A-18) 

a2 

The  angle  v  between  the  hyperbola  asymptote  and  conjugate  axis 
is  given  by 


sin  v  *  i- 
ez 


(A-19) 
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As  explained  in  Chapter  III,  the  effect  of  gravity 
assist  is  tc  rotate  the  inbound  HEV  vector  V„^  through  an 
angle  2v  in  the  plane  of  motion.  The  direction  of  this 
rotation  depends  upon  the  side  of  P2  which  is  passed  by  the 
spacecraft.  In  this  report,  since  it  is  desired  to  slow  the 
spacecraft  with  respect  to  the  Sun,  is  rotated  counter¬ 
clockwise  into  the  vector  7«a  as  shown  in  Fig.  64. 

For  navigation  purposes,  it  is  convenient  to  use  the 
impact  parameter  or  aiming  vector  which  is  the  perpendicular 
distance  from  a  point  on  the  hyperbolic  asymptote  to  the 
center  of  mass  of  P2.  The  distance  is  given  by  (Ref  2:152) 


The  impact  parameter  has  not  been  used  in  the  computations, 
since  this  study  is  not  concerned  with  navigation,  but  it  is 
given  here  to  show  that  it  can  easily  be  found  from  quantities 
which  have  been  derived  thus  far, 

Post-As3ist  Orbit 

The  post-assist  orbit  is  assumed  to  originate  at  the 
end  point  of  the  pre-assist  orbit,  so  the  heliocentric  radius 
at  intercept  is  again  given  by  Eq  (A-7) .  Characteristics  of 
this  orbit  are  then  computed  from  the  post-assist  velocity  V3, 
and  flight-path  angle  $3, 

Post-Assist  Orbit  Conditions  at  the  Assist  Planet.  The 
heliocentric  velocity  W3  is  determined  from  Fig.  64  as  follows* 
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V3  «  Vp2  +  V„o 


(3-2) 


The  angle  in  Fig.  64  may  be  used  to  find  the  magnitude  of 
V3.  From  the  Law  of  Sines, 

V, 


sin  (2v  +  tji)  =  — cos 
V„. 

i 


and 


«=  arc  sin 


vi 

_ COS  tfr, 

L*"i 


-  2v 


(A-21) 


(A-22) 


The  Law  of  Cosines  gives 


V|  «=  V|  +  V2  -  2Vp  V,*  cos  t 
2  0 


'P2  V“o 


where,  of  course 


(A-23) 


V«o  e  Vcoi 
0  1 


(A-24) 


The  direction  of  V3  is  specified  by  the  flight  path  angle 


cos 


*3 


S3 


sin 


(A-25) 


The  true  anomaly  of  the  post-assist  orbit  at  P2  is 
given  by 


rP2  vf  /'■»)  °°«  *>  Bln  *3  (A_26) 

1  -  (rP2  V!  /»*)  sil>2*3 

Post-Assist  Orbit  Characteristics .  The  orbit  conditions 
in  the  previous  section  are  used  in  the  following  equations 
for  the  post-assist  orbit  characteristics! 


tan  f 3  » 
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rP 

a.  s  . . . . . —  (b-27) 

2  -  fc2  V3  /».) 

e|  »  (rp  V|  /pQ  -  l) 2  sin2$3  +  cos2$3  (A-28) 

The  totwi  time  from  launch  to  perihelion  is 

T0P  “  Ti2  “  t  <a3#  e3,  f3)  (A-29 ) 

and  the  aphelion  and  perihelion  of  the  post-assist  orbit  are, 
respectively, 

ra^  *  a3  (1  +  e3)  (A-30) 


*  a3  (1  -  e3) 


(A-31) 
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< 

Appendix  B 


Three-Dimensional  Gravity-Assist  Equations 


Departure  Orbit 

Selection  of  Vb.  The  HEV  at  the  departure  planet  P  is 
found  with  the  same  expression  used  in  the  two-dimensional 
trajectories : 


(A-l) 


for  a  given  value  of  V^. 

Pre-Assist  Orbit 

Characteristics  of  the  pre-assist  orbit  ?,re  found  with 
the  heliocentric-transfer  orbit  equations  of  Chapter  II. 

These  equations  are  listed  again  for  completeness  of  the 
derivation. 

Selection  of  $0f  i.  The  XjY* Zj coordinate  system,  shown 
in  Fig,  1 ,  is  used  to  write  an  equation  for  the  heliocentric 
velocity  at  launch} 

V0  =  Vp^  sin  $0  cos  i  - 

sj  sin24>0  cos2i  +  (2-9) 

Pre-Assist  Orbit  Characteristics .  The  geometric  proper¬ 
ties  of  the  pre-assist  orbit  are  computed  as  outlined  in 
Chapter  II  and  Appendix  A: 


> 

i 
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(rP,  vl  /t*o)  c°s  ♦  „  Sin  *„ 
tan  f0  a  - j - - - 

1  -  (rP,  Vo  /'■J  sin  ^0 


(2-10) 


a.  =  - Hi - 

2  -  (rP,  vo  /•*«) 

®l  •  (rPj  /"o  ■  X)2  sin2*o  +  oos2*o 
with  the  constraint  that 

rp^  ^  a1  (1  -  ej)  (A-9) 

Pre-Assist  Orbit  Conditions  at  the  Assist  Planet.  As 
explained  in  Chapter  VII,  the  pre-assist  orbit  can  intersect 
the  SOI  of  P2  for  only  two  cases:  inclination  angle  small 
or  heliocentric  transfer  angle  0i2  equal  to  180°.  For 
either  of  these  cases,  rx ,  the  distance  of  the  spacecraft 
from  the  Sun  at  the  time  it  enters  the  SOI,  is  assumed  to 
be  rp^,  the  moan  orbital  radius  of  the  assist  planet.  Under 
this  assumption  the  heliocentric  conditions  of  the  pre-assist 
orbit  at  P2  may  be  expressed  as  in  Appendix  A: 


(2-11) 

(2-12) 


v? 


sin  $x 


(1  -  ef) 


|fP2  (2al  “  rP2) 


1/2 


(A-8) 


(A-9) 
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(rp2  V1  /wo)  C0S  *1  sin  *1 
.  tan  f,  =  — — — . — - - -  (A-10) 

1  “  {rP2  vi  /**©)  sin2*i 

The  transfer  angle  and  time  on  the  pre-assist  orbit  is 
given  by  the  following  expressions* 

e12  «  -  f0  (A-12) 

T12  =  t  (al#  elf  f:)  -  t  (a1#  n0,  ,  f0)  (A-13) 

XYZ  Components  of  r l  and  ^ .  The  vector  components  of 
Fj  and  Vl  in  the  XYZ  system  are  found  with  the  following 
coordinate  transformation  (see  Figs.  43  and  65) : 

1,  The  xyz  orbital  reference  system,  initially  aligned 
with  the  XYZ  ecliptic  system,  is  rotated  through  the  angle  i 
about  the  line  of  nodes  (the  Y  axis) .  This  places  the  x-y 
plane  in  the  pre-assist  orbit  plane. 

2.  The  xyz  system  is  rotated  through  ;n  angle  012 
about  the  z  axis,  so  that  the  y  axis  points  toward  the  assist 
planet  P2  at  the  time  of  intercept.  These  rotations  are 
expressed  in  the  following  matrix  notation?  where  for  conven¬ 
ience,  e12  =  0,  s  represents  sine,  and  c  represents  cosine: 


(  'I 
X 

/ 

C0 

ci 

S0 

—  C0 

si 

/  n 

X 

<y 

>  =  < 

— s  0 

ci 

C  9 

80 

si 

<J  Y 

z 

\  / 

si 

0 

ci 

j 

Z 

N  / 

Since  the  transformation  is  orthogonal,  the  inverse  of  the 
transformation  matrix  is  given  by  its  transpose. 
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Thus, 


/  N 

f  N 

X 

X 

J  Y 

>«  A  < 

y  1 

Z 

z 

(B-2) 


where  tlie  matrix  A  is: 


c9 

ci 

~s9 

ci 

A  «  < 

S0 

C9 

-ce 

si 

S8 

si 

Components  of  r’j  in  the  XYZ  system 


✓  N 

f  'l 

/  \ 

0 

0 

<r'Y 

>  ~  A< 

r‘ 

>s  A< 

rp  > 
*2 

.■v 

o 

0 

si 

< 

o 

ci 

are 


\ 


> 


given  by 


(B-3) 


(B-4) 


The  following  constraint  must  be  satisfied  in  Eg  (B-4) 

r>  $  pp  (B-5) 

2  r2 

that  is,  the  distance  of  the  spacecraft  above  the  ecliptic 
plane  must  not  exceed  pp^,  the  radius  of  the  SOI  of  Pz,  If 
this  constraint  is  not  satisfied  then  the  inclination  i  of 
the  pre-assist  orbit  plane  is  too  large  for  an  intercept, 
and  a  smaller  value  must  be  selected. 

The  velocity  ^ ,  as  shown  in  Pig.  65,  has  components 
in  the  negative  x  and  negative  y  directions,  so  its  components 
in  the  XYZ  system  are 
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✓  N 

/  s 

V1 

-Vj  sin  <j»1 

X 

<kl 

>  =  A  < 

-Vx  cos  <^1 

k. 

0 

Transformation  of  Vx  from  XYZ  to  X2Y2Z2.  In  order  to 
define  the  relationship  of  Vx  with  other  vectors  at  P2  it 
is  necessary  to  transform  it  from  the  ecliptic  system  XYZ 
to  the  assist  planet  system  X2Y2Z2.  This  transformation, 
shown  in  Figs.  66  and  67,  is  accomplished  as  follows* 

1.  The  X'Y'Z'  and  XYZ  ecliptic  systems  are  initially 
aligned  at  the  Sun.  The  X'Y'Z'  axes  are  rotated  about  the 
Z  axis  through  an  angle  n,  where 

n  =  y  “  */2  (B-7) 

s  e12  -  »/2  (B-8) 

so  that  the  X'  axis  is  directed  along  a  line  from  P2  to  the 
Sun,  and  the  Y'  axis  is  directed  opposite  to  Vp^,  the  velocity 
of  P2. 

2.  The  X'Y'Z'  system  is  translated  from  the  Sun  to  a 
position  02  near  P2,  where  the  inbound  and  outbound  assist 
hyperbola  asymptotes  intersect.  This  position  is  not  defined 
until  the  DOCA  and  an  angle  6,  discussed  later,  are  selected; 
for  the  present  degree  of  approximation,  it  is  sufficient  to 
know  that  02  lies  within  the  SOI  of  P2. 

The  relationship  between  the  XYZ  and  X'Y'Z*  systems  is 

given  by 
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✓  N 

✓  N 

X* 

X 

Y' 

>  =*  B< 

Y  l 

Z* 

Z 

N.  ✓ 

N  ✓ 

where  the  transformation  matrix  B  is  defined  as 


cn 

sn 

0 

B  » < 

-sn 

cn 

0  r 

0 

0 

l 

N 

✓ 

(B-9) 


(B-10) 


The  translation  from  X'Y'Z' 
so 


✓  N 

X, 


<  'i. 


X* 

>=<:  y« 


to  X2Y2Z2  involves  no  rotation. 


(B-ll) 


Finally,  the  vector  in  the  X2Y2Z2  system  near  the 
assist  planet  is 


'  V 

fvi  1 

*X 

lx2 

% 

lz 

>=  B  < 

VlY  ’ 

V1 

k  l7.  ' 

n  c  < 


-Vj  sin  «f»  j 
-Vj  cos 


(B-12) 


(B-13) 


where  use  has  been  made  of  Eq  (B-6)  and  the  definition: 
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C 


BA 

✓ 

(  cn  c8  ci  +  sn  se) 
i  (-sn  ce  ci  +  cn  se) 
(-c6  si) 


(-cn  s8  ci  +  sn  c0) 
(  sn  s8  ci  +  cn  c0) 
(  s0  si) 


(B-14)- 


(  cn  si) 
(-sn  si) l 


(B-15) 


(Ci)  J 


Assist  Orbit 

Calculation  of  Vo^.  The  inbound  HEV  at  P2  is  shown  in 
the  vector  diagram  of  Pig.  46.  The  vectors  ,  7ir  and  V*,^ 
are  translated  and  added  within  the  X2Y2Z2  reference  frame; 
the  fact  that  the  position  of  02  has  not  yet  been  defined 
does  not  affect  the  results,  since  vectors  may  be  translated 
to  any  point  in  space  with  no  change  in  magnitude  or  direction. 
Thus, 

V^i  =  -  Vp2  (3-1) 

In  matrix  notation,  Sq  (3-1)  becomes 


/  > 
v«i 

X2 

✓  > 
V1X 

/  N 

0 

<  v«i 
ly 

2 

>  =  < 

V1Y 
*  2 

>  -  < 

”VP  ? 
f2 

v„t 

k  H  J 

o 

V*Z 

k  *z) 

0 

2 


which  gives  the  magnitude  of  the  vector, 

1/2 


v»i  = 


v»4  +  vii  +  v®  j 

Xx2  *2  % 


(B-17) 


Also,  Ve>i  defines  a  semi-major  axis  for  the  assist  hyperbola; 
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a 


2 


(B-18) 


Right  Ascension  and  Declination  of  For  the  purpose 

of  defining  the  geometry  of  the  assist  hyperbola,  it  is 
necessary  to  determine  the  angles  a  and  6  shown  in  Fig.  46. 
Later,'  when  the  assist  hyperbola  has  been  completely  defined, 
it  will  be  seen  that  these  angles  are,  respectively,  the 
right  ascension  (measured  eastward  on  P2  in  the  ecliptic 
plane  from  the  direction  of  the  Y2  axis)  and  the  declination 
(measured  north  (+)  or  south  (-)  of  the  ecliptic  plane  along 
the  P2  meridian  determined  by  a) . 

Both  a  and  6  are  found  with  the  vector  V„.  which 

X0Y0 
2  2 

is  defined  as  the  component  of  Vwi  in  the  X2~Y2  plane.  As 
shown  in  the  figure,  the  magnitude  of  this  vector  is 


Va 


■x2Y2 


1/2 


(B-19) 


Thus,  for  the  right  ascension, 


cos  a  « 


~V”iY2 


V-i 


x2y2 


(B-20) 


and  for  the  declination, 


V“i 

cos  6  =>  _ 


V 


(B-21) 


155 


GA/AE/67-4 


Selection  of  3.  The  position  of  the  assist  orbit 
plane  is  defined  by  the  center  of  mass  of  P2.  Its 
orientation,  like  any  other  plane,  must  be  defined  by 
three  angles.  Two  of  these  angles,  o  and  5,  have  already 
been  calculated,  but  a  third  angle,  6  (defined  below),  must 
be  specified. 

After  3  is  selected,  a  coordinate  transformation  (see 
Fig.  68)  is  made  from  the  X2Y2Z2  system  to  the  x2y2z2  assist 
orbit  system  as  follows: 

1.  The  x2y2z2  axes,  initially  aligned  with  the  X2Y2Z2 
axes,  are  rotated  about  the  Z2  axis  through  the  angle  a. 

This  rotation  defines  the  primed  system  x2,y2,z2’  shown  in 
the  figure. 

2.  The  primed  system  is  rotated  through  the  angle  6 
about  the  x2'  axis.  This  rotation  defines  the  double  primed 
system  x2' 'y2r,z2' ' ,  and  aligns  the  y2''  axis  with  the 
incoming  HEV  vector. 

3.  The  double -primed  system  is  rotated  about  the  y2" 
axis  through  the  angle  3,  into  the  x2y2z2  system  so  that  the 
x2-y2  plane  lies  in  the  plane  of  the  assist  orbit  and  z2 
forms  a  right-handed  system. 

The  transfdrmation  is  expressed  as  follows* 
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(c3  ca  -s3  s<5  so)  (c3  sa  +  S3  s6  ca)  (-S3  cS) 

{-c5  sa)  (c6  ca)  (  s6) 

(s3  ca  +  c3  s6  sa) (s3  sa  -  c3  sfi  ca)  (c3  c6) 


"1 

f 

V 

A2 

>< 

*2 

z„ 

2 

^  / 

and  since  the  transformation  is  orthogonal, 


where  D  is  defined  as 

< 


D  = 


(c3  ca  -  s3  s6  sa) 
(c3  sa  +  s3  s6  Ca) 
(-S3  c6) 


(-C6  sa) 
(  c6  ca) 
(s6)- 


(s3  ca  +  c3  s6  sa) 
(s3  sa  -  c3  s6  ca) > 
(c3  c 6) 


(B-24) 


Selection  of  d.  The  position  and  orientation  of  the 
assist  orbit  plane  has  been  determined  above,  but  a  DOCA 
must  be  selected  in  order  to  define  the  assist  orbit  itself. 
This  procedure  is  analogous  to  the  two-dimensional  case  in 
which  the  assist  orbit  plane  was  automatically  defined  by  the 
ecliptic,  and  a  DOCA  was  selected  to  define  the  assist  orbit. 

The  value  of  d  determines  the  eccentricity  of  the  assist 


orbit 


e2  “  _i  (Rp  +  a2  +  dj  (B-25) 

a2  2 

The  angle  between  the  conjugate  axis  and  the  inbound  hyper¬ 
bolic  asymptote,  as  shown  in  Fig.  69,  is  given  by 
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sin  v  =  _i  (B-26) 

e2 

This  angle  is  used  to  find  the  components  of  V«0  in  the 
x2y2z2  system, 


f  V 
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V„0  sin  2\> 
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{  V»o 

>  «  < 

-V»0  cos  2v  ^ 
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> 

0 

<  z2; 

where,  of  course, 


(B-27) 


V. 


V»i 


(B-28) 


The  components  of  V«0  in  the  X2Y2Z2  system  are  found 
with  the  following  transformation: 


(B-29) 


Finally,  the  actual  position  of  02,  the  origin  of  the  X2Y2Z2 
system,  is  determined  by  the  vector  rm  from  02  to  the  center 
of  mass  of  P2.  For  the  components  of  this  vector  in  the 
X2Y2Z 2 system. 
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(B— 30) 
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where,  rm  “  a2  +  ^2  +  d 


(B-31) 


Post-Assist  Orbit 

Post-Assist  Orbit  Conditions  at  P2.  The  heliocentric 
velocity  of  the  spacecraft  as  it  leaves  the  SOI  of  P2  is 
given  by 


v3  =  V-0  +  Vp2  (3-2) 

as  shown  in  Fig.  47.  This  equation  may  be  written  in 
matrix  form: 
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(B-32) 


(B-33) 


The  flight  path  angle,  $3,  is  given  by 


cos  $3  = 


*X2 


V  3 


(B-34) 


and  the  inclination  of  the  post-assist  orbit  is 


i3  =  arc  tan 


(B-35) 


The  true  anomaly,  or  the  heliocentric  angle  from  P2  to 
perihelion  is 
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Lrp  V\  /v  ]  cos  sin 

tan  f  3  -  LlL— . . . _L—  <3~36> 

(rr-2  V|  Ae]  sin2  $  3  -  1 

Post-Assist  Orbit  Characteristics.  Finally,  the  post- 

i  n  ■«  .  ■  i».»  i  n  mm  i  M  mmm _  “* 

assist  orbit  characteristics  are  obtained  from  the  follow¬ 
ing  relations: 


„  rP2 

2  -  (rB2  V3  /%} 

(A-27) 

2 

V2  /VQ  “  l]  sin2$3  +  C0S24>3 

(A"  2  8) 

12  “  ^  ^©f  e3/  ^3^ 

(A-29) 

r  =  a3  (1  +  e3) 
a3 

(A-30) 

qm  =  a3  (1  -  e3) 

(A-31) 

Summary 

A  complete  summary  of  the  velocity  vectors  at  P2  is 
shown  in  Fig.  70.  This  diagram  is  similar  to  Fig.  16(a), 
and  may  be  used  to  separate  the  events  into  two  reference 
frames:  one  with  respect  to  P2,  and  one  with  respect  to 
the  Sun. 
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Appendix  C 


Pre-Assist  Orbit  Inclination  for  Category  II  Trajectories 


The  purpose  of  this  appendix  is  to  derive  an  expression 
for  the  pre-assist  orbit  inclination  i  required  for  given 
values  of  burnout  velocity  and  departure  flight  path 
angle  <|>o  for  Category  II  trajectories.  The  problem,  illus¬ 
trated  in  Fig.  71  ,  is  to  find  an  orbit  which  intersects  the 
orbits  of  Pi  and  P2  at  a  distance  of  rp^  and  rp  from  the 
Sun,  respectively. 

An  expression  for  i  is  obtained  from  Eq  (2-8) : 

V*  -  2Vr>  sin  <j>„  cos  i  Vfl  +  v|  -  vip  *  0  (2-8) 

o  Pi  o  0  pi  *1 


which  gives 


/V§  +  v|t  -  vjp  X 

i  a  arc  cos  _ 

\  2vp,  V0  sin 
1 


(C-l) 


The  quantities  Vp^,  <fr0,and  Veop^  (found  with  given  Vfa  in 
Eq  (A-l) )  are  all  known  in  this  expression,  but  V0  must  be 
calculated. 

From  the  energy  integral  for  elliptic  orbits, 


V 
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(C-2) 


where  ai  is  the  serai-major  axis  of  the  pre-assist  orbit. 

The  'r  ilue  of  is  found  by  first  writing  the  basic  orbital 
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equation  for  the  times  of  departure, -from  Pj  and  arrival 
at  P2: 
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ai  (1  -  e\) 

1  +  ex  cos  fQ 


(C-3) 


rP 


2 


=  al  <1  ~  el> 

1  +  ei  cos  fA 


(C-4) 


where  ej  is  the  eccentricity  of  the  pre-assist  orbit. 
Division  of  Eq  (C-3)  by  Eq  (C-4)  gives 


1  +  e  j  cos  f  j 


1  +  ex  cos  f0 


(C-5) 


For  Category  II  trajectories,  by  definition, 


*l 


(C-6) 


where  f0  is  assumed  to  be  a  negative  angle  for  pre-assist 
orbits  inside  the  orbit  of  Pj,  Substitution  of  Eq  (C-6) 
into  (C-5)  with  the  definition  R  =  rp^  Ap2  gives 


1  -  e,  cos  f 
R  *»  -  .  - _ . _ 1 

1  +  el  cos  f0 

which  can  be  rearranged  to  give: 


(C-7) 


ej  cos  f0 


(1  -  R) 
(l  +  R) 


(C-8) 


Now,  from  orbital  mechanics  (Ref  16:83)  it  can  be  shown 

that 
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sin  $0 


1  +  ej  cos  £0 
(1  +  ef  +  2ej  cos  £0) */2 


Definition  of  the  quantity 


(C-9) 


A  «  Lr-J 5  (C-10) 

1  +  R 

and  substitution  of  Eq  (08)  into  (C-9)  gives 


sin 


1  +  A 

U  +  e?  -r  2A]  V2 


(C-ll) 


An  expression  for  e.  is  found  by  squaring  and  rearranging 
Eq  (C-ll) s 


e.  »  _ (A2  +  (1  -  sin2$0)  (1  +  2A)]1/2  (C-12) 

sin  $0 

which  can  be  determined  from  the  known  values  of  A  and  $0. 

The  value  of  ex  from  Eq  (C-12)  can  now  be  used  in  the 
following  equation  (Ref  16s83)  to  find  a^ 


a\  (1  -  e?) 

sin2*°  =  r — HZ - TT 

r?i  (2ai  rp,J 

This  expression  is  rearranged  as  follows i 

(1  -  ef)  af  -  2rPi  sin2$0  al  +  r|^  sin2$0  > 

which  is  a  quadratic  equation  for  a1#  Solution 
gives 


(C-13) 


■  0  (C-14) 

of  Eq  (C-14) 
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ai  = 


rn  sin 


(1  -  e2) 


sin<{>0  *  (sin2$0  +  e2  -  1)  l/2 


(C-15) 


where  the  -  (+)  sign  is  used  when  the  orbit  of  P2  lies 
inside  (outside) -the  orbit  of  P1# 

In  summary,  the  sequence  of  calculations  for  the 
required  inclination  of  a  pre-a  ;sist  orbit  is  as  follows: 
1)  Given  values  of  rp  ,  r?2,  Vp^,  $Q,  and  Vfa 


2)  _  K  -  2UPi 


V, 


“Pi 


(A-l) 


3) 


4) 


1 -  [a2  +  (1  -  sin2*0)  (1  +  2A)]  1/2  (C-12) 


ai  =  — 


sm  $0 
rPl  sin  <}> 


(1  -  e2) 


[sin  <{»0  -  (sin2 4>0  +  e \  -  1)  V2] 


5) 


V2  =  u_ 
o 


2.1 

rPl  al 


6)  /  V 2  +  V2  -  V2r,  \ 

i  =  arc  cos 


o  'Pj  ywPi 


2VPi  V0  sin  <l>( 


(C-15) 

(C-2) 

(C-l) 


/ 
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Appendix  D 
Launch  Opportunities 

The  initial  assumption  of  circular  and  coplanar 
planetary  orbits  means  that  any  calculation  of  launch 
opportunities  is  only  an  approximation.  However,  this 
approximation  is  consistent  with  other  work  in  the 
study,  so  the  results  may  be  followed  as  a  guide  for  more 
precise  studies. 

Fig.  72  illustrates  the  heliocentric  longitudes  110 
and  120  of  planets  ?1  and  P2,  respectively,  with  respect 
to  the  vernal  equinox  direction  T  at  some  epoch  tQ.  For 
the  longitudes  at  any  time  t,  after  tQ: 

lx  =  110  +  (t  -  t0)  (D-l) 

12  =  l2o  +  «2  -  t0)  (D-2) 

where  e 1  and  oj2  are  the  mean  daily  motions  of  P:  and  P2, 
respectively. 

To  accomplish  some  given  transfer  (pre-assist  orbit) 
in  time  Ti2  through  an  angle  e^,  from  planet  P2  to  P2, 

“ 6 ^  =  u2  T^2  “  ®12  (D-3) 

which  gives  the  lead  angle  eL  measured  from  Earth  longitude 
at  launch,  as  shown  in  Fig.  73. 

To  find  a  launch  time  for  the  given  trajectory,  let 
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Fig.  72 


Position  of  and  P2  at  Epoch 
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~8l  =  lj  -  12  (D-4) 

and  from  Eqs  (D— 1)  and  <D-2) ,  for  t  =  fcL ' 

~®ja  =  ljo  *  ui  (^l  “  t0)  -  12q  “  w2  “  fco)  (D-5) 

Finally,  after  rearranging  Eq  (D-5) , 


tL  =  - i -  [dio  “  I20)  +  (u2  “  ul>  t0  -f  0L]  (d-j6) 

(u2  “  wl) 


which  gives  the  first  possible  launch  time  near  t0.  In  this 
thesis  the  epoch  tQ  was  taken  as  1.5  Jan  1960  which  is  Julian 
day  number  2,436,935. 

Successive  launch  dates  for  the  given  transfer  orbit 
may  be  generated  by  adding  multiples  of  the  synodic  period 
to  tL.  Thus,  for  n  launch  dates  in  addition  to  the  date  tL 
from  Eq  (D-6) , 


t 


L 


l 


tL  +  S 


^2  =  +  2S 


tL  +  ns 


(D-7) 


where  S  is  the  synodic  period  of  Pj  and  P2. 

The  angle  ep,  measured  from  P2  to  P2  at  the  time  of 
intercept,  is  given  by 


9p  =  0l2  -  Ul  t12 


(D-8) 
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Launch  dates  are  generated  with  Program  8  for  ecliptic 
and  out-of-ecliptic  orbits.  The  method  for  calculating 
6 12  and  T12  is  similar  to  that  used  in  Program  6  for  the 
ecliptic  and  that  used  in  Program  7  for  the  out-of-ecliptic 
trajectories.  Eqs  (D-l)  through  (D-8)  apply  to  both  types 
of  trajectories. 
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Appendix  E 
Computer  Programs 
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Program  1 

»j08  0.4. 10000  67-012<LT  M»ERS<AF|T-Sr 

3I8JOS  MAP 

JJ6FTC  MAIN  H94/2<XR7 

c  DIRECT-TRANSFER  TRAJECTORIES  FOR  SOUR  PROSES 

TIMC(A«0<E»T)«l2.«ATANISO«TI(l.-EI/ll.*Ell»TAN(TAR/l.M-E«SORm.- 
1E<<2 )*S!NtT»R I /< 1<<E<COS<T<RI1 1  <S0RT<<A<431 /6l/«6400< 
REA013<11GCS<GCP1<RP1<VP1<RA0P1 
1  FORMAT <  3E10<0 1 
YB*S4000< 

Pl<3<14139 

R.PI/100. 

M<72 

N»9 

L<19 

00  2  I»1<M 

T<VBM2-2<*GCP1/RADP1 

IF<T1102<103<103 

102  WRITEI6<104)VS 

104  FORMAT<1KA<10X<3HVB<<F10<1<21HBEL0W  ESCAPE  VELOCITY) 

00  TO  31 

103  VINFP1<S0RTITI 
WRITE(6tl051 

103  FORMAT  1 1H  <2X<66H<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<M<M<MM<<<<MM< 
IWHIHHimiMI) 

PH0«9O. 

00  3  J-l.N 
PH0R-PH0»R 

WRITEI6.100IVB>VIHFP1.PM0 

100  FORMAT  1 1HA»10X«3HVB*  *F10»l*3X*7HVINFPl<»F9*l» 3X«3KPHB*F3»l//4Xt2HX 
1 1 »8X»2HF0*SX*2HA1»8X*2HE1 *7X»5HTIH0P  »6X  t2HRPl 
XI*0. 

00  4  X«1.L 
XIR«XI*R 

IFI(VP1»SIM(PHORI»COSIXIR)|4*2*V!NFP14»2-VP14*2)S*3.6 
3  WR1TE(6»7 JVBiPHOtXI 

7  FORMAT  1 1HA#  4X«17HNO  ESCAPE  FOR  VB<<F10<1<7HAT  PH0<<F7<2<6HAM0  t>< 
1F6.2I 
GO  TO  30 

6  VO-VP1<S1N(PHORI<COSIX1R1-50RT  I<VP1<SINIPHOR1<COSI  XIRM<<24VINFP1< 
1«2-VP1«2) 

RP1VO<RP1*V04*2/OCS 

TANFO*RP1VO<CQS(PKOR)<SIN(PHOR)/(1.-RP1VO*SIN<PHORI”2) 
IF<TANF0U<9<9 
•  FO-ATANITANFOI/R 
00  TO  10 

9  F0<<ATAN(TANF0)-PI1/R 


10  A1<RP1/I2.-RP1V0) 

A1AU-A1/<92901000.»3200<) 

El-SO«Tt<<RPlVO-l.  I  <SIN1PH0R)  1M24COS<PHOR1<<21 
RPER1-A1*! 1<-E1) 

RPAU<RPERl/< 92901000<<3280<  I 
TIHOP— TIME<A1<0CS<E1<F0I 
WRITE<6<101IXI<FO<A1AU<E1<T1MOP<RPAU 
101  F0RMAT<6F10<41 
Xt*XHS. 

4  CONTINUE 

30  PH0«PH0-10< 

3  CONTINUE 

31  VB<VB<1000< 

2  CONTINUE 

STOP 

END 


3  DATA 

4<6868E21  1.4077E16  4<9032E11  9<7731E04  2.0902E07 
f  EOF 


175 


GA/AE/67-4 


Program  1  Output 


V8=  42000.0  V INFP1  =  20421.7  PH=  90.0 

XI  FO  A1  El  TIMOP  RP 

0.  -179.9999  0.7277  0.3742  113.2736  0.4554 

5.0000  -179.9999  0.7409  0.3498  116.3526  0.4817 

10.0000  -179.9998  0.8030  0.2454  131.2852  0.6059 


NO  ESCAPE  FOR  VB=  42000.0AT  PHO=  90.00AND  1=  15.00 


VB=  42000.0  V  IN  FP 1=  20421.7  PH=  88.0 

XI  FO  A1  El  TIMOP  RP 

0.  -176.6150  0.7297  0.3719  109.3978  0.4583 

5.0000  -176.2300  0.7432  0.3471  112.1644  0.4852 

10.0000  -173.6313  0.8082  0.2398  125.1468  0.6144 


NO  ESCAPE  FUR  V8=  42000. OAT  PHO  =  88.00AND  1=  15.00 


VB=  42000.0  V INFP1=  20421.7  PH=  86.0 

XI  FO  A1  El  TIMOP  RP 

0.  -172.9760  0.7359  0.3648  106.2147  0.4674 

5.0000  -172.1182  0.7506  0.3388  108.7804  0.4963 

10.0000  -165.5901  0.8264  0.2209  120.4026  0.6439 


NO  ESCAPE  FUR  VB=  42000. OAT  PH0=  86.00AND  1=  15.00 


VB=  42000.0  VINFPi=  20421.7  PH=  84.0 

XI  FO  A1  El  TIMOP  RP 

0.  -168.7244  0.7474  0.3520  103.5780  0.4843 

5.0000  -167.1456  0.7647  0.3235  106.0049  0.5173 

10.0000  -149.7241  0.8747  0.1767  114.6357  0.7202 


NO  ESCAPE  FOR  VB=  42000. OAT  PH0=  84.00AND  1=  15.00 


VB=  42000.0  V  INFP1=  20421.7  PH=  82.0 

XI  FO  A1  El  TIMOP  RP 

0.  -163.1628  0.7671  0.3313  101.3462  0.5129 

5.0000  -160.1405  0.7900  0.2978  103.5839  0.5547 


NO  ESCAPE  FOR  VB=  42000. OAT  PH0=  82.00AND  1=  10.00 
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Program  2 

0.3.4OCO  4T-012.  NYERS.APIT-SE 

SIF.MB  HAP 

1IBFTC  MAIN  M94/2.XRT 

C  VIHUS  A 55 1 i TCO  TRAJECTORIES  IN  THE  ECLIPTIC 

rrME(A.G,£.T|.l2.*ATANISORTl(l.-E|/(l,AE»l*TAN(TM/2.>l-e»SORT(l.- 

lE»*2l*SlN(T*RI/<l.+E»COS<T4Rm»SORT((A4M|/GI/B4AOO. 

C  CALCULATION  OP  DEPARTURE  ORBITS 

REA  5(5.1 )GCS,GCP1.GCPZ,RP1»RP2.VP1.VP2.RADP1.MA0M2 

1  FORHATISElO.O/AElO.Ot 
P!*>.1A1S9 
R-PI/180. 

C  DIVERSION  FOR  ONE  GREATER  THAN  REQUIREO 

DIMENSION  VDI  162  I  .PHI  142>10).FO(1A2.10).VOI  142.10)  .A  ( 142. 10)  >A0(  It 
I2i  KIM rr 1162. 1 0  J >V INFP1 1162) .01221 ,£2122) .IMVOC4I22I .PI IOCO (22 1 (VS 
2(22  .F3I22) .A»(22).ES(221.0H(22I*J1( 142). TOPI  221 
C  Villi )  HOST  3E  GREATER  THAN  ESCAPE  VEL  •  S4T00  FT/SEC  FOR  EARTH 
VB(1.)  ■59000. 

N*3’i 

N-7 

L-21 

oo  i  ■  l  .m 

VI NI  P 1 1  1 1*504 f  IV8! II»*2-2.*OCPl/RAOPl> 

PHI  I .1 ) «R0. 

oo  ;i  j-i.n 

PHR.PHI I.J)*R 

IF(VIRF«H  I  )•  »2-lvPl*COS(PHRn**212A.2*.2J 
2 A  Jill  1-J-l 
GO  10  29 

25  TANi:TP*VPl»COS  IPHRI/ISORT  IVINFP1(II**2-IVP1*C0S  tPHRII**2)l 
Jill  >*J 

ETI1  .J).|PI-AfAN(TANETPI1/R 
At  I . Jl*P|/2«-ETI 1 • J I »R+PHR 
ADI  I • Jl *A( I . J I /R 
IFI1HI I  .JI-90. 119. 20,19 

20  VOII »J|«VP1“VINFP1( II 
GO  10  21 

19  V0( 1 .JI-VINFP1I I 1*SIN  (AII.J1I/C0S  IPMRI 

21  RPVC'*RP  1*  (  VO  I  1  » J 1 **2 I  /GCS 

TANFO-RPVO»COS  I PHR 1*51 N  I PHR 1 / 1 1 .-RPVO»S I N  (PHRl»»2) 

FOM.JI«(ATAN  ITANFOt-PIt/R 
PHI  I  .J+l 1 - PH I  ! . J 1-5. 

3  CONTINUE 

29  VBII+ll-VDIIl-lOOO. 

2  CONTINUE 
WRITEI4.A) 

A  FORMAT I  1H1.50K.21 HEARTH  DEPARTURE  ORRITI 


DO  6  I.l.H 
J2-JKII 

NR I1EI4.5IVBI I  I .VI NFP 1 1 1 1 . (PHI I.JI.FOII.JI « V0( I . J I .ADC l.Jl.ETll.Jl 

1.J.1.J2I 

5  FORMAT  1 1HA.25X.3HVB*.F10.0.MX.7MV1NFP1«.F10.0/ITX.2HPM.I»A»2HFO.I 
14X.2HV0.1BX.IHA, l9X.2HtT//(3H0.2ll 

6  CONTINUE 

C  CALCULATION  OF  PRE-ASS  1ST  ELLIPSE 
00  13  1*1, M 
J2-J1 1 1 1 
D 1 1 1 *0, 

00  13  J*l.  J2 

A1*RP1/ (2.-RP1* I VOI I , J) **2 I /GCS) 

A1AU-A1 2 (92901000. *3250.1 

E1*S0RT  ((IRPl»(V0(l,Jt**2l/GCS-l.l*SIN  (PHI  I ,JI*R| I =»2tC0S  IPHII. 
1J)»R»R»2I 
RPER1*A1»( l.-El) 

IF (RPER1-RP2 113.13,14 
16  NRlTEIt.lTI 

IT  FORMAT UHA. 3 1HVB  INSUFFICIENT  TO  INTERCEPT  P2t 
00  TO  IS 

13  Vl’SQRT  (GCS* 12, 2RP2-1.FA1 1 1 

SPHIl.SQRT  ( (A1**21*I1.-E1**2I/(RP2RI2.*A1-RP2»>) 

CPHU.SORT  (l.-SPHIl*»2l 

TANF1.(RP2»(V1*«2|20CS1*CPHI1*SPHI1/(1,-IRP2*(VI»«2I/0CTI»SPHI1»»2 

11 

IF ITANF1 IA0.A0.A1 
AO  F1-ATANITANF1I2R 
GO  TO  A2 

A1  F1*(ATAN  (TANF1I-P1 1/R 
A2  THT12*F1-F0( I .J) 

T IH12«TINE( Al.GCS.El.Fl 1 -TIMEIAl .CCS, ei.FOl I. Jl I 
C  CALCULATION  OF  HYPERBOLIC  PASS  TRAJECTORY 

VI NF*SORT  ( VI *»2*VP2**2-2.»Vl»VP2»SPMI 1 » 

A2*GCP2/(VlNF*»2l 
IF (VP2-V1»SPH!1) 8.7.9 
T  WR1TEI4.1A) 

1A  FORMAT! 1HA.26HVF2*V1*SPH11  FOR  THIS  CASH 
GO  TO  11 

8  TANNPH— V1*CPHI1/SQRT  «VINF»»2“1VI«CPHI  1I»»2I 
GO  TO  10 

9  TANNPH*V1*CPH1 1/SORT  (VINF«2-IV1»CPHI1)»»2I 
10  DO  11  K«1,L 

E2(XI*(RAOP2+D(XI»408O.27*A21/A2 
XNU-ATAN  (l./SORT  ( E2 ( X l **2-1 . It 
XNUOEOIKI-FNU/R 
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Program  2  (Contd.) 


IFtVP2-Vl*SPHIllAA,AA,A3 
AA  PSI*ATANtTANNPHI-2.»XNtt+PI 
00  TO  SO 

AS  PSUATAN  (TANNPHI-2.PXNU 
SO  PS10EO<E)*PS!/X 

VJtXUMNT  IVPJ»»2*VlNr*»2-z,«vP2*ylHP*COS  IPJMI 
CPHI  J.VINFMIN  (PS!)/V3(K) 

SPHIJ-SONT  (l.-CPHIJ»»2> 

APV  J-RP2*  t  V  J  ( X  U»2 )  /ACS 

TANF J-RPV J»CPH I J»SPHI S/ ( 1 . -NPV J» JPN J J»»2 ) 

IPtTANPJU5.A4.AA 
AS  FJ<KU(ATAN<TANFS>1/R 
00  TO  A 7 

AA  FJtKUtATAN  ( TANPSI-Pl  I /* 

A7  ASIKU(RPi/(2«“)ffVJ ))/' 7290 J 000. *5210. 1 

EJtKLSONT  <<tRPVJ-l.t*SPHIJU»JACPWIJ942) 

AA-RP2/(2.-RPVJ> 

T0PIKUTIN12-TINE(AA.GCS.EJ(K».FJ<K)I 
OMtKUAJtKUt  l.-EJIK) ) 

0IKA1U0IKU100. 

U  CONTINUE 

WRITEI4.12IVBIII.PHII  .J1  .  I THT  12.F1 . AlAU.El . Vl.VINF.Df  X)  ,FJtX)  .AJU 
1).EJ(XI.VJ(X).XNUOSO<XI  .E2<X>.PJIDE0tK|*TIHll»T0PIK).Q«IXI.K»l.U 

12  FORHATt 1HA. 2 XX .SHVB* >F10 .0.A7X.SHPH* .F10.2//lX«AHTH12.tX«2HFl .JX.2 
1HA1 .SX.2HC1 .7X.2HV1 »6X» AHVINP .AX.AHOOCA.AX.2HPJ.AX.2HAJ.9X.2HC J.4X 
2 .2HVJ  .SX.2HNU.SX.2HE2.SX.  JHPSI  .AX.5HTIN12.JX.JHT0P.AX.2MM//IP7.2. 
JFS.2.2F7.A.2P9.1.PS.0.P».2,2P7.A.P».1.PA.2,P4.J.P«.2.P7.1.PT.1.PA. 
AS)  I 

13  CONTINUE 
STOP 
END 

SOATA 

A.6S6SE21  1.A077E16  1.1A44E14  A.9052E11  J.SAAOEU 
9.7751E0A  1.1A9JE0S  2.0902E07  2.O01JE07 
JEOP 
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Program  3 

SJ08  0.5.8000  67-01 2* MYERS .AFIT-SE 

SIBJGB  MAR 

S IBFTC  MAIN  H9»/2.XR7.NOD£CX 
C  VENUS-HERCURY  PERIHELION  ASSISTED  TRAJECTORIES 
C  CALCULATION  OF  DEPARTURE  ORBITS 

TlHE(A>G.£.T)-(2.*ATANtSQRTIU.-£)/!l..E))*TAN(T*R/2.)!-E*SORm.- 
l£»*2)*SIN(T»fi 1/(1. ♦£*«>$( T*R|))*SOAT((A**3)/G)/86*00. 

READ! 5.1  IOCS. 0CP1.0CP2.PP1.RP2.VP1.VP2.RADP1.RADP2 
1  FORMAT<5£10.0/*£10.0) 

101  READ  1 5 . 100IGCP3.RA0P3.RPS. TA3 
100  FORHATIAEIO.O) 

C  DIMENSION  FOR  ONE  GREATER  THAN  REQUIRED 

DIMENSION  VB 1 162 1 .PH 1 162. 101. FO 1 162. 101. VO  1 162. 101. AT  162.101. ADI 16 
12.10) *ET( 162.10) .VlKFPl 1 1621 .0(221 *E2 (221 »XNUDEG( 22) .PS1DEGI 221 .V2 
2(22 1 .F3I22) .A3I22) »E3(22) .QM(22) .J1S 162) .AAI221.T122) 

DIMENSION  D3(10).E*(10).XN*DEG(10).PSI*DG(10).V*(10).F7(10).AS(10) 
1 *A5AU( 1 3) *E5( 10) »OKM( 10).THT3P(10).T3P(10).THTOT(10)  .TIMOR! 10) 

C  V 3 ( 1 )  MUST  BE  GREATER  THAN  ESCAPE  VEL  •  36700  FT/SEC  FOR  EARTH 
VB( 1 ) *38000. 

PI-3. 1*159 
R-P.I/180, 

VP3-S0RT ( 2 .*GCS*RA3/ ( RP3* ( RA3+RP3 ) 1 1 

M-2S 

N-A 

L-6 

Ll-6 

DO  2  I-l.M 

VINFP1 ( I l-SORT  (V8( I I»*2-2.»GCP1/RA0P1 ) 

PHIt. 11-90. 

DO  3  J-l.N  .. 

PHR-PHI  I.J)*R 

IF<  VINFPK I)»*2-(VP1*C0S(PHR)!»*2)2*.2*.23 
2*  J1III-J-1 
GO  TO  29 

23  TANETP-VP1*C0S  (PHR)/(SORT  (VINFPK 1 ) **2- <VP1*C0S  (PHR))**2M 
J1ID-J 

ET( I .J)-IPI-ATAN(TANETP) l/R 
A(I.J)-PI/2.-ET(l • J)*R.PHR 
ADI  I  .  J 1  —  A ( I  .  J)  /R 
IF(PH(I.J)-90. 119.20. 19 

20  V0( I . J) -VP1-VINFP1 ( 1 1 
GO  TO  21 

19  V0(  I  .  J)  -  VINFPK  1 1 -SIN  (AII.JII/COS  (PHRI 

21  RPV0-RP1*(V0( I .JI**2I/GCS 

TANFO-RPVO*COS  (PHR)*SIN  (PHR) /( l.-RPVO*SIN  (PHfi|»«2) 

FOd.Jl-IATAN  (TANFOI-PIl/R 


3 

29  VB(I+1)-VB(I)*1000. 

2  CONTINUE 
WRITE-6.*) 

*  FORMAT! 1H1. SOX. 21HEARTH  DEPARTURE  ORBIT! 

DO  6  I-l.M 
J2-J1III 

WRITE(6.5)VB(l).VINFPl(II.(PM(t.J|.F0(I.J).V0(I«J)>AD(I»J|»ETtI.J) 

I.J-1.J2) 

5  FORMAT(lHA.29X.3HVe-.(-10.0.33X.7HVl(IFPl-,F10.0/lTX.2HPH.18X.2HF0.1 
18X.2HV0. 18X.1HA, 19X .2HET// (5F20.2  1 1 

6  CONTINUE 

C  CALCULATION  OF  PRE-ASSIST  ELLIPSE 
DO  13  I-l.M 
J2-J1II) 

0(11-0. 

DO  13  J-1.J2 

Al-RPl/(2.-RPl*(V0( I.J)«»2)/GCS! 

A1AU-A 1 / ( 92901000. >5280. ) 

El-SORT  ( ( (RP1*(V0( 1 . J)**2 l/GCS-l. )*SIN  (PHI  I • J)*R) )**2.C0S  (PHI  I > 
1JI»R)**2) 

RPER1-A1*( 1.-E1) 

Ir'.RPERl-RP2)15.15.16 
16  WRITEI6.1T) 

IT  FORMAT! 1HA.31HVB  INSUFFICIENT  TO  INTERCEPT  P2) 

GO  TO  13 

15  Vl-SQRT  IGCSM2./RP2-1./AI ) ) 

SPHI l-SORT  ((A1»*2)»(1.-E1»»2)/IRP2*(2.*A1-RP2)() 

C PHI  1-SORT  (1, -SPHI 1**2) 

TANF1-1RP2*! V1**2'/GCS1 *CPHI 1*SPHI 1/1 1«-IRP2*IV1**2)/0CS)*SPHI 1**2 
1) 

IFITANF1!*0.*0.*1 
*0  F1-ATAMITANF11/R 
GO  TO  *2 

*1  Fl-IATAN  (TANFl)-PI l/R 
*2  THT12-F1-F0II.J) 

TIM12-T I ME (Al.GCS.EltFl l-TIMEI Al.GCS.El.FOI I • J>) 

C  CALCULATION  of  HYPERBOLIC  PASS  TRAJECTORY 

V1NF-S0ST  !V1**2*VP2**2-2.»V1*VP2*SPHI1I 
A2-GCP2/(VINF**2) 

IFI VP2-VI*SPH1 1 )8*7 .9 
T  WRITE-6.1*) 

l*  FORMAT! iHA.26HVP2-Vl*SPMIl  FOR  THIS  CASE) 

GO  TO  13 

8  TANNPH— V1*CPHI  1/SORT  (VINF*»2-(V1*CPHI  1  )*»2I 
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so  to  10  Program  3  (Contd.) 

9  TAK*PH-V1«CPHI  1/SORT  (VINF**2-(  V1*CPHI  H«»2) 

10  00  11  X-l.L 

E2{Kl"<RADP2*OUJ«60t0.27+A2»/A2 
XNU-ATAH  11. /SORT  I E2(X>**2-1. M 
XNUDEG(Xl-XNU/R 
!F(VP2-V1*SPHI1>44.44.43 
44  PSl«ATA*ICTAXHPm-2.*XNU4PI 
GO  TO  30 

43  P31-ATAH  (TAMKPHJ-2.4XMJ 
30  PSlD£CtX»«PSI/R 

V2  (X)  "SORT  (VP2«»2*V!AM»2-2.*VP2*VINMC0$  (?$!)! 

CPHI 3"VlNF*5IN  (PSt)/V2(X) 

SPHI3-SORT  (l.-CPHI3»*2» 

RPV3«RP2*(V2(X)4*2|/OC$ 

TANFS-RPVS'CPHlS^SPHISm.-RPVSiSPHI?*^) 

IF(TAKF3143»46.46 
43  F3IX) •( ATANITANF31 ) /R 
GO  TO  47 

46  F3UI-CATAN  ITARF3J-PD/R 

47  A3(X)*CRP2/C 2.-RPV3 It 

AA(F.)  "AS  (X  )/(  9290 1000 .*3280.  ) 

E3(X)-SOR7  mRPV3-l.J»SPHl3)«*2+CPH!3»*21 
QNtX)»AA<X)»(l.-€3(XJ> 

T2P— TIN£lA3lXj.GCS.ES(X|.F3(Xn 
T(X»**T2P4TIM12 

dcx+ii-oixi+ioo. 

11  CONTINUE 

VR ITE (6. 12 ) VB( I I.PHI I  •  J ) • ( THT 1 2 »F1 » A1AU.E1 »V1 »T( X ) «0<X ) .F3IX) »AA(X 
11.E3IX).V2<X),XNUDEG(X>*TIM12.PSIDEG(X)  .QN(X)  .X*1»L) 

12  FORMAT ( 1HA  »23X.3HVft*»Fi0»0»47X  »3HPH* .F10.2//1X.4HTH12 .3X.2HF1 »3X»2 
lhAl »3X .2HE1.7X.2HVI »6X»4HTT0T »5X .4HD0CA.4X.2HF3.6X. 2KA3  *5X. 2HE3.6X 
2.2HV3*8X»lHN»3X»3HT12'5X'3HP$!«6X'2HOH//IF7.2»Fa*2*2F7»4*2FloU»F7 
>.1.F8.2.2F7.4.F11.1.F8.2,F7.2.F8.2.F8.SM 

102  DO  13  K*l »L 

I FIQMCX)#92901000. •3280. -RP3) 104.104.121 
121  WRITE<6.)C3> 

103  FORMAT ( 1HA.31HV0  INSUFFICIENT  TO  INTERCEPT  P3> 

GO  TO  13 

104  03(11*0. 

Y3»SORT(GCS»C2./RP>-l./A3UH> 

5PMI5-SORT((A3IX)«*2l»Cl,-t3(X>*«2)/(RP3*(?.*ASIA)-RP3m 
CPHI  3 "SORT  1 1 .-SPHI3**21 
C  HYPEROOUC  TRAJECTORY  AT  P3 

VlNF3"SORTIV3«*2+VP3»»2-2.*V3»VP3fcSPHl3| 

TANF3*<RP3*lV3#*2)/GCS)«CPHl3*SPHI3/a.-(RP3*(V3«°2l/GCS»#SPHl3»*2 


1) 

IFITANF31103. 103.106 
105  F3-ATAN(TANF3J/R 
GO  TO  107 

10*  F5-IATANCTANF3)-.PI)/R 

107  THT23-F3-F3IX) 

T23*TIMC(A3(X)  .GCS.E3tX),F5>-TlHt(Aj{X).GCs»E3<X»,F3(XU 

A4«GCP3/(V1NF3**2I 

IFCVP3-V3»SPH13I108.109,110 

109  WRnE(6.111> 

111  FORMAT  1 1HA.26HVP3*V3*SPHI3  FOR  THIS  CASE) 

GO  70  13 

108  TANNP4— V3«CPMI3'SQRT<VINF3**2-IV3KPMISI««2| 

GO  TO  112 

110  TANNP4»V3*CPHl3/rORM VlNF3#*2»(V3*CPHI 51  **21 

112  DO  120  X1-1.L1 

E4(Xl ) • (RAOP3+D3 ( K1 1 •6v80. 2T4A4 I /A4 
XNU4«ATAK(  1  ./SORT  (£4tKl)**2-*l«)| 

XN4DEGC  XI ) "XNU4/R 
IFIVP3-V3»$PHI 31113. 113. 114 

113  PSI4«ATANITANKP4»-2.*XNU4*PI 
GO  TO  113 

114  PSU-ATANITANNP4)-2.*XNU4 
113  PSIADGtKl ) "PSI4/R 

V4(X1 )  "SORTI  VP3»*’2*VINF3»«2-2.*VP3*VINF3*C0S(  PS|4 1 ) 
CPH17-VINF3»SIN(PSI4)/V4(X1| 

SPH17«S0RT ( 1 •-CPHI 7**2 ) 

RPV4«RP3»IV4(X1)**2I/GCS 

TaNF7-RPV4*CPHI7*SPMI7/|1,-RPV4»(SPHI7««2>) 

IF(TANF7)116. 117.117 

116  F7IK1 )«ATANC TANF7I /R 
GO  TO  111 

117  F7IX1 !■ ( ATANI TANF7) -PI ) /R 

118  AS  1X1 l«RP3/( 2.-RPV4 ) 

A5AUIK1 )*A5(K1 )/ (92901000. *3280*) 
E3(Xn-SORTIURPV4-l.)*SPH17)«*2*CPH|7«»2) 
0MM(X1)-A3AU(X1)«(1.-E3(X1M 
THTSPCXll— F7(X1) 

T3P(K1)-*TIHE(A3(K1).GCSiE5IK1).F7(X1)I 
THTOT (XI  )•  TMT12>THT23*THT3P(X1‘ 

TIM0PIXn«TIM12^T23»T3P(Xn 

DJ(Xl4l)*D3(Xn4lOO. 

120  CONTINUE 

WRI TEI6.122 ) VB( l ) .PH( I  * J) »D(X)  .QM(K)  .ITMT23.F3.V3.VINF3.T23.03IKM 
1  »F7(X1)  .AJAU(Xl).E3(Xn,V4(Xl)  .XN40EGI X 1 )  »PSI4DG(  XU  .THTOT  (X1J  »  T  |M 
2 OP  I XI ) »QMft( K1 )fXlal.Ll> 


122  FORMAT (lHA.,4X.3HVa-.F10.0.3X»3HPH>»F7*2.3Xt2HD*»F7«1.3X»3KQM*»F8« 
13//2X.4HTH2>.4X.2HF3.6X.2HV5.7X,5MV1NF5.4X*3HT23i3X»2HOI»3X»2HF7*6 
2Xf2HA3*3X.2HE3.6X.2HV7.6X.2HN4*4X»4HP$l4.4X.4HTH0P.4X»3HT0P»3X.3MO 
3MM//(2F8.2.2F10.i.F7,2.F7.1,F«*2,2F7#4.F10.1.F3*2.2F8.2.F7*2.F8,3l 
4) 

13  CONTINUE 
STOP 
END 

fOATA 

4.6448E21  1.4077E16  U1408tl6  4.9032E11  3.34I0E11 
9.7T312C4  1.1493E03  2.0902E07  2.0013E07 
7.T290E14  7.9393E«6  1.3083E11  2.2893E11 
f  EOF 
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Vfl- 

50000. 

Program  3  Output  fi 

90.00 

THU 

n  41 

£1 

VI  TTOT 

OOCA  F3  A3 

£3  V3 

H 

722 

PS!  QH 

44*26  - 

133.74  0.6353 

0.5741 

106675.2 

92.2 

0. *  -152.78  0.3283  0 

.6464  91283.8 

6.74 

55.41 

52,50  0,18679 

44.26  - 

135.74  0.6353 

0.5741 

106675.2 

92.2 

100.0  -152.42  C.5305  0 

.6444  91689.0 

6.57 

55,41 

52.86  0,18863 

44.26  - 

135.74  C.6353 

0.5741 

106675.2 

92.2 

200.0  -152.07  0.5325  0 

•6425  92073.5 

9.40 

55.41 

53.79  0.1)040 

44.26  - 

135.74  0.4353 

0.3741 

106675.2 

92*2 

$00,0  -151.74  0.5346  0.6407  92438.8 

6.24 

55,41 

53.51  0.19209 

44.26  - 

133.74  0.6353 

0.5741 

106673.2 

92.2 

40S.0  -151.4?  0.5365  0 

•6389  92786,5 

6,09 

55,41 

53.61  0.19372 

44.26  - 

135.74  0.6333 

0.3: 

41 

106675.2 

92.2 

5C0.0  -131.11  0.5383  0 

.6373  93117.6 

5.94 

55,41 

54,10  0.19528 

V8- 

5OC0C. 

PH-  VO. 00 

Or  0.  OH*  0.18679 

?H23 

F6  V3 

VINF3  T23 

D3 

F7  A5  £5 

Vr  « 4  PS14  THOP  FOP  QHH 

62.79 

-69.99  209905. 1 

t 15243.2  21.40 

0. 

-91.66  0.5102  0.6391 

206S73,3  0.42 

80,55 

198,72 

92.24  0. 16411 

62.79 

-89.99  209905.1 

115243.2  21.40 

100.0 

-91.55  6.5114  0.6396 

<>C(*6$.e  0*39 

60.61 

198,60 

92.26  0,19430 

62.79 

-89.99  209905.1 

115243.2  24,40 

200*0 

-91.44  0.5125  9.6401 

208575.9  0,36 

80,66 

198.49 

92.25  0.28447 

62.79 

-69.99  209905.1 

115243.2  23.40 

300.0 

-91.35  0.5135  0,6405 

208658,4  0,34 

80.70 

191.40 

92.25  C. 18461 

62.79 

-89.99  209905.1 

115243.2  26.40 

400.0 

-91.27  0.5143  6.6408 

209731,2  0,32 

80,74 

198,32 

92*24  0.18474 

62.79 

-89.94  209905.1 

'1 

2243.2  28.49 

500.0 

-91.20  0.5151  0.6411 

208796.0  0,30 

80,76 

158,25 

92.24  0.18483 

V8- 

30000. 

J>M«  90.00 

0«  100.0  QMu  0.18663 

VH23 

F3  V5 

V1HF5  T23 

03 

F7  A3  65 

V7  N4  P$14  THOP  TOP  OHH 

67.29 

-89.22  210081.6 

114266.7  26.42 

0. 

~ 90.92  0.5121  0.6369 

208537,2  0,» 

80,65 

198,37 

92.27  0,18594 

63.20 

-89.22  210081.8 

114216.7  21.42 

100.0 

-90.80  0.5133  0.6374 

206646,5  0.*4 

00.91 

198.25 

92,27  0,18613 

63.20 

-69.22  210011.6 

114216.7  26.42 

200.0 

-90,69  0.5144  0,6379 

208741,4  0,37 

80,96 

196jlS 

92.26  0,18629 

63.20 

-89 .22  210081.6 

114286.7  28.42 

300.0 

-90.60  0.5154  0.6383 

206824,6  0,34 

•  1.00 

199. OC. 

92,26  0,18644 

43,20 

-89.22 

210081.6 

114216.7 

28.42 

400.0 

-90.52  0.5162  0.6366 

206896,0 

0.32 

81,04 

197,90 

92.25  0.10657 

63.20 

-89.22 

210081.6 

114266.7 

26.42 

500.0 

-90.45  0.5170  0.6369 

208963.3 

0.31 

61.06 

197.90 

92.25  0.1861K 

V6-  .50000.  PM- 

90.00 

0- 

200.0  OH-  0.19040 

TH23 

F3 

VS 

V1NF5  T23 

D3 

F7  A3  15 

vr  H* 

PS24  THOP  TOP  OHM 

43.59 

-86.48 

2 10249.7 

113367.5 

26.44 

0. 

-90.20  0.5139  0.6346 

206692.9 

0.4) 

61.13 

198.05 

92.29  0.18769 

43,59 

-88.48 

210249.7 

113367.5 

28.44 

100.0 

-90.07  0.5151  0.6353 

208803,2 

0.40 

61,20 

197.92 

92.20  0.10708 

63.59 

-88.46 

210249.7 

113367.5 

28.44 

200.0 

-69.97  0.3163  0.6337 

206898.8 

0.37 

01.25 

197.02 

92.27  0.18105 

43.39 

-88.48 

210249.7 

P 3267.5 

28.44 

300.0 

-89.18  0.5172  0.6362 

208982,6 

0.35 

81.29 

197.73 

92.27  0.16819 

63.53 

-88.48 

210249.7 

L,*S07.$ 

28.44 

400.0 

-89.80  0.5181  0.6365 

209056,6 

0.13 

81.34 

197.65 

92.1*6  0.188)2 

43.59 

-31,48 

210249.7 

113*47.5 

28.44 

500.0 

-69.72  0.8109  0.6361 

209122.4 

0.31 

01.37 

197.57 

92.! 6  0,18844 

VI-  50000.  H- 

90.00 

0- 

300.0  OH-  0.19209 

TK23 

f  5 

V3 

VINF3  T23 

03 

FT  A5  83 

vr  N 

PSI4  THOP  TOP  QHH 

43.97 

-87.77 

210409.9 

112483.5 

28.46 

0. 

-89.50  0.3136  0.6327 

208841.1 

0.44 

81,41 

1.7# 73 

92.30  0.109)7 

63.97 

-17.77 

210409.9 

112483.5 

100.0 

-•9.38  0.5169  0.6333 

208952.2 

0.41 

•1.47 

1*7.61 

92.29  0.18956 

63.97 

-87.IT 

210409.9 

112483.3 

20.46 

200.0 

-89.27  0.5180  0.6337 

209048.6 

0.38 

•1.5) 

1*7.30 

9^,20  0.1897) 

43.97 

-17.77 

210409.9 

112483.5 

28.46 

300.0 

-69.18  0*5190  0.6341 

2091)3.0 

0.36 

01,58 

1*7.41 

92.28  0.16988 

43.97 

-87.77 

210409.9 

112483.5 

20.46 

400.0 

-69.10  0.3199  0.6343 

209207,6 

0.34 

•1,62 

1*7.33 

92.27  0.19001 

43.97 

-87.77 

210409.9 

112413.3 

28.46 

500.0 

-19.02  0.5207  0.4248 

209273,9 

0.32 

•  1,65 

1*7,2* 

92.27  0.I9OU 

VB-  50060.  PH- 

90.00 

D- 

400.0  OH-  0.19372 

TM23 

F5 

VS 

VIKF3  T23 

03 

FT  A5  65 

V7  W 

PSI4  THOP  TOP  OHM 

44.34 

•67,08 

210562.9 

111632.7 

28.48 

0. 

-81.13  0.5172  0.6306 

208982.4 

0.44 

81,60 

197.43 

92.31  0.19096 

“  64.34 

-87,01 

210542.9 

111632.7 

20,48 

100.0 

•88.71  0.3185  0.6313 

209094,3 

0.41 

01,74 

197,31 

92.30  0,10117 

44*34 

-87*08 

210562.9 

111632.7 

28.48 

200.0 

-66.60  0.3197  0.6318 

209191,4 

v.)9 

61,10 

197.20 

92.29  0.19134 

44.34 

-67,16 

210362.9 

111632.7 

26,46 

300.0 

-88.51  0*3207  0,6322 

209276,4 

0,36 

01,05 

197.10 

42.28  0.19149 

64.34 

-87.08 

210562.9 

111632,7 

28.46 

400.0 

-08,42  C.5216  0.6326 

209351.5 

0.34 

01,09 

197.02 

92,28  0.19167 

64.34 

-87.00 

210562.9 

111632.7 

26.46 

300.0 

-88.35  0.5224  0.63)0 

209418,4 

0,32 

01,9) 

194.95 

92,26  0.19176 

VB-  50000.  PH- 

90.00 

0- 

500.0  QM>  0.19328 

IH23 

F3 

V5 

VINF5  723 

03 

FT  A5  63 

VT  M4  PSI4  TMOP  TOP  QH* 

64.69 

-86.42 

210709.0 

110113.1 

28.50 

0, 

-18.19  0.5168  0,6269 

209117.0 

0,45 

01.94 

1  17.14 

92.32  0.14253 

64.49 

-86.42 

210709.0 

110913.1 

28.50 

100.0 

-68.04  0,5201  0.6295 

209229.7 

0,42 

02.U1 

197.01 

92.31  *.1427 2 

64,69 

-84.42 

210709.0 

110813.1 

2I.SC 

200.0 

-£7,93  0,5213  0.0300 

209)27.5 

0.39 

02,06 

196.91 

42.30  0,14289 

64,69 

-66.42 

210709.0 

110813,1 

28.50 

300.0 

-r  ,86  0.5223  0.6304 

209413.2 

0,57 

02.11 

196.01 

92.29  C. 19104 

64,69 

-66.42 

210709.0 

110813*1 

28,30 

400.0 

•87.77  0.52)2  0.6306 

209468.8 

0.35 

02,15 

146.73 

42.29  0.19317 

64.49 

-16.42 

210709.0 

110113,1 

28,50 

500.0 

-87,70  0.5240  0.6312 

209556.2 

0.33 

02.19 

196.65 

92.20  0.19329 

i 

i 


i 
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Program  4 

3303  0.8.4500  67-012.LT  HYERS.AFlT-SK 

NAP 

8I6FTC  KAIH  M»4/2.XR7,ftO0E« 

C  DOUBLE  ASSIST  AT  VERU? 

TlHE(A.O.E.T)>(2.*ATAMS<SaTI(l.-£)m.4E)>*TANtT*R/2.))“E*S«R»«l.- 
lE**2)*SIN(T»R)/(1.4E»C0SIT«R))>*SQRT(CA*»3)/(i>/864Cia. 
REA0(5»1)«CS.GCP1»GCP2.RP1.RP2.VP1.VP2.RA0P1.RADP2 
1  FOIIMAT(5E10*0/4E1G.O> 

C  DIMENSION  FOR  ORE  GREATER  THAN  REQUIRED 

DIKEHSION  V£(162).PH(162.10) .F0U62.10)«VOt  162.10) .AI162.10I.ADU6 
12.10) >ET(162.10).VIMFP1 (162) .311 162) 

DIMENSION  02(25) .£4(25).XN4DEG(25) «PS14BG(25)»V4I23) »F5(25)*ASI25) 
1  .A5AU(25).E5(2S>. OH2(25). THTOTI25) .TIH0PI231 
C  V8(l)  MOST  BE  GREATER  THAN  ESCAPE  VEL  ■  50700  FT/SEC  FOR  EARTH 
VRU)  >58030. 

Pl«3.14159 

R-PI/I80. 

N«33 

N>1 

L2*21 

L>1 

TAUP2>2.e?l»S0«T(RP2»«3/GCS)/t4<k''0. 

C  CALCULATION  OF  DEPARTURE  ORBITS 
DO  2  l*l.H 

VINFPK  I  )<'SQftT  (VB(  ll*»2-2.*GCPl/RA0Pl) 

PH(l.l)*90. 

DO  3  3*1  .N 
PHR*PHI1.3)»R 

JF(VINFP1(II**2-(VP1*C05(PHR))»»2)24.24.2S 
24  J1(I)>J-1 
GO  TO  29 

23  TANETF>VP1*C0S  (PHR)/(S0RT  (VINFPKI  )*«2-(VPl»C0S  !PKR))*«2)) 
3l(l)*3 

ET(I.J)-(PI-ATAN(TANETP))/R 
A(!.3)*Pl/2.-£T(1.3)*R*PHR 
AD(I.3)*A(I.3)/R 
IFIPKt 1.31-90. 119,20. 19 

20  VO(I.J)-VPI-VINFPKI) 

GO  TO  21 

19  V0C.J)>VINFP1(I)*SIN  (A(  1 .  J) )  /COS  (PHRI 

21  RPVO-RP1»IVO(I.J)«2)/GCS 

TANFO*RVVO»COS  (PHR)*SIN  (PHRI/(1.-RPV0»SIN  IPKR)**2) 

F0(I.J)-(ATAN  (TANFOI-PD/R 
PHI  t.3H5*PH(t. 31-2. 

3  CONTINUE 

29  V0( 1+1 )*VB( 1 1+1000. 


2  CONTINUE 
WR1TEI6.4) 

4  FORMAT 1 1H1 .30X.21HEARTH  DEPARTURE  ORBIT) 

DO  6  t-l.H 

J2-J1U) 

WR:TE(6.9IVO(II.VINFPllI>.(PH(I.3>,FOtl.3)»VOt!»3).ADtI.3).ETII.3) 

1.3>1.32) 

5  FORMAT(lHA»25X.3HVB**F10.0.33X*7HVINFPl*.Fl0.0/17X.2HPH»18X.2HF(‘»l 
18X.2HV0.16X. 1HA. 19A.2HET// ( 5F20.2 ) ) 

6  CONTINUE 

C  CALCULATION  OF  PRE-AS5IST  ELLIPSE 
DO  13  I«1.H 
32*31(1) 

DO  13  3*1.32 
0»0* 


Al«RPl/(2.-RPl*(V0( I »3)**2)/GCS) 

AlAU"Al/( 9290 1000**5280.) 

El-SORT  (((RPl«(V0II.3)**2)/GCS-l.)*SIK  (PHI  I *3)*RI )**2»C0S  (PHII. 
13)*R)**2) 

RPER1>A1*( 1.-E11 
IFIRPER1-RP2) 15.15.16 

16  NRITE(6.17) 

17  FORMAT! 1HA.31HVB  INSUFFICIENT  TO  INTERCEPT  P2) 

GO  TO  13 

15  Vl-SCRT  (GCS*(2./RP2-1./A1)) 

SPHIl»SO®T  ((Al**2)*U.-Ei**2)/(RP2««2.*Al-vV2))> 

CPHI1«SQRT  (l,-SPHll«*2) 

TAN?>(RP2*IVl*»2)/GCS)*CPHIl*SPHll/(l.-(KP2*l\'l**2|/6CS)*SPHll«*2 

1) 

1F(TARF1I4L'.40.M 

40  F1*ATAHCTANF1|/K 
GO  TO  42 

41  Fl-IATAR  (TANF11-PD/R 

42  THT12““1-F01 I .31 

TIH12«TIHE(A1.GCS.E1(F1)-TIN£( A1.GCS.E1.F0I 1.3) ) 

C  CALCULATION  OF  HYPERBOLIC  PASS  TRA3ECT0SY— FIRST  PASS 
VINF»SSRT  (V1**24VP?.**2”2»*V1*VP2*SPHI1I 
A2»GCP2/«¥IHF»*2) 

10  E2>(RADP24D*6080.27tA2)/A2 
XNU*ATAM ( 1 ./SORT ( E2  *»2-l . ) ) 

XNUDEa»XNUFR 
1F!VP2“V14SPHI118.7 ,9 

7  WRITEI6.14! 

14  FORMAT ( 1KA.26HVP2*VI*SPHI 1  FOR  THIS  CASE) 

GO  TO  13 

8  TAHNPH—V1KPHI1/S0R7  IVINF»*2-tVi*CPHU)*«2l 
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Program  4  (Contd.) 

44  PSI*ATAN!TANNPHI-2.*XIIU+PI 
CO  TO  30 

9  TAKHPH*V1«CPH1 1/SORT  IVIKF**2-IV1*CPHI1I**21 
A3  PSI-ATAN  ITAKXPHI-2.*XNU 
50  PS1DE0  OP51/R 

V2  *S«RT  !VP2**2+V!NF**2-2.*VP2*VINF*C0S  IPS! II 
CPHI3*V!NF*SIN  (P3I1/V2 
SPHIS*SCRT  (l.-CPHiJMI) 

RPVS*RP2*iV2  **2I/«CS 

TANr5^IPVS*CPHl»»SPHt3/<l.-«PV$»SPH13»»2l 

IP1TAK?S)43.44«4* 

45  F3  •IATAJUTARF3))/R 
CO  TO  AT 

46  FS  *!ATAM  (TANF3I-P1I/R 
AT  A3  ■IRPl/li.-RPVSlI 

AA  *AS  JI92901000.P3240.I 

E3  *SCPT  «HR»V3-Ss!*$Pn;siM2ACPHt3«21 

AP3*A3*il.+E3) 

lPiAP3-RP2l214.215.21S 

214  WRITEI6.217) 

217  FORMAT  IlHA.3Xf34KAPKELI  OR  TOO  SKAU  FOR  SCCORO  PASS) 

GO  TO  13 

215  OM*AA*!l.-E31 
T2P*-TIKEIA3*CCS*E3*F3I 
T-T2P+TIM12 

TAUP*2.*P!*SeRriA3»*S/6CSl/46400. 

00  200  HP1*1.5 
00  200  N1P*1.? 

IF1ABSIFI.OATINSP1»TAUP-FI.OATINPU*TAUP21-.1I201. 200.200 

200  CONTINUE 
0*0+5. 

IF( 0*5000* 110 i 10*202 

202  WRITEI4*209IVB1I I *PH( I • J) 

203  FORMAT  I  IKAOX-TiMNO  MORE  TRAJECTORIES  FOR  YB*,F10.1*2X.7MAN0  PH* 
1.F4.2.2X.13K7EL0W  3000  Nil 

CO  TO  13 

C  CALCULATION  OF  SECOND  ASSIST  TRAJECTORY 

201  V1NF3"SQRT!V2**2+VP2**2-2,*V2*VP2*SPM13| 

A4*CCP2/VIHF3**2 

02111*0. 

00  212  K2*1*L2 

£4 ! K2 1* ( RADP2+02 1 K2 1*4010. 2T+AA I /AA 
XRUA*ATANI1./S0UT(EA(K21S«2~1.]I 
XN4DECI X2 1 "XKUA/R 
IP I VP2-V2+SPH13) 204 *203*204 

204  TARNP4— V2»CPHI3/SQRTIVINF3**2-!V2*CPHIS1**2I 


PC14.ATANITANRPL1-2  .»XN'J4+P  I 
CO  TO  204 
243  VR1TEI4.20TI 

20T  FORMATI 1HA*24HVP2*V2*SPHI3  FOR  THIS  CASE1 
CO  TO  212 

204  T/.NNP4+  V2*CPHIl/SQRT!VIKF3*+2-!V**CPNI3l**21 
PSI4*ATAN!TANNP41-2**XKU4 
204  PSIADOIR21-PSIA/R 

V4IX2l*SORTIVP2**2+vmF3**2-2.*VP2*VlNF3*COSIPSI4»l 
RPVA.RP2* ( VA (K2 I **2 I /CCS 
IFIRPV4-2. 1219. 222. 222 

222  L*X2+l 

CO  TO  223 

21?  CPHI3*VINF 1*SIN!PSI41/V4(K21 
SPNI3*SaRKl.-CPHl3**2l 
TANF5*«PVA^PHI5*SPHI9/I1.-RPVA*SPHI5**2I 
IF5PSI4 I 327*324 *324 
327  F3IX2I*IATAR(TAKF31+PI1/R 
00  TO  211 

324  IF CTARF51 203*210*210 

209  FSIR2I*(ATAN(TANF3II/R 
CO  TO  211 

210  F5(R21*(ATAM(TARF3I-PI1/R 

211  ASIX21-RP2/I2.H1PV4I 

A3AUIX2 I *A3(K2 1/192901000**3240*1 
E5IK21*SO*TIIIRPV4-l.t*SPHI5l**2+CPHI3**21 
OM2 1 K2 1  *A3  AUt  K2 1  *  1 1 *-E3 ( K2 1 1 
1FIPS1AI329. 310*330 
32?  THTOTIX21-THT12-P5IX21+360. 

,TIM0PIK21*TIM12+PL0ATINSP1*TAUP“TIHE!A3!F.2I  .OC3.E5IX2I  *F3tR2)  l+TAO 
IP 

SO  TO  223 

330  THTOTIR21«THT12-F3ltt2l 

7IHOPIX2l*TlM2+FlOAT!ll4PI*TAUP-7IME!A5IX2l,eC5.E5iX21*P5!X2>l 

223  02 1X2+1 l«D2IK2 1+10000* 

212  CONTINUE 

WRlTEI6.2131VBtII.PH! I*J}*D*TAUP2*HP1*M3P*THT12*P1*A1AU*E1*V1*TINI 
12,T*F3.AA*E3.V2.XNU DEO. P5I0E0.TAUP.cn, 1021X21. V!NF3.VA(X2l*XNAOEGt 
2K2)*P5IA00(K2I,F5(K2I*THT0T(K2I •TIM0P,C2),A3AU!X21 *E3(K2}*GM2tK21 * 
3X2*L.L21 

213  FORMAT! lHA>SX*3HVe*iF10*l*5R*3HPH*,F10.2*5X*4H0OCAl*,F5*2*5X*4HTAU 
lP2**F7.1v5X.4KHPl+.14*9X,4HN5P**IS//2X.4HTHl2.3*i2KPl*4X.2HAl.3X,2 
2HE1*4X*2HV1*4X*3HT12*9X*1HT*4X*2KF3,4X*2HA3*3X*2KE3*4Xi2KV3»7X*1HN 
3,7X*3HPSt«AX.AHTAUP.AX.2HCN//2F4.2,7FT.A,F10*l*2FT.l,F0.2*2F7*A>Fl 
40*1 *2F4*2.FT*1.F4*3///3X.3HOOCA2.9X*3MV INFS *TX.2HV3*0Xt2HNA,4X. AMP 
3SI4*6X*2HF9*4X*3HTHT0T*lX*3HTlK0P*3X*2HAS*4X*2K£5*6X*tH0M2//IF9*0* 


42Fll*l*AF9*2*F0*l,2Fe*A,F9*SII 

0*0+13* 

00  TO  10 
13  CONTINUE 
STOP 
END 


4.65S4E21  1.A0T7E14  1.14ME14  A.9032EU  3.5A40EU 
t.TTSieOA  1.1A93E03  2.C902E07  2.001SE07 

4ECP 
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V0*  51CCC.0  PM- 


Program  4 

CGCA1*  145.00  TAU?2° 


Output 

214.tr  M 


MP1*  5 


NSM  I 


IH12 

M  A1 

El 

VI 

T12 

T  ?* 

A3 

E3 

V3 

H 

PS! 

42,46 

•137.54  0.6276 

C.5933 

IC5611.0 

54.4 

90.5  -152. 6'.  0.5287  0.6603 

91367.9  6. 

14  52.64 

CCCA2 

VIKF3 

VJ 

M4 

P514 

F5 

7H707 

riMOi* 

A3 

» 

0M2 

0. 

67782.7 

7647C.3 

6.38 

99.88 

-163.53 

205.99 

1219.0 

0*4645 

0.7303 

0.12527 

iuce»;«j 

67)32.7 

79594.3 

5.02 

42.60 

-161.62 

204.09 

1212.9 

0.4757 

0.7163 

0.13495 

2000.00 

67782.7 

6164C.1 

4.14 

44.37 

-160.28 

202.74 

1212.8 

0.463? 

0.75*9 

C. 14173 

3000,00 

67782.7 

63CI2.0 

3.52 

45.60 

-159.30 

201.76 

1212.8 

0.4896 

0.7002 

C,i\H77 

4CCO.CO 

67T62.7 

84152.5 

3.06 

<6.52 

-151.55 

201.01 

1212.8 

0.4941 

0.4952 

0.130*2 

SOCO.OO 

67782.7 

84978.3 

2.71 

47.22 

-157,95 

2C0.41 

1212.8 

0.4977 

0.6912 

0.15U? 

*300.00 

67782.7 

65634,8 

2.43 

47.78 

-157.47 

199.99 

1212.8 

0.5906 

0.4881 

0.1S614 

/COO. 00 

67782.7 

86169.0 

2.20 

48.23 

-157.07 

199.53 

1212.8 

0.3030 

0.6850 

0.15818 

£000.00 

67782.7 

86612.3 

2.01 

41.61 

-156.74 

199.20 

1212.8 

0.5051 

0.6834 

0.15990 

9000.CC 

67782.7 

86985.9 

1.66 

48.93 

-156.45 

191.91 

12*2,8 

0.5068 

0.6CJ* 

9cl*13* 

1CCC0.00 

67782.7 

87305.2 

1.72 

49.20 

-156.21 

198.67 

1212.8 

0.5083 

0.680V 

0.16762 

11CCO.OO 

67782.7 

67581.1 

1.6Q 

49.49 

-:55.99 

191.45 

1212.8 

0.5096 

0.6787 

0.1 6372 

12000.00 

67782.7 

87821.9 

1.50 

49.44 

-195.80 

198.26 

1212.8 

0.5108 

0.6776 

0.16468 

130C0.C0 

67702.7 

88033.9 

1.41 

49.82 

•155.64 

191.10 

1212.8 

0.5118 

0.6766 

0.16554 

14000.00 

67762.7 

•8222.1 

1.33 

49,97 

-155,49 

197.95 

1212.8 

0.5127 

0.675* 

0.16630 

16000 .UO 

67762.7 

88390,1 

1.26 

50.12 

•153.36 

197.81 

1212*8 

0.5135 

0,*748 

0.14698 

UCOO.QO 

67782.7 

68541.1 

1.2C 

50.24 

-155.23 

197.69 

1212.8 

0.5143 

0.*741 

0.16760 

27CCO.CO 

67782.7 

68677.5 

1.14 

SC. 36 

-153.13 

197.58 

1212.8 

0.5149 

0.6734 

0.16816 

16000.00 

6775:2.7 

68101.4 

1.09 

50.47 

-135.03 

197.48 

1212.8 

0.5155 

0.6728 

0.t6867 

19010.00 

67782.7 

61914.3 

1.04 

30.56 

-154.99 

197.19 

1212.8 

0.51*1 

0.6723 

0.16913 

20000.00 

67762.7 

89017.8 

0.99 

50.65 

-154.85 

197.31 

1212.8 

0.5166 

0.6718 

0.16956 

ve- 

SlUOO.O 

PM* 

90.00 

OOCAt*  164C  #  CO 

TAUP2* 

224.3 

MPi» 

2  MSP*  3 

TH12 

FI  XI 

fct 

VI 

712 

7  F3 

A3 

VS 

M 

PS! 

42.44  - 

137.54  0.6276 

C.5933 

lcsni.o 

54.4 

90.5  -149.16  0.5520  0.6405 

95438 

•  9  4. 

41  56.08 

CCCA2 

VIKF3 

V5 

K4 

PS  14 

F5 

7HT07 

TIWP 

A5 

ES 

QM2 

0. 

67782.7 

80430.4 

6.31 

43.99 

-161.09 

203.54 

539.4 

0.4789 

0.7125 

0.13766 

1000.00 

67782.7 

8360C.8 

5.02 

46.03 

-158.94 

201.40 

539.3 

0.4917 

0.6978 

0.14862 

2000.00 

67782.7 

65669.3 

4.14 

47.81 

-157.45 

199.90 

539.3 

0.5008 

0.6679 

0.15*27 

TAU?  CM 
140.2  0.17940 


TAUP  CM 
149. A  0.19149 


30CO.C0 

67782.7 

87124.0 

3.32 

49.04 

-156.35 

148.81 

539.3 

0.50  n 

0.6(10 

0.161*0 

4000.00 

07782.7 

88202.2 

3.04 

49,96 

-155.50 

197.9* 

539.3 

0.5126 

0.6757 

0.16622 

SOCO.OO 

67782.7 

89033.3 

2.71 

50.66 

-154.84 

197.30 

539.4 

0.5167 

0.6717 

0.16*65 

6000.00 

67782.7 

89693.2 

2.43 

51.22 

-154.30 

196.7* 

539.4 

0.5200 

0.66(5 

0.1725* 

7000.00 

67782.7 

90230.0 

2.20 

51,67 

-153.85 

196.31 

539.4 

0.5228 

0.6(5* 

0.176(6 

8000,00 

67782.7 

90675.2 

'2.01 

52.05 

-153.41 

195,94 

539.4 

0.5251 

0.6657 

0.17*5* 

90L0.C0 

67782.7 

91050.2 

1.86 

52.37 

-153.16 

195.62 

539.4 

0.5270 

0.661* 

0.17*10 

looco.oo 

67782.7 

91370.6 

1*72 

52.64 

•152.88 

195.34 

539.4 

0.52(7 

0.6605 

0.17*61 

11000.09 

67782.7 

91*47.3 

1 .40 

52.87 

-152.64 

195.10 

539.4 

0.5502 

0.65*0 

0.190(5 

12000.00 

67782.7 

91888.9 

1.50 

53.08 

-152.43 

194.89 

539.4 

0.5515 

0.657( 

0.1*150 

130C0.00 

67782.7 

92101.5 

1.41 

51.26 

-152.24 

194.70 

539.4 

0.5527 

0.6567 

0.192(5 

14000.00 

67782.7 

92290,1 

1.33 

53.42 

-152.01 

194.54 

539.4 

0.5517 

0.6556 

0.1*570 

15009.00 

67782.7 

92458.5 

1.26 

53.5* 

-151.93 

194.39 

539.4 

".5567 

0.6550 

0.1(666 

16900.00 

67702.7 

92609.8 

1.20 

53.69 

-151.79 

194.25 

539.4 

-.5555 

0.6565 

0.16516 

17000.00 

47782.7 

92746.5 

1.14 

53.80 

-151.67 

194,13 

531.4 

0.5565 

0.6556 

0.1(578 

11000.00 

67712,7 

92870.6 

1.09 

53.91 

-151.5* 

194.02 

539.4 

0.5570 

0.6550 

0.1(655 

14000.00 

67782.7 

92983.8 

1.04 

54.00 

-151.4* 

193.92 

439.4 

0.5176 

0.6526 

0.16689 

20000.00 

67782.7 

93007.4 

0.99 

54.09 

-151.36 

193,82 

539.4 

0.55(2 

0.651* 

0.1(752 

»#■ 

5100C. 0 

PH* 

9C.00 

OOCAl*  4960.00 

TAUP2* 

22*.  5 

MPi« 

.5  MSP*  7 

THU 

FI  A1 

El 

VI 

T12 

7  FJ 

AS 

83 

VS 

N 

PSt 

42.46  - 

137.54  0.6276 

0.5933 

105811.0 

54.4 

90.6  -145.08  0.5780  0.6214 

99*3*,*  2. 

72  54.47 

00CA2 

V1NF3 

V3 

N4 

PSI4 

F5 

THTOT 

TIMOP 

A5 

E3 

0M2 

0. 

6?m.  r 

84362*7 

6.38 

46.71 

-158.38 

200.8* 

1212.8 

0,4931 

0.69*1 

0,19146 

1000.00 

67782.7 

87582.2 

5.02 

49.43 

-155. 49 

19I«*S 

1212.1 

0.509* 

0.6787 

0.1*372 

2000.00 

67712.7 

89663.3 

4.14 

51.20 

-154,32 

194.78 

1212.9 

0.3199 

0.6686 

0.1722* 

3000.00 

67782.7 

91124.1 

3.52 

52.43 

-153,09 

195.55 

1112.9 

0.527* 

0.6615 

0.17853 

4000.00 

67782.7 

92205.5 

3.06 

53.34 

-132,15 

19*.61 

1212.9 

0.5333 

0.6562 

0.18)32 

50(0.00 

67782.7 

93038.2 

2.71 

54.05 

-151.41 

193.17 

1212.9 

0.3379 

0.4922 

0.18710 

6000.00 

67782.7 

93699.0 

2.43 

54.61 

•150,81 

0.43*26 

1212.9 

0,5*17 

0.6*89 

0.19019 

70C0.00 

67782.7 

94236.3 

2.20 

55.06 

-150.31 

192.77 

12U.t 

0.5**l 

0.6*63 

0.192*7 

8000*00 

67782.7 

94681.5 

2.01 

55.44 

-149.89 

192.33 

1213.0 

0.3*7* 

0.6*42 

0.19*78 

9000.00 

67782.7 

95056.6 

1.86 

55.76 

-149,53 

191.99 

1213,0 

0.5*9* 

0.6423 

0.19659 

10000.00 

67782.7 

95376.8 

1.72 

56.03 

•1*4.22 

191.48 

1213.0 

0.531* 

0.64CI 

0.19113 

11000.00 

67782.7 

95653*4 

1.60 

56.26 

-141,9* 

191.42 

1213.0 

3333 

0.6395 

0,199*8 

12000.00 

*7712.7 

85894.8 

1.50 

56.46 

-1*8.72 

191,18 

1213,8 

0,35*8 

0.4383 

0,200*6 

U0C0.0O 

*7782.7 

96107.2 

1.41 

56.64 

-1*8,51 

190.97 

1213.0 

0.55*1 

0.4373 

0.20170 

14000,00 

67782.7 

86295.5 

1.33 

56.80 

-1*1.33 

190.79 

it:*. o 

0.3572 

0.63*4 

0.202*4 

15000. CO 

67782.7 

96463.7 

1.26 

56.95 

•1*8.16 

190.62 

1213.0 

0.3583 

0.6359 

0.203*7 

16000.00 

67762.7 

96614.8 

1.20 

57.07 

•1*8.01 

190.*7 

1213.0 

0,5393 

0.6348 

0,20*23 

IToOO.OO 

67782.7 

96751.3 

1.14 

57.19 

-1*7.17 

190.33 

1213.0 

0.5*01 

0**342 

C.20491 

18000.00 

67782.7 

9*875.2 

1.09 

57.29 

-1*7.75 

190.21 

2213.0 

0,5*09 

0.4336 

0,20553 

19CC0.C0 

67782.7 

8(  988,2 

1.04 

57.39 

-1*7.43 

190.09 

1213.0 

0.5*1* 

0,4330 

0.20*10 

20000,00 

67782.7 

97091.6 

0.99 

57.41 

-1*7.53 

189.99 

1213.9 

0.9*23 

0.4325 

0.206*2 

TAUP  CM 
160.9  0.21810 


MO  MCAE  TUAJECTCAIES  FOR  V6»  51000, C  AMO  PH-  90.00  BELOW  5000  Nf 
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Program  5 

(JOS  0.2S.7000  67-012*LT  KYERSi AFIT-SS 

*19 JOS  MAP 

(IBFTC  MAIN  MJA/2.XR7 

C  TRIPLE  ASSIST  AT  VENUS 

TIKEtAfG,E,T>»t2«*ATAN(SORT<(l.-E»m.+eSl«TAKIT*H/2tM-S»SOam.- 
lEM2)«5lNCT»R)/<1.4£*C0StT»R)))#S0RT<<A»»S)/6l/fSA6i. 
READI5.1)GCS.GCPI.GCP2.RP1.RP2.VF1.VP2.RADP1.RACP2 
■  1  FORMAT! 3E10.0/AE10.0) 

C  DIMENSION  FOR  ONE  GREATER  ,:iAX  REQUIRED 

DIMENSION  VBI 162) <PM 162*10) «F0( 162» 10) *V0( 162*10) *A 1162*10) *AD! 16 
12*10) »ET( 162 *10) iVINFPl! 162) *J1 (162) 

DIMENSION  D4I29)*E6(2B)«XN60EG(25) *PS160G(23)*V6(2S)*F7(2S.**A7t2S) 
1.ATAU12S).E7I2SI.ONA(2S).THTC2(28I.TIK02(25) 

C  VB(l)  MUST  BE  GREATER  THAN  ESCAPE  VEIL-  S67C0  FT/SEC  FOR  EARTH 
VB!1) >3(000* 

FI-J.1A159 

R«PI/1I0. 

M*SS 

N«1 

12*21 

L«1 

TAUP2o£*«Ft«SQRTIRP2P*3/GCS)/(6A00* 

C  CALCULATION  OF  DEPARTURE  ORBITS  • 

DO  2  I-l.M 

V1NFPK 1  >»SORT  (V81 1  )»»2-2.PGCP’/RAOPl ) 

PHU.D-90. 

DO  3  J»l*N 
PHR*PH(I.J)tR 

IFIVINFP1(I)R*2-(VP1«C0S(PHR))«»2)2A.2A*23 
2*  JliI)«J-l 
GO  TO  29 

23  TANETP-VP.v«COS  (PHR)/«SQRT  IV1NFP1  (I )*«<2'  ( VP1«C0S  1PHR,)«*2)) 
J1ID-J 

ETU.J)»!PI-ATANITANETP))/R 
AC I*JI«PI/2#-ET( I* J)*R*PHR 
AO(I.J)-A!I.J)/R 
IF(PHCC*J)-9U.«  119*20*  19 

20  VOCt*JI*VPl-VINFPl<I) 

GO  TO  21 

19  V0CI.J)*VINFP1CII«5IH  (AU.JD/COS  IPHRI 

21  RPV0>RP1*(V0(I*J)«*2E/GCS 

TANF0.RPV0»C<JS  (FHR)»SIN  (PHR)/(1.-RPVG«SIN  (PKR)»»2) 

F0C1.J)«CATAN  CTARFOI-PD/R 
PH(I.J+l)*PH(I.J)-2. 

9  CONTINUE 

29  V0U9l)«VB(  1)«1000. 


C 


2  CONTINUE 


MC1TEC6*A> 


A  FORMAT ( 1H1 *S0X.21HEARTH  DEPARTURE  OR61T) 


DO  6  I"1*M 
J2-J1U)  . 

WRlTE(S.5)Vem.VlNFPllM»(PHIl.J).F0(I.Ji.VP(I.J).A0M!J}«Crtl.J) 

l*Jal»J2) 

S  FORMAT! 1HA*2JX*1HVB>*F10.0*33X*THV1NFP1«.F10.0/27X*2HPH*1«X.2HF0.1 
2BX»2HV0*16X* lHAr 19X*2KET//(3F20»2)) 

S  CONTINUE 

CALCULATION  OF  PRE-ASSIST  ELLIPSE 
DO  13  I*1*M 
J2*J1II) 


0*0* 

AlaRPl/C2.-RPl«CVSCI*J)M2)/GCS) 

A1AU*A1/(92901000.»32B0.) 

El-SORT  ( ( (RP1*IVS(  I*J)**2)/GCS-1.)*SIN  IPHI  I*J»RR))»2K0S  IPHII. 
1J)*R)**2) 

RPERlaAl**(  1.-E1) 

IFCRPER1-RP2 119*19*16 
IS  WRITEIS*!?) 

17  FORMATI1HA.S1HVB  INSUFFICIENT  TO  INTERCEPT  P2) 

GO  TO  IS 

18  Vl»SORT  CGCSM2./RP2-1./A1)) 

SPMU«SQRT  (IAl»t2)»(l.-El»*2)/IRP2»(2.»Al-*P2>») 

CPHSlaSQRT  (l.-SPHIi»»2) 

TAMFl»IRP2»lVl*»2)/8C8)KPHIlRSPHH/Il.-lRP2»CVl«92)/GCS)»SPMIlPR2 

1) 

IFCTANFDA0.AO.Al 
AO  Fl-ATANITANFD/R 
BO  TO  A2 

A1  Fl-IATAM  ITAKFl)?PI*/R 
A2  THT12"Fl-F0tl*J)  ,  ' 

TIM12»TIMIIA1*GCS.E1.F1)“T1K£I A1.GCS.E1.FOI I . J>) 

CALCULATION  OF  HYPERBOLIC  PASS  TRAJECTORY— FIRST  PASS 
VI KF-SQRT  CVT-»2*VP2**2-2 *«V1»VP2»SPM 1 ) 

A2*SCP2/CVINf!  *2) 

10  E2e.IRADP24D«ik"»«2T»A2)/A2 
XNU»ATAMI1*  FVC!«T\T2**2-1») ) 

XMUOCO.XNU/R 
1FCVP2-V1«SPM:«».T.9 
T  VRITEI6.1A) 

1A  FORMAT < 1HA*26HVP2»V1-IBHI1  FOR  THIS  CASE) 

GO  TO  IS 

8  TANNPH,’”Va*CPHIl/SQRT  1VINF««2-(V1*CPMH)»#2) 
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44  f>jjI|gT|tj(TANHPH)-2.4XHU4pt 


Program  5 


9  TAHS»H*V1»CPHI1/S0RT  (V1HF«»2-(V1»CPHI1)»»2) 

43  P5I«ATAN  ITANMPH)-2.«KU 
90  P5I0CG  "PSJ/R 

V2  *SORT  IVP20«2*VIKF»»2-2.*V?2»VIKF«COS  IPStl) 
CPH13*V!KF»SIR  IPSI1/V2 
SPH13»S«r  (1.-CPHI3M2) 

RPV3*RP2»<V2  8C21/GCS 

TANF3»RPVS*CPHI3»SPH!3m.-RPVS»SPHIS«*2) 

1F(TANPSI43*46*44 

43  F3  *lATANfTAHFJ))/R 
00  TO  47 

44  F3  *(ATAN  <TANFSt-Pl)/R 
47  AS  «<RP2/l2.-RPVSn 

AA  «AS  /( 92901000 **5280o) 

E3  •SORT  If  IRPV3-l.)*SPHm«»2*CPH13<»*2) 

AP3«A3»(1.4E3) 

IFIAP3-RP2I 216*213.213 
214  WRITE16.217) 

217  F0RKAT<1HA.3X.3AHAPHELIC«  TOO  SHALL  FOR  SECOND  PASS* 
GO  TO  13 

213  Qn»AAWl,-E3) 


(Contd.) 


T2P«-TH»E1AS*GCS«ES.FS> 

T-T2P+TIM12 

TAOP-2. *PI»SQRT <  A3*»3/GCS ) /« 4400 . 

DO  200  MPIMtS 
DO  200  N5P“1*9 

IF(  ABS(  FLOAT  ( NSP 1  PTAUP-FLQATt  HP  U«TAUP2)-.  1)201  *200*200 
ZOO  CONTINUE 
212  D«IHS. 

IF(0-S000.tl0. 10.202 

202  NRITE(G.203IVB(I).PMI*J» 

203  FORMAT*  1HA.5X.30HM0  MORE  TRAJECTORIES  FOR  VB».F10.1*2X.7HAND  PH* 
1*F3«2*2X*13KB£L0W  9000  Ml* 

GO  TO  13 

C  CALCULATION  OF  SECOND  ASSIST  TRAJECTORY 
201  VIMF3*SORT(V2»»2+VP2»»2-2.»V2»VP2»SPH!3l 
A4»GCP2/VINF3**2 


02*0. 

310  E4  «<RAOP2*D2>40C0.27  4A4I/A4 
XKU4*ATANt l./SORT ( E4  ••2-1.11 

XN4DEG  -XNU4/R 
!F(VP2-V2*SPHI3)204*209i204 
204  TAWIP4— V2»CPHI3/S0RT1VINF3«42-(V.’»CPHI31*»2J 
PSI4.ATANITAMNP4)-2.«XMU4+PI 


j 


I 


GO  TO  200 

203  URITE-4.207) 

207  FORMATI  1HA.26HVP2“V2*SPHI3  FOR  THIS  CASEI  , 

GO  TO  212 

204  TAMMP4*  V2»CPHI3/SQRTIVINF3»»2-(V2»CPHI3I«»21 
PS I 4»ATAN 1TAHNP4) -2 .•XNU4 

208  PS140G  °PSI4/R 

V4  •SORTIVP2»»2+VIKFJ«»2-2.»VP2»VINF3»COSIPSI4l I 
RP»4*RP2»IV4  M2I/0CS 
IFIRPV4-2. 1219*223, 223 
219  CPHI3>V1NF3*SIN(PSI4I/V4 
SPHI3.SQ«TI1.-CPHI5»»2) 

TAMF3-RPV4»CPHI 3»SPHI 3/ <1 . -RPY4»SPHI 3»»2 1 
iriPSIA I 327*328 *328 

327  F3*(ATANITANF3l4PI)/R 
GO  TO  211 

328  IF (TAMF3 1209*210*210 

209  FS  *IATAN(TANF3II/R 
GO  TO  211 

210  F3  •IA7ANtTANF31-PII/R 

211  A3  •RP2/I2.-RPY4) 

ASAU  «A3  /1 92901000. • 9280.) 

E3  "SORT! < CP.PV4-1.  )»SPHI3)  »»2+CPHI  S»»2I 
OH2  *A3AU  Ml.-ES  I 
IP (PSI4 1329* 330* 530 

329  THT9T*THT12-F34380. 

TIKOP  •TIH124FL0AT1NSP1»TAUP-TIKE(A9*GCS*E3*F3)4TAUP 
GO  TO  923 

980  TMTOT  •THT12-F9 

TIHOP  •TIH124FL0AT(NSP1»TAUP-T|KE(A5*0CS*E9*F3I 
323  TAUPP>2.*Pl«$ORTIA3M3/GCS)/t6400. 

DO  300  HM»1*S 
DO  300  N43*l*9 

IPIAB3 1  FLOAT (NS2l*TAUPP-FL0ATIMP2)4TAUP2)-.2 1  SOI *300*300 
300  CVRTJRUZ 
123 

IPID2-90SO* 1 310*310*202 
CALCULATION  OF  THIRD  ASSIST  TRAJECTORY 
381  V!RF3*«»TiV4*J*4VP2*»2-2.»VA»VP2»5PHI9) 

A4»6C32/VINF3**2 

C4I1I»0. 

■  00  312  K2*1«L2 

C4IK21*IRA0P24D4IK2I*40I0.2T4A»I/A» 

XRU4>ATANt l«/SORTIE4(K2|P*2-l* 1 1 

XHJ0C0(K2)»XKU4/a 

IF«VF2-Y4»*»HI3)306.309*304 
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Program  5  (Contd, ) 

306  TARNP6.-V**CPHI3/S0RTIVINF5»»2-tV*«CPHr5)«»2» 
PSl6*ATANtTANNP6)-2.*XJ!U6+PI 

00  TO  303 
305  KRm<6»*07» 

307  F0«*ATflHA.26HVP2»V*»5PHIS  FOR  THIS  CASE) 

00  TO  312 

SO*  TAWJP6-*  V*»CPHI5/S«mVIKe5*»2-;V*«CPHI3)*»2» 
PSIt*ATAHt7ARKP6)-2.»XHU6 
S0»  PS!60GtK2)*P516/R 

V6<K2>*S»!TtV?2«92+V!NFS»»2-2.«VP2»VtNF5«COSIPS16)J 

Rm»RP2»IVS<K2J*«2)/0CS 

IFIRPV6-2. >319.322.322 

322  l’K2+l 
00  TO  323 

319  CPHI7«.VINF5»SIN<PSI6l/V6tX2) 

SPHI7«SORT(l.-CPHI7**2) 

TARF7«f8>VS*CPH!7»SPHI7/li.-RPV6»SPHI7<H>2> 

IF(PSI6)32*. 330.330 

329  IFtTAKF7IX29.829.309 

<2*  F7<K2>*tAfAN<TANF7)+PI)/R 
GO  TO  311 

330  IFtTARF7l309.325.325 
309  F7<X2I*<ATAH(TA*P7»!/R 

GO  TO  311 

323  F7<K2)*tATAN(TARF7)-PI J/R 

311  A7<K2>-RP2/t2.-RPV6> 

A7AUIK2 >«A7( K2 )/( 92901000. *5280. ) 
E7CK2>«S0RTttlRPV6-l.>*SPHI7l»#2+CPHI7»«2> 
QN*(K2)«A7AUtR2)*(l.-ET(IC2)t 
IFtPSI6l629.630.630 

629  THTMtK2J«THT12"F7tK2>+3*0. 

TIH02tK2)«TIH126Ft.OAT<NSP)*TAUP-TlME(A7tK2>.GCS.E7(K2).F7tK2n*FLO 
1AT  t HS2 ) »TAUP?+TAUPP 
GO  TO  323 

630  THT02-K2)*THT12-F7IK2J 

T!K02<K2)«T!H12+FlGAT<NSP)*TAUP-T!MEtA7(K2).GCS.E7fK2).F7(K2)HFLO 

lATtR33l«TAUPP 
'  323  0*(K2+l)«D*(K2)*100o 

312  CONTINUE 

WRITEt6c213)Vfi( I ).PHf I.JI.O.TAUP2.NP1.NSP.THT12.F1.A1AU.E1.V1.TIN1 
12. T.F3. AA.E3.V2.XNUDEO.PS I DEG. TAUP.QH.02. VI NFS. V6.XN60CG.PSIADG.F5 
2.THT6T.T1N0P.m9AU.E5.TAUPPiMP2.NS2.ON2.(D*(K21  .V2NF3.V6IK2I  .XN6DEG 
3<K2I.PSI60G(K2).F7IK2).THT02tK2).TIH02(R2>.A7AUIK2)iE7tK2).QM*(K2) 
*.K2°L.L21 

213  F0RNATtlHA.5X.3KV8°.F10.1«3X.3HPH».F10.2.5X.6H00CAl«.F8.2.3X.6HTAU 
lP2-.F7.1.9X.*HMPl>.I3t5X.*HNSP-»I32/2Xt*HTH12.3X.2HF1.6X.2HA1.5X.2 


2HE1.6X.2HV1.6X.3HT12.5X.1HT.6X.2HF3.6X.2HA3.3X.2KE3.6X.2KV3.7X.1HN 
3.7X.3HPSI .*X.*HTAUP.*X.2HQH//2F8.2.2F7.*.F10.1.2F7.1.F8«2.2F7.*.F1 
*0.1.2F8.2.F7.1.F8.3///3X.5HOOCA2.5X.3HVINFS.7X.2KV5.8X.2HN*.6X.*HP 
3SI6.6X.2HF3.6X.5HTHT0T.3X.3HTIH0P.5X.2HA5.6X.2HE5.AX.3HTAUPP.2X.3H 
6HP2.1X.3KNS2.3X.3HOHZ^7F9.2.2F21.1.*F9.2.F8.1.2F8.*>F7.1f2I*.F9.5F 
7/3X.5HOOCA3.3X.3HV2NF3.7X.2HV7.8X.2HN6.6X.6HPSI6.6X.2HF7.6X.5HTHTO 
62.3X.3HTIH02.5X.2HA7.6X.2HE7.6X.3HQH3//tF9.0.2Fll.li*F9.2.F8.1.2F6 
9.6.F9.5I) 

DntH23« 

GO  TO  10 
13  CONTINUE 
STOP 


6.6S68E21  1.6077E16  1.1688E16  6.9032CU  3.3680E11 
9.7T31EO*  1.169SE05  2.0902E07  2.0013E07 
3EOF 
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ve» 


♦J444.4  FP“ 


-■> <S76ir 


Program . 5 

Tssxmii.da 


Output 

twm*  u*;s - 


TOP 


«IM  J 


-am 


JU_ 


_A1_ 


_.u_ 


jn _ Til .  ...t_ 


_£1 _  *>_.  .  (> _ V»  -N _ »SJ.  ..TAUP  « 

Hdt  •!»).»!  s.tm  o.apm  niiw.)  42.4  104.2  -144.40  e.Mto  0,5145  98441.4  t.tr  sj. j  -  1*4.7  0.25471 


occa'i 


"v®j 


wT 


“ST 


"P'614“ 


TT 


“YKT07 — 71H0P - *r 


£5  rAOPP“HM:  N$2  W2 


- iTtjres  5BTj7>  iitibTi  T55  4i;a~=Mo7f4“  'nr.vp  '  »i.v  o.'aifsTAor  'np.r"  s  ■  •>  e.woso 

“OOCA3  VjBfi- 


n - tkToj— ypst" 

■nr.nrr«iw 

*9,44  224.03  1475*5 

ratio  224.7V  t£7*V5" 
47.43  223.74  1476.3 


~Vf~ 


TT' 


-atti- 


’OIWO‘1  '67443*4  0Ya*43S”‘ 
0.4440  0.4721  0.14434 
'074'52'4— bT«n—  O'.  15210“ 
0.4564  0.4847  0.13474 
0VVW9“0TWI4~«71WI<-' 

_  _  0.4427  0.6472  0.16324 

«iT»rn-67*70~0^*?S— 676V34“TrnW43“- 
143.19  221.34  1473.9  0.4440  0.4400  0.16743 

W.'B — 121.19 '  1473.0  0.4444  0.434*  011 6937 

144.33  220.00  1473.4  0.4700  0.6343  0.17043 

- T«“T<T*Tr^74y«— 07*3VY-#nT2l4- 

220.42  1473.4  0.4723  0.4331  0.17327 

«07f2  “  Y47S74"  0.“4W  OV4117 

■143.70  220.04  1473.7  .0.4741  0.4303 

T43.34  aV.IO  <673.7  0.4744'  0.4294 

143.40  219.74  1673.7  0.4733  0.4244 

~n*7«—rm.r-07 

210.30  1473.7  0.4747  0.4267 


34.33 

3474T 


rarrrr 

143.13 

TiKW 


0.17427  ' 
0.17316 
0.17596" 
0.174*4 

BTtmrv— 

0.17793 

“0717444" 


-162.93  210.29  1473.7  0.4777  0.4233  0.17494 

'-1£27»'6“  " 219. 20  '1673.7  0V4T(1“'6.*24l'  0.17444 


THU 


_  .oiooc^a.. 

n _ *i 


-tHi 


-00.00- 


j»aii-U4ajSfi _ i*yf2S.  ?j4,5._ 


712 


fi 


43 


£3 


34.33  -123.43  0.6441  0.4614 _ 111300.3 


Ml* 

VI 
00430.4 


HIP- 


TAUP 


OH 


“I0f42'  ' 

“  viWi 

“  H4-  ' 

P$14  P5  *■; ' 

’’THTOT 

TINOP 

~'A5  ' 

'85 

TAUPP  NP2 

NS2  0K2 

~ 210754 

— Tift  177“ 

'  42497.4 

. 9.10 

41.62  -160 .*49  " 

217.04 

1226.7 

0.4*666 ’"0*6100 

126*7  5 

9  0.19042 

oCWi  *" 

'  VtHPS  * 

V7"“* 

N6 

4314  £7 

THT02 

TIH02 

A7‘ 

87 

0H3 

0." 

■  jjsrt.r 

61522.4' 

'9.59 

22.44  -171.44 

227.63 

2351.1 

0.4396 

0.6836 

0.11904 

1000. 

55593.7 

71004.5 

7.64 

24.13  >169.47 

226.02* 

2350.5 

0.6670 

0.4720 

0*16651 

2064. 

-mrsr.y 

“772791.2 

£1)3 

26.92  -166.39 

'224.74 

2130.1 

0.4324 

0.4*13 

0.18212 

3000. 

55593.7 

74127.6 

5.43 

30.76  -147.43 

228.74 

2349.9 

0.6546 

0.6566 

0*15679 

"  4000. 

33393.7 

73137.0 

4.73 

32.13  -164.46 

223.03 

2349.6 

0.6600 

0.4314 

0.160)6 

5000. 

53393.7 

75973.0 

4.22 

33.19  -166^6 

222.43 

2349.7 

0.4*24 

0.4472 

0.14126 

4000.' 

‘  33393.7 

76634.9 

3.79 

14.04  -U8Y60 

221.94 

2349.6 

0.4*80 

0.6630 

0.1*3** 

7000. 

53593.7 

77142.2 

3.44 

34.73  -165.19 

221.54 

2149.8 

0.44*9 

0.4407 

0.147*7 

i443T~ 

32*9172 

“77642.1 

3.14 

11. ll  -164.(4 

221.19 

2349.3 

0.4*46 

o.»«? 

0*14939 

WOO. 

.  33391.7 

760J3.9 

2.91 

13.40  -144.85 

220.69 

2849.3 

0.4700 

0.6366 

0.17007 

i55oo. 

31591.7 

73171.3  '2170 ' 

'  ;«.2l“'-T64V2V 

'220 .’64 

*  2349*4'" 

”0.4712* 

0.6144 

0.1721* 

nooo. 

_ 33591.7 

_ 77663.5 

2.52 

96.56  -164.07 

220.41 

2349.4 

0.6723 

0*6331 

0.17129 

12Qq5T 

ili«7» 

”  2116 

ITTlI-COTlY 

”0.4117' 

"0.17429 

os  ■ 

*  ‘  ’  . 

uoooT 

53393.7 

*  79162*.  3* 

2.22 

37.17  -163.69 

220.04 

2349.4 

0.67'  1 

0*6303' 

'0*17519 

_ IASWL,  3339}, 7  . 

.  .77} 3.6.1 

2.10 

37.42  .  -143.53 

219.66 

2349.4 

0.67<>3 

0*4296 

0.17599 

15000* 

53593.7 

7*536.0 

1.99 

37.0$  -163.39 

219.74 

2349.3 

0.6715 

0.6266 

0*3*671 

1.6000  • 

31393,7  . 

79703.6 

1.69 

37.64  -163.24 

219.41 

2349.3 

0.6762 

0*4275 

0*17756 

17000. 

31391.7 

79153.3 

1.60 

36.02  -.dS.15 

219.49 

2349.3 

0.6767 

0*6247 

0*17794 

leooo. 

.3139}. 7 

79969,4 

1.72 

36.19  -169.04 

219.39 

2349.3 

0.6772 

0*6260 

0.17480 

19000. 

333*3.7 

40111.9 

1.64 

84.14  -Itf.W 

219.29 

2349.3 

0.6777 

0*6253 

0*17900 

_-J«989*- 

...,5^91.7. . 

.60^26,?. 

.  .1*?* 

36.47  -162.65 

219.20 

2349*3 

0.6762 

0*4247 

0*17966 

ve« 

43000.0 

PH* 

90.00 

00CA1*  3415.00 

TAUP2- 

224.5 

API* 

3  M5P*  6 

THU 

PI  A1 

'"61 

VI 

112  T  P3 

A3 

S3 

V3 

N 

PSI 

TAUP  OH 

5*6.35  - 

123.65  0.6941**0.4616 

11.390.3 

42.6  104.1  -119.4*  0.8970  0.8924 

102057.0  5. 

12  62*42 

168*3  0.29704 

00CA2 

VINF3 

V5  . 

H4  .  . 

MI4  •  P5 

. TH707 

TIHOP 

A5 

65 

TAUPP  WP2 

NU  0X2 

141,00. 

53593.7 

65^61.1 

M» . 

44.11  -156.04 

115.16 

777*6 

0.6961 

0*3979 

128.1  4 

7  0*20026 

DOC  A3 

VIHF5 

W7  . 

. . _  NO.  .  . 

PSI4  P7 

THTC2 

TIH02 

AT 

17 

«».  . 

mmrn 

53593.7 

70671,0 

9*59 

.24,.9)  -1T0,15_ 

224.70 

1677*1 

0*6662 

0*4767 

0.14170 

■  l'rl'fl 

53393.7 

72725*9 

7.66 

26. 63  -146*44 

224*79 

1676*5 

0*6522 

0.4484 

0.11210 

■  111  ri 

.  53593.7. 

76616*6. 

6,35 

11.41  -147,94 

223*42 

1675.2 

0*6582 

0*4542 

0.18447 

Bwvfl 

53393.7 

76013*6’ 

5*6) 

11.24  -1*4.0* 

222*40 

1676*1 

0*6629 

0*4470 

0*14340 

60C0* 

33593.7 

77087.2 

6.75 

14.41  -141.24 

221*61 

1675*9 

0*6466 

0*4414 

0.14742 

5000. 

33593.7 

77915.7 

v;a 

11.44  -1WV42 

“220.97 

1675*9 

0.6496" 

■tfi'4844' 

0.17050 

6000* 

33393.7 

76622*6 

3*79 

14.82  -1*4,10 

220.45 

1675*6 

0*6721 

0*6)3) 

0*17)13 

7000. 

53593.7 

79169*2  ' 

)V66 

17.22  -141. Pi 

220*01 

1675*6 

0*4742 

'0.4801 

.:tJii 

4000.. 

.91891.7 

79646*7 

3*16 

37.79  -U3U9 

219.44 

1675*7 

0.47*0 

0*4277 

0.17721 

9000. 

91*91.7 

66069*3 

2*91, 

30*26'  -162*94 

219.82 

1478.7 

0.4778 

0.4288 

0’ 17412 

10000. 

33593.7 

60617*7 

2*70 

36*70.  -162*70 

219.08 

1*75.7 

0,47(9 

0.4234 

0.14028 

uooo. 

91191.7 

80720*6 

2*52 

39*06  -162*46 

214.41 

1475.7 

0*4101 

0.4220 

0*16144 

12000* 

33593*7 

60966*6  .. 

2*34 

39*36.  -162*25 

216*60 

1675*6 

0*4611 

0.4208 

0*16255 

uooo. 

53393.7 

61222*2 

2*22 

39.66  -162*06 

216*61 

1*78.* 

0*48*0 

0.4191 

0.14152 

16000. 

53393*7 

61631*9 

2*10 

19.91  -1*1.49 

216*26 

1*79.* 

0.4429 

0.4141 

0.14440 

18000. 

51593*7 

61419*9 

1*99 

40.11  -141.74 

216*09 

1675*6 

0.413* 

0*4171 

0.11514 

1*000. 

53593*7 

61789*6 

1*6« 

40.11  -141. 40 

217*95 

1475*6 

0.4(41 

0*4141 

0*18589 

iTsoa. 

5)393*7 

61963*1 

1.(9 

40.81  -141.41 

217.(2 

1*78.4 

0.4449 

0.4181 

0*14454 

16000* 

53593*7 

62083*0 

1.72 

40.47  -141.8* 

217.71 

1475*6 

0,4(18 

0.6149 

0*1671) 

19000. 

3!S9J*T 

62210*6 

1*64 

40*82  -161*26 

217*60 

1478.4 

0*4569 

0*6138 

0*16747 

20W. 

11591,7 

62)28*2 

1.17 

40.94  -141.14 

217*51 

1478.4 

0.4I4S 

0*61)2 

0.14417 

HO 

«Mi  TMJECTOftllS  PO*  V6 

•  41000.0  MO  PH*  90.00  4ELQH  8004  HI 
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8  JOS 


0,3,9000 


Prog 

67-012.LT  MYERS.AFIT-SE 


ram 


6 


S1SJOB  MAP 

8I3FTC  MAIM  M94/2.XR7 

C  VENUS  CATEGORY  1  OUT-OF-ECLIPTIC  TRAJECTORIES 

TIME(A.S»EfT|.tE.»ATANISMT(<I.-E)/U.AE>l*TAM(T«R/2.»)-E»S9STa.- 
l£9*2)9StNIT®R)/(l«*£*C0S(T9R> )  )9&0RT( (A**$)/S|/866G0, 
READtS,l>GCS«6CPl»GCP2,RPl.RP2.VPl»VP2.RADPl*RA0P2*RH0 
1  FORMATISEIO.O/SSIO.O) 

DIMENSION  D(2St ,E2< 231, XRU<2S!,V!NFOX(23>,VINFOY12S)»VINFOZ(23l,DX 
1<2S>.DY(2S),DZ<2S>*V2XI2S)«V2Y(2S)*V2Z<2S}.V2I23).PH3(2S)»XI2<23). 
2F3I2S1.A3AU(2S1*E3(2S1.T0P(25!,0M(2S) 

VB-JSOOO. 

PI«3, 14159 
RoPI/180. 

H*33 

N-l 


L-l 

N8«lO 

MO-il 

C  DEPARTURE  (WITS 
DO  2  I»';*M 

VIMFHoSORT(va»«2-2.*GCPl/RADPll 

FH0-90. 

DO  3  J»1»N 
XI-O. 

00  4  K*1.L 
PHOR»PHG»R 
XIR»Xt«R 

!FUVPl«SIM(PH0R)*C0S(XIR»)»»2+ViMFPI»»2-VPI»»2)SjS.6 
3  WR1T£(6,7)V8,PH0,XI 

7  FORMAT (IHA*20X*17H?IO  ESCAPE  FOR  V8*.F10.1»7HAT  PH0*,F7.2,6HAJ10  I*. 
1F6.2I 
GO  TO  31 

i  VO-VPl»SiniPHOR>«COSCXIR>-SQRTI(VPl»S!N(PHoai»eOStXIRIJ«2+V!HFPl» 
1«2-VP1*»2! 

V0X*V09SI<4.rH0K)«C0S<XIRJ 

VOY*VO®COSCPHOP.) 

VOZ*VO»SIS45FHOR)«SINIXIR) 

C  PRE-ASSIST  TRAJtCTORY 

RP1VO«RP1»VO»02/GCS 

TANF0»RP1 VO»COS( PHOR 1 *J IN (PSiOR ) / ( 1 .-RP1 VO«S IN I PHOR ) **2 ) 
IF(TANFO)3J,9,9 
S  FO-ATANCTANFOI/R 
GO  TO  10 

9  FO-(ATANITANFO)-PI)/« 

10  A1-RP1/I2.-RP1V01 


^iisORTl}?8?lVO-lI*J|?8jiHORJ  )  *»24C0i'  I PHOR I  r»2 1 
"PER1»A1«(1.-E1» 

IF(RPER1-RP2)11. 11.12 

12  VAITE<6,13)VB,P(I0«XI 

13  F 0RMAT(1HA,1 OX »3HVB",F10,1#3SH INSUFFICIENT  TO  INTERCEPT  P2  AT  PKO- 
l.F7.2»6HARD  I-.F7.2) 

GO  TO  31 

11  Vl«SORT<OCS*(2./RP2-l./Alll 

PHlR-ARSniSORTIAl«2»(l.-El»»2l/IRP2»(1.9Ai-RP2l>l» 

PH1-PH1R/R 

RP2V1-RP2»V1»»2/GCS 

TANFl«RP2Vl»C0S(PHlR)«SIN(PHlR)/ll«'ftP2Vl*SlNIPHlR|*»2) 

IFiTARFl)l4*14,15 

14  P1»ATAH«TARF1)/R 
GO  TO  16 

13  Fl«f ATANITANFl I-Pll/R 
.16  THT12-F1-F0 
THR«THT12*R 

TIN12«T1ME«A1?GCS.E1*F11-TIME(A1,OCS.E1.FOI 
R1XP— RP2»SINITNRI»C0SIXIfi) 

R1XPAU*R1XP/I92901000.»32e0.) 

R1YP>JIP2*C0SITHR> 

RlYPAU«niYP/S 92901000, *3280. ) 

R12P*RP2»SIN<THRl»SIRIXIRt 
R12PAU«R12P/i 92901000,93280. > 

IF(Pl',’P-XKC'17,Ra,28 

28  «ITt-V»2f)X,.Vi.rrtO 

29  FORMA’  lKA,10X,2‘iI»,F3.2»SOHTOO  LARGE  FOR  INTERCEPT  AT  VB>,F10.1,8 
1HAK0  i  «0-,F8.;») 

00  TO  SO 

C  CALCULATION  OF  E  MATRIX 
17  ETR-1MR-PI/2, 

en»C0SfETR»P:CiCTHS|9C0SIXIRl451METR|93IMtTK*V 

El2r-COSIETR)»SlM(TKR)9COSIXIR)4SiaiETR)9COSITN8* 

E1X«C0SSETRI*SINIXIRI 

E21<-SIHIEVP1»COS(TRA»9COSIX!5!I4COSISTR|9SIN(TKRI 

f42«'SSBI£78)*8INtTHSI#CG5<XIR)9COSIETRt,COSITHR> 

E1S»-SIR(E70>*SIM(XIRI 

m*-C0S<T.'«)9SlHtX!RI 

02.SIN(THAl»SINtXIRl 

R».3"C0Sm«l 

Y  1X““VU-S  IR  t PH1R)  8£  1 1-V19C0S  ( PH1R 1  *E  12 
VlY"-Vi9SlfHPH,R>*E21-V19C0StPHlR)»E22 
V12«TVl9SINIPHlR)*E3l-Vl»C3Stpmft»rIS2 
VINFIX*V1X 
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Program  6-  (Co&tlL  i 

vinpiy-viy+vpz  a 

VINFIZ-V12 

viNFi-MtnviNFiXME+viHFiYMZ+viNFiz--*.) 
VJHFXY*SCflT(VINFIX»»24VfHFIY»»2l 
PS1R-ARC0SI-VINFIY/VINFXYI 
PSI-PSIK/R 
lP!yiZ)lt*lt.l9 
It  ALPR-ARCOSIVINFXY/VINFI! 

00  TO  20 

19  ALPS<s-ARC05<VINFXY/VlNFM 

20  ALP-ALPR/R 
VINF0-VIHF1 
A2-GCP2/VINFI«2 
A2P"A2/RADP2 
VRITEI6.32) 

32  FORMAT  I  lHAtl0X«100H***R«+»*»»*4»s**»*»»®4,»***»»A***»P*P*,»,,,,,»»#» 
lUtllCIIHMIlHIKfiaHMMMOtllllllfMIKIOmtl  • 

HR  I TE (6*21 )V8#PH0»X J»VINFPl.V6fVOX#VOY#VOZ|VllVlXiVlY#VlZ»THYlS*FO 
l.AlAU.El.PHI.RIXPAUfRIYPAU.RIZPAUtTIHIZ tVINFI fVINFXYsPSI *ALP 

21  FORMAT(lHA.49A»3hVB-.F10.l«3X.4HPHO-»F6.1///SX»im.9X.6HVtNFPl,4X. 
12HV0#8X»3HV0Xf7X.3HV0Y#7Xf  3KV0Ztl7X*2HVl*8X*3HVin*7X*3HVlY.7X»lJ(Vl 
22//F7,3.9F10.1.10X.4F10.1///;,X*4HTH12»3X<2HFO!6X«2HA1.4X»2H£l»j!I8.3 
3HPHl»4X*4HRlXP.4Xt4HRlYP#4X*4HR12Pt4X.3HTIN12»4X*9HVlttPI#4X.6HVIXF 
4XY.9X.3MPSI.3X.3HALP//2F8.2.ZF8.4.F8.2.3Fa.4.FJ.2,2F10.1.2F8.2) 

BEt-0. 

00  22  IB-1. M3  ’ 

BETR-BET»R 

C  CALCULATION  OF  f  MATRIX 

F)l-C0S(BETfi)«C0StPSIR)-SIN(BETR)»SINIALPR|»SINlP5IR> 

Fl2o-COS(ALPRI«S1NIPSIRS 

F13«SIN(BETR)*C0SIPSIR>+C0S(BETR)»SINIALPR|»5IH(PSIR) 

F21-COS(8ETR)*SiN(PSIR)+SIN(BETR)iSINIALPR)»COSCPSIR> 

F22-COS(ALPRI«COSIP$IR) 

F23-SINIBETR)*SIN<PSIR)-C0SIBETRI»SINIALPR)»C0SIPS!R> 

F3l— S1NIBETR)*C0SCALPR) 

F32-SINIALPR) 

F33-COSIBETR)»COS<ALPR) 

C  HYPERBOLIC  PASS  TRAJECTORY 

0(11-0. 

00  23  I0-1.M0 

£2<ID)-(RAOP2+0(ID)»6080.27+A2)/A2 
XNUR-ARSINI 1«/E2< 101 ) 

XNUt I0)aXNUR/R 

VNF0X-VINFO«C0St  P I/2.-2.-XNUR ) 

VNFOY— VIXFO»SlN(PI/2.-2.«XNURI 
VNFOZ-O. 


VINF0X<lpi-VNF0X»Fll+VHF0Y«F12+VNF02»F13 
VIHFOY(IOI-VNFOX».-21+VNFOY*F22+VNFOZ«F23 
VINFOZ(IOI-VNFOX«F31+VNKOY»F32+VNFOZ»F33 
02X-IA2+DI ID)f6080.27+RADP2)*COS(XNUR)/RAOP2 
O2Y«(A2+0(ID)«6080.27+RADP2)»SIN(XNUR|/RA0P2 
022-0. 

DX( IDI-02  •FllK>2Y»F12tD2Z«F13 
0Y(I0)-02X»F21*D2Y»F22*02Z»F23 
0Z( IDI-02X«F31+02Y»F32*D2Z»F33 
V2X( I0I-VIKFOXI IOJ 
V2Y(IO)-VINFOYIIO)-VP2 
V2Z(IOI-VIHFOZ(ID) 

V2(IO)-SORT(V2X(IOI*»2+V2Y(IO)«»2*V2ZIIDl»«2) 

PH3R-ARCOSIV2XI IOI/V2I 10) 1 
PH3( 10) -PH3R/R 

XI2(ID)-ATANtV2Z<IDI/V2Y(IO))/R 
RP2V2-RP2»V2t ID)»»2/GCS 

TANF3-RP2V2»C0S<PH3R)»5!N(PH}?)/U.-RP2V2'  SIHlPHJr,  »»2> 
IFITANF3I24.29.29 
24  F3(10I-(ATAN(TANF3) l/R 
GO  TO  26 

29  F3(I0)-tATANnANF3l-PI)'R 

26  A3-RP2-' '.-RP2V2) 

A3AUIU  A3/ 1 97901000  **3760.) 

F3I IDI-SORTI C  (NP2V2-1.)  •>4INIT-H3R)  l••^♦C0S<PH3RI*»^l 
TOPI IOI-TIH12-T IM£(A3»GCS.E3I 10) .F3I 101) 
OM(ID|-A3AU<IOI*Cl.-E3<IO>l 
0(10411-01  ID) +  10000. 

23  CONTINUE 

WRlTE(t*27)6ET.A2P.(D(lD).XNU(lD).V2( IDI .V2XI IDI.V2YI 10) .V22I ID) *0 
lX(IO).OY(iO).OZ<IO).TOPUOI«PH3(lC)>F3UO)»A3AU(IOI.E3(ID)»XI2(!0) 
2.0M1 10) » 10-1.N0) 

27  FORMAT  1 1HA.49X.9HBETA-.F7.2.3X.3HA2- .F7.2///4X. 4HD0CA»4X»2HNUt6Xf2 
lHV3»6X«3HV9*»6T«3HV3Y»6X»3HV3Z»?Xt2H0Xs4X«2HDY»4X»2K02t4X»3HT0P.4X 
2.3HPH3.3Xt2HF3,6X.2HA3.3X.2HE3*3X.2Hi3.5X»2HOH//<F9.1.F7.2i4F9.1.3 
IF6.2.F7.1.F7.2.F8.2.2F7.4.F7.3*F8.9)> 

BET-BET+10. 

22  CONI INUE 
XI-XI+.l 
4  CONTINUE 

30  PHO-PHO-2. 

3  CONTINUE 

31  VB-VB+1000. 

2  CONTINUE 

STOP 

END 


4.6864E21  1.4077E16  1.1488EI6  4.9092EU  3.94B0E11 
9.7791E04  1.1493E09  2.0902E07  2.001SC07  2.0200E09 
tEOF 
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Program  7 

8J0S  0*2.7000  47-012.LT  HYERS.AFIT-SE 

SI8J08  HAP 

8I8FTC  WAIN  H94/2.XR7 

C  VENUS  CATEGORY  2  OUT-OF-£CHPTIC  TRAJECTORIES 

TIN£<A.G.£.TI-(2.*ATAKIS0RT{<l.-£)/(l.*£ll«TAN(T*R/2.M-E*$CAT  1.- 
lE**2l*SIN(T«R)/(l.*e*COS(T*R>))*SORT((A**9)/O)/8440O. 

REAP! 5*1 IGCS.CCPl.GCP2.RPl.RP2.VPl.VP2.RAOPl.RAOP2.RHO 
1  FORNATI5E10.0/5E10.0) 

DIMENSION  0(251 .E2I25).XNU(25) .VIKF0XI25I .VIKFOYI2S) .V1NF02I25) .OX 
1(251.07(251.02(25} >V2Xt 25 > .V27( 25) .V22( 25) *V2(25) «PH3(23) *XI2(25) t 
2F5(25).A5AU(25).E3(25).TOP(25).ON(2S) 

vb-ssooo. 

PI-3. 14159 
R-»l/180. 

M-33 
R*t 
HB*10 
MO- II 

R12-RP1/RP2 

A*(l.-R12)/(1»+R12) 

C  DEPARTURE  ORBITS 

00  2  t-l.M 

V INFP 1-SORT (VB*»2-2 .*GCP 1 /RAOP 1 ) 

PH0-90. 

00  S  J-l.H 
PHOR-PHO»R 

El«SQRT(A*»2-Ml.-StN(PH0R>»*2)»tl.*2.*A))/SlM(PH0R> 

40  Al-RPl*StN(PH0R)«(SlN(PH0R)-S0RT(£l*»2+SIN(PH0R)**2-l.))/(l.-El**2 
1) 

VO-SORT (CCS* ( 2./RP1-1 ./Al ) ) 

VEL-(VO**2+VPl««2-VINFPl**2>/(2.*VPl*VO*SiN(PHOR)) 

If(A6S(VELI-l. 142.42. 30 
42  A1AU-A1/ (92901000**5280. ) 

XIR-ARCOSMVO»*2*VP1»»2-V!NFP1*»2)/(2«*VP1»VO*SIN(PHOR))) 

Xt-XIR/R 

102  VOX-VO*SINIPHOR)*COS(XIR) 

VOY-VO»COS(PHOR) 

V02-VO*SIN(PH0R)*SIN(XIR) 

C  PRE-ASSIST  TRAJECTORY 

RP1V0-RP1*V0**2/GCS 

TANFO-RP1VO*COS(PHOR)*SIN(PHOR>/(1.-RP1VO*S1N(PHOR>**2) 

IF(TANF0)8.9.9 

8  FO-ATAN(TAAFO) /R 
GO  TO  11 

9  F0*(ATAN(TANF0)-P1)/R 

11  Vl-SORT  <GCS*( 2./RP2-1 ./Al I ) 


PH1R-ARS1N(S0RT(A1**2*(1.-E1**2)/(RP2*(2.*A1-RP2)))I 

PH1-PH1R/R 

RP2V1-RP2*V1**2/GCS 

TANF1-RP2VI*C0S(PH1R)*S:H(PH1R1/(1.-RP2V1»SIN(PH1R)**2) 
tFITANFl 114.14.15 

14  F1-ATAN(TANF1)/R 
00  TO  14 

15  F1-(ATANITANF1)-PI)/R 
14  FI— FI 

PH1R-PI-PH1R 
PH1-PM1R/R 
THT12-F1-F0 
34  THR-THT12*R 

TIM12-TIHE(A1.GCS.E1.F1)-TIME( A1.GCS.F1.F0I 
R1XP— RP2*SIN(THR)*C0S(X1R) 

R1XPAU-R1XP/(92901000.*9260« ) 

RIYP-RP2*C0S(THR) 

RlYPAU-RlYP/( 92901000. *9260.) 

R12P-RP2*SIN(THR)*SIN(XIR) 

R12PAU-R12P/( 92901000. *5280.) 

IF (R12P-RH0) 17.28.28 

28  VRI7EI4.291X1 .V8.PH0 

29  FORHAT(lHA.10X.2Hl-.r3.2.30H700  LARGE  FOR  INTERCEPT  AT  VB-.F10.1.8 
1HAN0  PH0-.F8.2) 

CO  TO  30 

C  CALCULATION  OF  E  MATRIX 
17  ETR-THR-PI/2. 

E11-COS(ETRI*COS(THR)*COS1X1R)4SIN(ETR)*SIN(THR| 

E12— C0S(ETR)*SIN(7HR)*C0S(XIR)+SIN(ETR)*C0S(THR) 
E13-C0S(ETR)*SIN(XIR) 

E21— SIN(ETR|»C0S(THR)«C0S(XIR)*C03(ETR)*SIH(THR) 
E22-SIN(ETR)*SIN(THR)*C0S(XIR).C0S(ETR)*C0S(THR) 

E23— SIN(ETR)*SIN(XIR> 

E31— C0S(THR)*SIN(X1R) 

E32-SIN(THR)*SIN(XIR) 

e33-C0SIXIR) 

V1X— VIASIPIPHIRI^EII-VHCOSIPHIRIPEIE 

ViY— 51*SIN(PHIRI*E21-V1«C0S(PHIR)«E22 

VIE— Vt«SIN(PHlRI*E31-Vl«COSIPHlR)*E32 

VINFIX-V1X 

VIHFIY-V1Y*VP2 

VINFIE-V12 

VIHFI-SOR7(VIHFIX**2*VIHFIY**2*VINFl'2**2) 

VIHFXY-SORT(VIMFIX»*2*VIKFIY**2) 

IFIVtNFIYI800.800.801 
800  PSIR-O. 
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psI>0>  Program  7  (Contd,) 

CO  T0*302 

301  PSIR*P1 
PSt-180. 

302  IFIV1Z)18.18.19 

18  ALPR«ARCOSIVlNFXY/VIRFtl 
GO  TO  20 

19  ALPR*-ARCOSIVINFXY/V'RFn 

20  ALP-ALPR/R 
VINFO«VINF! 

A2*GCP2/VJRFt»*2 

A2P"A2/RA0P2 

NRITEI6.32) 

32  FORMAT!  1KA.10X.  . 

KR!T£I6.Z1)V8.PK0.XI.V1NFP1.V0.V0X.V0Y.VCZ.V1 IV1X.V1Y.V1Z.TKT12.F0 
1.A1AU.E1.PH1.R1XPAU.R1YPAU.R1ZPAU.TIM12.VINFI.VIHFXY.PS1.ALP 

21  FOfiKATIlHA,*SX.3HVB-.F10.1.3X.*KPH0».F6.l///3X.lHI.SP4»6HVlRFPl,6X. 
12HVO.8X.3HVOX.7X.3HVOY.TX.3HVOZ.17X.2HVt.8X.3KVlX.TX.3HVlY.7X.lHVt 
2Z//F7.3.3F10.1.10X.4P10.1///2X.4HTH12.3X.2HF0.6X.2HA1.6X.2HE1.6X.S 
3MPH1  .4X.4HR1XP.4X.4HR1YP.4X.4HR1ZP.4X.3HT1K12.4X.SHV1RF! .4X.6HVIRF 
*XY.5X,3HPSI,SX.3HALP//2F8.2»2F8.4.ifB.2.3Fa.4.FS.2.2F10.1.2PS.2l 

8ET-0. 

DO  22  I8»1.M8 
BETR-BETAR 

C  CALCULATION  OF  F  MATRIX 

Fll.COS16ETR)»COStP5!R>-SINt0£TRI«SlNIALPR)»SIHIPSlR> 

F12«-C0S(ALPR)®SIHIPS1R) 

F13"SINI8ETR>»C0S!PSIR)+C0S{SETR)«S1N1ALPR)»SI«IPSIR) 

F21.COSI8ETR)»SIMIPS!R)+S!N{6ETR)«SINIALPR)*COStPStR> 

F22«COSIALPR1»COSIPSIR) 

F23.StNI8£TR)»SINtPSlltl-C0SI8ETRI«S!NIALPR>*C0SIPS|R) 

F31—  SINI0ETRI«C0S(ALPRI 
F32«SIH(ALPR1 
F33»C0SIBETR)*C.0S(ALPRI 
C  HYPERBOLIC  PASS  TRAJECTORY 
0(1)«0. 

DO  23  IO«l.MD 

E2 1 101- IRA0P2+DI 101 *6080.27+A2)/A2 
XNUR»ARSIMtl./E21ID>) 

XNU(SD)*XHUR/R 

VNFOX°V!NFO«COSIPI/2.-2.*XNUR) 

V'4FOY*-VINFO»SINIPI/2.-2.»XNUR) 

VHFOZ-O. 

VINF0XIIDI-VNF0X«FU*VKF0Y»F12*VNF0Z»F11 

VIHFOYIIO)«VMFOX«F21+VNFOYYF22*VNFOZ»F21 


VINFOZI 101 ■VNFOX®F3l+VMFOY»F32+VHFOZ*Pll 
D2X-(A2+OIIDl»6080.27»RA0P2)«£OSIXRUR)/RA0P2 
D2Y- I A2ADI 101 A6080.27+RADP2 I *SI Ml XRUR ) /RA0P2 
D2Z-0. 

0X1 lD|.02X»Fll+D2YOr!2A02Z»F11 
DYIID)-D2X«F21+02f«F22+D2Z«F23 
OZI I01.02X»F31+02Y»F32A02Z»F33 
VZXItOl-VINFOXtIO) 

V2YIIOI«VlNFOYIIDI-VP2 
V2ZI 10)*V1NF0ZI 10) 

V2IIO)«SQRTtV2X!!D>«»2»VZYIlD)»»2*V2ZIIOI»»21 

PH3R»ARC0S( V2X I  ID) /V21 101 1 

PH1II0)*PH3R/R 

X12I I0)*ATAMIV2ZI 1D1/V2YI 10) )/R 
RP2V2«RP2»Y2 1 101  *#2/<iCS 

TAHF3.RP2V2»C0S1PH3R)«S1MIPH3R)/!1.-RP2V2»SIM(PH3R1»*2I 
1FITANF1I24, 23.23 
24  FltlDl.lATANITANFllt/R 
GO  TO  26 

23  FI  1 1 01  MATAR(TANFl) -PI  )/R 

26  A3«RP2/*2.-RP2V2) 

AlAUIIDI*Al/t9290>C«5s*S280.» 
ElODl*S0RTIIiRP2V2-l.)»$IH(PHlRn»*2AC0St»MW)M2l 
TOPI I01-T1M12-T1MEIA1.GCS. Ell  101 .F3IID1) 

CHI  10) "A3  Mil  101*1 1.-E3I 101 1 
01 10*11*01 1DI+100. 

21  CONTINUE 

MR  I TE 16 .27 1  SET »A2P. 101 lOl.XXOI 101.V21 10) *V28I ID) .V2YI I0I.V2ZI 101. 0 
1X110) »0Y( I01.0ZI 10 1 .TOPI  101 .PH11I0I » FI 1 10) .A1AUI IDI.CIIIOI.XIZIID) 
2. WHO). 10-1. MO) 

27  FORMAmHA.ASX.3M9ETA.,F?.2»iX.SHA2«.FT.2//m.*HOCCA»6X»IHfK).TX.2 
lHV1.7X.lHVlX,7X.lHVlY.7X.«!VSZ.3X.2HOX.AX»2KSY.»X»2)S)Zs*X.5HTeP.*X 
2<lHPHl)3X.2HF3.6X<2KA}t3X.2HE3i6X.2HI3.3X»2HCM//lFt.l.F7.2.6FXS.l. 
13F6.2.F7.1.F7.2.F8.2.2F7.A.F6.6.F3.31) 

.8ET*BET*10. 

22  COMTIRUE 

10  PH0*PH0-2. 

1  CONTINUE 

11  V8*VS*1000. 

2  CONTINUE 
STOP 
END 

IOATA 

6.6S6BE21  1.A077E16  1.16I1E16  4.903ZE11  1.S420E11 
9.7731E0*  1.1491E03  2.0902E07  2.0013C07  2.0Z0OS89 
1E0F 
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Program  8 

*JOB  0(4*9000  67-012. LT  NYER5.AFIT-M 

3I6JOB  MAP 

3IBFTC  MAIM  M94/2.XR7 

C  LAUNCH  OPPORTUNITIES  FOR  VENUS-ASSISTED  TRAJECTORIES 

TIM£!A»G.E.TI-<2.»A7AN(S0RT!ll.-£l/(1.4Ell4TANtT»«/2.ll-E»SQRTtl.- 

lE«»2l»SINIT»R)/<l.+£»C05tT»R}|l»SORTHA»«3>/8|/»6400. 

READ!  5*1  IGCS*GCPl»CCP2*RPl*RP2(V!,X*VP2tRA9Pl*RAOP2iRHO 
1  FORMAT (5E10.0/5E10. 01 
READI5.601W1.W2.XL10.XL20.S.T0 
60  FORMAT (4F12* 6 *F10*2  »FI2e2) 

DIMENSION  TL! 30) 

V8-38000. 

PI-3.14159 

R-PI/180. 

M»33 

N«4 

1*4 

ITT-13 

C  DEPARTURE  0R8ITS 
00  2  t-l.M 

VINFP1-50RT(VB»»2-2.4GCP1/RA0P11 

PH0-90. 

WRITE16.32I 

32  FORMAT!  1HA*10X«100H**4***.******4>4***4»*44444****4***4*44*4*4»4*4 

DO  3  J-liN 
XI-O. 

DO  4  K-l.L 
PH0R-PH04R 
XIR-XI»R 

IFIIVP1»SINIPHOR»»COSIX1RII»»2*VINFP1»42-VP14»2)».5.6 

5  WRITE! 6*7) VB»PH0»XI 

7  FORMAT! 1HA.20X.17HN0  ESCAPE  FOR  VB-.F10.1.7HAT  PH0-.F7.2.6HAND  I-. 
1F6.2) 

GO  TO  31 

6  VO«VPl*SIN(PHOR)*COSIXIR1-SORTi(VPl*S!N<PHOR)*COSIXIRl|M24VINFP14 
1»2-VP1»»2I 

C  PRE-ASSIST  TRAJECTORY 

RP1V0-RP1*V0**2/GCS 

TAMFO-RPlVO»COStPHORl»SIN!PHOftim.-RPlVO»SIN!PKORI**2) 

IFITANF018.9.9 

8  F0-ATAN1TANF0I/R 
GO  TO  10 

9  FO-tATANITANFOl-PIl/R 
10  Al«RPl/!2.-RPlV0) 

E1-SORT1IIRP1VO-1.1«SIN1PHORI14424COSIPHOR|«»21 


RPER1-A1*! 1.-E1I 
IFIRPERl-RP2lll.il. 12 

12  WR1TE16.131VB.PH0.XI 

13  F0RMATUHA.10X.3HVB-.F10.1.36H1NSUFFICIENT  TO  INTERCEPT  P2  AT  PKO- 
1.FT.2.6HAND  1-.F7.2) 

GO  TO  31 

11  Vl-SQRT!GCS*t2./RPZ-l./Alll 

PH1R-ARSINIS0RT1A1»»2»I1.-E1M2I/IRP2»(2.»A1-RP21I11 

PH1-PH1R/R 

RP2V1-RP2»V1»»2/GCS 

TANF1-RP2V1«C0S!PH1R|4SIM1PH1R1/I1.-RP2V1»S!N1PH1R|*»2> 

IFITANF1 114.14.15 

14  F1*ATAM!TANF11/R 
GO  TO  16 

15  Fl-IATANITAHFlI-PIl/R 

16  THT12-F1-F0 
THR-THT12*R 

T IM12-T IMEIA1.GCS.E1.F1 I -TIMEtAl.GCS.El.FO I 

R1ZP>RP24SIN!THR1«SINIX1RI 

IFIR12P-RH01 17.21.28 

28  WRITE-6. 29IXI .VB.PHO 

29  FOfiMATHHN.;?X.2HI-.F3.2.SOHTOO  LARGE  FOR  INTERCEPT  AT  V6-.F10.1.8 
1HAND  PH0°«F8.*< 

GO  TO  30 

17  THL-THT12-W2«TIH12 
THP-THT12-W1»TIM12 

TL!l)»IXL10-XL204(W2-Wll»T04THL)/!W2-iai 

DO  50  IT-2. ITT 

TCI ITl-TLI 114FL0ATI IT-II^S 

50  CONTINUE 

WRITEI6.511VB.PH0.XI.TIM12.THT12.THL.TKP.ITLII7I.IT-1.ITT1 

51  FORMAT! lHA.48X«SHV8-.FlO.i»2X»4HPH0-*f 7*2 .2X.2HI-.F6.2//S9X.6HTIN1 
12-.F9.2.2X.6HTHT12-.F7.2.2X.4HTHL-.F7.2.2X.4HTHP-.F7.2//49X.12HJUL 
2 IAN  DATES. 10X. 14HCALEN0AR  DATES//IF60.2) I 

XI-X1+.1 
4  CONTINUE 

30  PH0-PM0-2. 

3  CONTINUE 

31  VB-VB41000. 

2  CONTINUE 

STOP 
•  END 
8  DATA 

4.6868E21  1.4077E16  1.1488E16  4.9052E11  3.5480E11 
9.TT51E04  1.1493E05  2.0901E07  2.001SE07  2.0200E09 

.913609  1.602131  100.158150  174.29431  513.92  2436915.00 
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;  Prograln  8  Output 

<  VB*  49000.0  PH0»  90.00  !»  0. 

i 

l  TIM12*  56.52  THT\2=  46.21  THL*  -44.33  THP»  -9.49 

* 

i  JULIAN  OATES  CALENDAR  OATES 

* 

2436742. T 3 

i  2437226.69 

t  .  2437910.59  See  Ref  1 

2438494.53 

2439078.44 

2439662.38 

2440246.28 

*  2440830.22 

’  2441414.13 

2441998.03 

2442581.97 

}  2443165.88 

2443749. 81 
2444333.7? 

2444917.66 
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13.  ABSTRACT 


An  investigation  of  direct-transfer  trajectories  for  solar  probe 
missions  indicates  that  expensive  launch  vehicles  such  as  the  $125- 
million  Saturn  V  are  required  to  achieve  perihelia  less  than  0,25  AU 
and  inclination  angles  above  20°,  Methods  are  developed  for  two-  and 
three-dimensional  gravity-assisted  trajectories  (trajectories  which 
pass  through  one  or  more  planetary  gravitational  fields) ,  and  are 
applied  to  solar  probe  missions  in  an  attempt  to  reduce  launch  vehicle 
costs.  The  analysis  is  based  on  the  pieced-conic  approximation  and 
the  assumption  of  circular  coplanar  planetary  orbits.  It  is  found 
that  a  perihelion  of  0.16  AU  can  be  obtained  with  the  $16-million 
Atlas/Centaur/TE-364-3  by  using  a  Venus  assist.  Even  greater  reduc¬ 
tions  in  perihelia  are  attained  with  multiple  passes  at  Venus;  in 
addition,  these  missions  allow  the  exploration  of  several  regions  near 
the  Sun  with  a  single  launch.  Venus-Mercury  combination-assisted 
trajectories  are  of  little  value  in  solar  probe  missions.  A  Venus 
assist  is  useful  for  10°-out-of-ecliptic  trajectories  when  the  p.2- 
assist  orbit  lies  in  the  ecliptic.  For  a  pre-assist  orbit  out  of  the 
ecliptic,  an  inclination  angle  of  25°  and  a  perihelion  of  0,53  AU  can 
be  achieved  with  the  Atlas/Centaur/TE-364-3  launch  vehicle.  Solar 
impact  and  90°-out-of-ecliptic  trajectories  can  be  attained  with  a 
Jupiter  assist;  however,  traversal  of  the  asteroid  belt  and  a  three- 
year  mission  time  decreases  spacecraft  reliability. 
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